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Abstract 

We review different dark energy cosmologies. In particular, we present the ACDM cosmology, 
Little Rip and Pseudo-Rip universes, the phantom and quintessence cosmologies with Type I, II, 
III and IV finite-time future singularities and non-singular dark energy universes. In the first part, 
we explain the ACDM model and well-established observational tests which constrain the current 
cosmic acceleration. After that, we investigate the dark fluid universe where a fluid has quite general 
equation of state (EoS) [including inhomogeneous or imperfect EoS]. All the above dark energy 
cosmologies for different fluids are explicitly realized, and their properties are also explored. It is 
shown that all the above dark energy universes may mimic the ACDM model currently, consistent 
with the recent observational data. Furthermore, special attention is paid to the equivalence 
of different dark energy models. We consider single and multiple scalar field theories, tachyon 
scalar theory and holographic dark energy as models for current acceleration with the features of 
quintessence/phantom cosmology, and demonstrate their equivalence to the corresponding fluid 
descriptions. In the second part, we study another equivalent class of dark energy models which 
includes F(R) gravity as well as F(R) Hofava-Lifshitz gravity and the teleparallel f(T) gravity. 
The cosmology of such models representing the ACDM-like universe or the accelerating expansion 
with the quintessence/phantom nature is described. Finally, we approach the problem of testing 
dark energy and alternative gravity models to general relativity by cosmography. We show that 
degeneration among parameters can be removed by accurate data analysis of large data samples 
and also present the examples. 

PACS numbers: 04.50.Kd, 95.36,+x, 98.80.-k 
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I. INTRODUCTION 



Cosmic observations from Supernovae la (SNe la) 111, cosmic microwave background 

nn n n 

(CMB) radiation |2H4(, large scale structure (LSS) [5j, baryon acoustic oscillations (BAO) |6[, 
and weak lensing |7| have implied that the expansion of the universe is accelerating at the 
present stage. Approaches to account for the late time cosmic acceleration fall into two 
representative categories: One is to introduce "dark energy" in the right-hand side of the 
Einstein equation in the framework of general relativity (for recent reviews on dark energy, 
see The other is to modify the left-hand side of the Einstein equation, called as a 

modified gravitational theory, e.g., F(R) gravity (for recent reviews, see 12H17]). 

The various cosmological observational data supports the A cold dark matter (ACDM) 
model, in which the cosmological constant A plays a role of dark energy in general relativity. 
At the current stage, the ACDM model is considered to be a standard cosmological model. 
However, the theoretical origin of the cosmological constant A has not been understood 
yet fis| . A number of models for dark energy to explain the late-time cosmic acceleration 
without the cosmological constant has been proposed. For example, a canonical scalar field, 
so-called quintessence 19[, a non-canonical scalar field such as phantom 20(, tachyon scalar 
field motivated by string theories 2l| . and a fluid with a special equation of state (EoS) called 



as Chaplygin gas 



22j-[24(. There also exists a proposal of holographic dark energy 25N27|. 



One of the most important quantity to describe the features of dark energy models 
is the equation of state (EoS) wde, which is the ratio of the pressure P to the energy 
density pde of dark energy, defined as u^e = -Pde/pde- We suppose that in the background 
level, the universe is homogeneous and isotropic and hence assume the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) space-time. There are two ways to describe dark energy models. 



One is a fluid description 28M30| and the other is to describe the action of a scalar field 



theory. In the former fluid description, we express the pressure as a function of p (in 
more general, and other background quantities such as the Hubble parameter H). On the 
other hand, in the latter scalar field theory we derive the expressions of the energy density 
and pressure of the scalar field from the action. In both descriptions, we can write the 
gravitational field equations, so that we can describe various cosmologies, e.g., the ACDM 
model, in which wde is a constant and exactly equal to —1, quintessence model, where w~de 
is a dynamical quantity and —1 < we>e < — 1/3, and phantom model, where Wde also varies 
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in time and u>de < — 1- This means that one cosmology may be described equivalently by 
different model descriptions. 

In this review, we explicitly show that one cosmology can be described by not only a fluid 
description, but also by the description of a scalar field theory. In other words, the main 
subject of this work is to demonstrate that one dark energy model may be expressed as the 
other dark energy models, so that such a resultant unified picture of dark energy models 
could be applied to any specific cosmology. 

This review consists of two parts. In the first part, various dark energy models in the 
framework of general relativity are presented. First, we introduce the ACDM model and the 
recent cosmological observations. At the current stage, the ACDM model is consistent with 
the observational data. We then explain a fluid description of dark energy and the action 
representing a scalar field theory. In both descriptions of a fluid and a scalar field theory, 
we reconstruct representative cosmologies such as the ACDM, quintessence and phantom 
models. Through these procedures, we show the equivalence between a fluid description and 
a scalar field theory. We also consider a tachyon scalar field theory. Furthermore, we extend 
the investigations to multiple scalar field theories. In addition, we explore holographic 
dark energy scenarios. On the other hand, in the second part, modified gravity models, 
in particular, F(R) gravity as well as F(R) Hofava-Lifshitz gravity and f(T) gravity with 
T being a torsion scalar 3lM33j. i.e., pictures of geometrical dark energy, are given. It is 
illustrated that by making a conformal transformation, an F(R) theory in the Jordan frame 
can be moved to a corresponding scalar field theory in the Einstein frame. It is also important 
to remark that as another modified gravitational theory to account for dark energy and the 
late-time cosmic acceleration, F(R, T) theory has been proposed in Ref. 34|, where T is the 
trace of the stress-energy tensor. We use units of fee = c = K — 1 and denote the gravitational 
constant 8?rG by k 2 = 8vr/M P1 2 with the Planck mass of M P1 = G~ 1/2 = 1.2 x 10 19 GeV. 
Throughout this paper, the subscriptions "DE", "m", and "r" represent the quantities of 
dark energy, non-relativistic matter (i.e., cold dark matter and baryons), and relativistic 
matter (e.g., radiation and neutrinos) respectively. 

The review is organized as follows. In the first part, in Sec. II we explain the ACDM 
model. We also present the recent cosmological observational data, in particular, in terms 
of SNe la, BAO and CMB radiation, by defining the related cosmological quantities. These 
data are consistent with the ACDM model. 
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In Sec. Ill, we investigate a description of dark fluid universe. We represent basic for- 
mulations for the EoS of dark energy. We also introduce the four types of the finite-time 
future singularities as well as the energy conditions and give examples of fluid descriptions 
for the ACDM model, the GCG model and a model of coupled dark energy with dark matter. 
Next, we explore various phantom cosmologies such as a coupled phantom scenario, Little 
Rip scenario and Pseudo-Rip model. Furthermore, we show that the fluid description of 
the EoS of dark energy can yield all the four types of the finite-time future singularities. A 
fluid description with realizing asymptotically de Sitter phantom universe is also examined. 
In addition, we investigate the inhomogeneous (imperfect) dark fluid universe. We study 
the inhomogeneous EoS of dark energy and its cosmological effects on the structure of the 
finite-time future singularities. Moreover, its generalization of the implicit inhomogeneous 
EoS is presented. 

In Sec. IV, we explore scalar field theories in general relativity. We explicitly demonstrate 
the equivalence of fluid descriptions to scalar field theories. In particular, we concretely 
reconstruct scalar field theories describing the ACDM model, the quintessence cosmology, 
the phantom cosmology and a unified scenario of inflation and late-time cosmic acceleration. 
We also consider scalar field models with realizing the crossing the phantom divide and its 
stability problem. 

In Sec. V, we examine a tachyon scalar field theory. We explain the origin, the model 
action and its stability conditions. 

In Sec. VI, we describe multiple scalar field theories. First, we examine two scalar field 
theories. We investigate the standard type of two scalar field theories and the stability of 
the system. We then introduce a new type of two scalar field theories, in which the crossing 
of the phantom divide can happen, and also explore its stability conditions. Next, we extend 
the considerations for two scalar field theories to multiple scalar field theories which consist 
of more scalars. and clearly illustrate those equivalence to fluid descriptions. 

In Sec. VII, we study holographic dark energy. We explain a model of holographic dark 
energy as well as its generalized scenario. In addition, we examine the Hubble entropy in 
the holographic principle. 

In the second part, in Sec. VIII, we consider accelerating cosmology in F(R) gravity. 
First, by using a conformal transformation, we investigate the relations between a scalar 
field theory in the Einstein frame and an F(R) theory in the Jordan frame. Next, we 
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explore the reconstruction method of F(R) gravity. We explicitly reconstruct the forms of 
F(R) with realizing the ACDM, quintessence and phantom cosmologies. In addition, we 
study dark energy cosmology in F(R) Hofava-Lifshitz gravity. We first present the model 
action and then reconstruct the F(R) forms with performing the ACDM model and the 
phantom cosmology. Furthermore, we explain F(R, T) gravity. 

In Sec. IX, we describe f(T) gravity. To begin with, we give fundamental formalism and 
basic equations. We reconstruct a form of f(T) in which the finite-time future singularities 
can occur. We also discuss the removal way of those singularities. Furthermore, we represent 
the reconstructed f(T) models in which inflation in the early universe, the ACDM model, 
the Little Rip scenario and the Pseudo-Rip cosmology are realized. In addition, as one of 
the most important theoretical touch stones to examine whether f(T) gravity can be an 
alternative gravitational theory to general relativity, we explore thermodynamics in f(T) 
gravity. We show the first law of thermodynamics and then discuss the second law of 
thermodynamics, and derive the condition for the second law to be satisfied. 

Next, in the following sections, we develop the observational investigations on dark energy 
and modified gravity. In Sec. X, we discuss the basic ideas and concepts of cosmography in 
order to compare concurring models with the Hubble series expansion coming from the scale 
factor. In particular, we derive cosmographic parameters without choosing any cosmological 
model a priori. 

In Sec. XI, we examine how it is possible to connect F(R) gravity by cosmography and it 
is possible to reproduce by it the most popular dark energy models as so called Chevallier- 
Polarski-Linder (CPL) 35|, |36J, which is a parameterization of the EoS for dark energy, or 
the ACDM. 

In Sec. XII, we show, as examples, how it is possible to constrain F(R) models theoreti- 
cally. However, the approach works for any dark energy or alternative gravity model. 

In Sec. XIII, we discuss constraints coming from observational data. It is clear that the 
quality and the richness of data play a fundamental role in this context. 

Finally, conclusions with the summary of this review are presented in Sec. XIV. 
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II. THE A COLD DARK MATTER (ACDM) MODEL 



The action of the ACDM model in general relativity is described as 

S = J d 4 x^^ (R - 2A) + J d A xC M (g»u, *m) , (2.1) 

where R is the Ricci scalar, g is the determinant of the metric tensor g^ v , A is the cosmo- 
logical constant, and Cm with a matter field \I/m is the matter Lagrangian. The EoS of the 
cosmological constant, which is the ratio of the pressure Pa to the energy density pa of the 
cosmological constant, is given by 

w A = ^ = -1 . (2.2) 
Pa 

We assume the 4-dimensional FLRW metric which describes the homogeneous and 
isotropic univese 



ds 2 = -dt 2 + a 2 (t) 



dr 



2 



+ rW 



(2-3) 



1 -Kr 2 

where a(t) is the scale factor, K is the cosmic curvature (K = +1,0, —1 denotes closed, flat, 
and open universe, respectively), and dVt 2 is the metric of 2-dimensional sphere with unit 
radius. The redshift z is defined as z = a^ja — 1 with ao = 1 being the current value of the 
scale factor. 

In the FLRW background (12. 3p . the Einstein equations are given by 

k 2 A K . . 

H 2 = -p M + --- 1 2.4 
S S a z 

H = -^-(pm + Pm) + 4> ( 2 - 5 ) 

2 GT 

where H = a/a is the Hubble parameter and the dot denotes the time derivative of d/dt. 
The fractional densities of dark energy, non-relativistic matter, radiation and the density 
parameter of the curvature are defined as 



O - PDE n _Pm n _Pr n — K 

Perk Perk Pcrit \ aH ) 



where p^ it = 3Hq/k 2 is the critical density with Hq being the current Hubble parameter. 
By combinning Eq. (12. 4ft with the quantities in (12. 6ft . we find 

^DE + fi m + fi r + fix = 1 , (2.7) 
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If there only exists cosmological constant A, i.e., />m = and Pm = 0, from Eq. (12. 5p we 
have H = H c = constant, so that de Sitter expansion can be realized. We also note that by 
comparing Eq. (12 .4p with Eq. (12. 2p . we obtain p\ = A/k 2 = —Pa. 

In addition, the scale factor is expressed as 

a = a c e Hct , (2.8) 

where a c (> 0) is a positive constant. 

In this section, we present the observational data of SNe la, BAO and CMB radiation, 
which supports the ACDM model (for the way of an analysis of observational data, see, 



e.g. 



3* 



A. Type la Supernovae (SNe la) 



With SNe la observations, we find the luminosity distance di as a function of the redshift 
z. We define the theoretical distance modulus as /ith(^i) = 5 log 10 Di(zi) + fi . Here, 
D L = H dL is the Hubble- free luminosity distance and di = \/L s / (An J 7 ), where L s and T 
are the absolute luminosity of a source and an observed flux, respectively, is the luminosity 
distance. Moreover, /i = 42.38 — 51og 10 /i, where h = i7 /100/[kmsec _1 Mpc" 1 ] 
express Dl as 



38 



. We 



D L (z)=H d L = (l + z)f K (y) 



(0) 
K 



y 



dz' 
EJz 7 ' 



(2.9) 
(2.10) 

(2.11) 



with 



E(z) 



H(z) 



(2.12) 



3(1+Wde W 



nP(l + zY, (2.13) 



' (i + zf + ni 0) (i + zf + exp 

where fir ' 1 = ^7 (1 + 0.2271A' r e fj) with f2 7 3 ' ) being the present fractional photon energy 



density and N e g = 3.04 the effective number of neutrino species ^|. In what follows, 
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the superscription "(0)" represents the values at the present time. Moreover, in deriving 
Eq. (I2.13P we have used the continuity equations 



p j + 3H(p j + P j ) = 0. 



(2.14) 



where j 
by 



"DE", "m" and "r", and P m = 0. Furthermore, f K (y) in Eq. fl2TTUD is described 



f K (y) 



siny for K 
y for K 

sinh y for K 



+ 1. 

0, 

-1. 



(2.15) 



By using Eqs. (12.91) . (12 . 1 3[) and (I2.15p . for the flat universe, we have Dl(z) = 
(1 + z) J* dz' /E(z'), where E{z) in Eq. (I2.13P with Vt^ = 0. From this relation, we find 
H(z) = {(d/dz) [di(z)/ (1 + z)]}~ . Accordingly, for z < 0(1), the cosmic expansion his- 
tory can be obtained through the measurement of the luminosity distance di{z). We expand 
Eq. (12.91) around z = as 



D L {z)=H d L = z+[l 



1 dE(z = 0) 

2 ~dz 



z 2 +0(z 3 



H-- ( 1 - 3w DE Q 



(o) 

DE 



n 



z 2 +0(z 3 ) 
(2.16) 



where in deriving the second equality we have used Eq. f !2. 13|) . It is clearly seen from 
Eq. (I2.16P that suppose there exists dark energy, > ari d "^de < 0, and hence the 
luminosity distance becomes large. We note that since the universe is very close to flat as 
—0.0179 < fig 1 < 0.0081 (95% confidence level (CL)) 3], even though the universe is open 
(tlx > 0), the change of the luminosity distance is small. 

By applying the Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Obser- 
vations data [4J, the latest distance measurements from the BAO in the distribution of 
galaxies, and the Hubble constant measurement, for a flat universe, the current value of 
a constant EoS for dark energy has been estimated as wde = —1-10 ± 0.14 (68% CL) in 
Ref. [4]. Moreover, as an example of a time-dependent EoS for dark energy, for a linear 



form wr, E (a) = w DE o + w DEa (1 - a) 



35 



36j, where u>deo an d w^Ea are the current value 



of wde and its derivative, respectively, by using the WMAP data, the BAO data and the 
Hubble constant measurement and the high-redshift SNe la data, u>deo and wn Ea have been 
analyzed as u>deo = —0.93 ± 0.13 and wr,Ea = — 0-411q>i (68% CL). This form is called as 
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the CPL model 35 



|. Consequently, for the flat universe, the various recent observational 
data are consistent with the cosmological constant, i.e., u>de = — 1- 

We also mention the way of analyzing the \ 2 °f the SNe la data, given by Xstss = 
[fJ'obs(zi) — ^th( z i)] 2 /of with /i obs being the observed distance modulus. In the following, 
the subscripts "th" and "obs" mean the theoretical and observational values, respectively. 



We expand Xsn as 

[39] 



Xsn = ^sn - 2/io-BsN + /^oCsn , (2-17) 
A SN = ^2 ^ obs ^) ~ ^ th ( z " = °)J ^ 18) 



-dsn = 2 ' (z.iyj 

i a * 

£ i ■ < 2 ' 2() ) 



CsN 

a, 



Since we do not know the absolute magnitude of SNe la, we should minimize x|n with respect 
to /io related to the absolute magnitude. We describe the minimum of Xsn with respect to 
as x|n = ^sn — -Bsn/Csn by using, for example, the Supernova Cosmology Project (SCP) 



Union2 compilation with 557 supernovae, whose redshift range is 0.015 < z < 1.4 401 ]. 



B. Baryon Acoustic Oscillations (BAO) 



Baryons couple to photons strongly until the decoupling era, and therefore we can detect 
the oscillation of sound waves in baryon perturbations. The BAO is a special pattern in the 
large-scale correlation function of Sloan Digital Sky Survey (SDSS) luminous red galaxies. 
Hence, we can use the BAO data to explore the features of dark energy. 

We measure the distance ratio d z = r s (zd)/Dv{z) through the observations of the BAO. 
The volume-averaged distance Dy(z) Q and the proper angular diameter distance Da(z) 
for the flat universe are defined by 

1/3 

'l + zYD\(z) 



D v {z) = 
D A (z) = 



Hz) 



dz' 



l + zj H{z')- 



(2.21) 
(2.22) 



13 



The comoving sound horizon r s (z) is expressed by 

;l/(l+z) 

r s (z) 



1 

7s 



-dz' 



^1+ (3ftf/4fif) /(1 + 2') 



(2.23) 



H(z' 



with £1^ = 2.2765 x 10~ 2 h~ 2 and = 2.469 x 10~ 5 h~ 2 being the current values of baryon 
and photon density parameters, respectively ^|. The fitting formula of the redshift Zd at 
the drag epoch 4JJ, when the sound horizon determines the location of the BAO because 



,(0) 



baryons are free from the Compton drag of photons, is represented as 



42| 



1291(0^/i 



(°) l,2\0.251 



. + 0.659(^ } ^ 2 ) 
h = 0.313(^/i 2 )- - 419 

h = 0.238 {n^h 2 ) - 223 



\0.828 



6l (0?V 



1 + 0.607 (ft^h 2 ) 



6 2 



0.674 



(2.24) 

(2.25) 
(2.26) 



For 



0.276 and h = 0.705, we have z d » 1021. 
By using the BAO data from the Two-Degree Field Galaxy Redshift Survey (2dFGRS) 
and the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) 41|, we measure the distance 
ratio d z at two redshifts z = 0.2 and z = 0.35 as rf°= i 
0.1097 ± 0.0036 with the inverse covariance matrix, defined by 

30124 -17227 
-17227 86977 



0.1905 ± 0.0061 and dfj 035 



M B l 



(2.27) 



~obs "\ 



We express the x 2 of the BAO data as xIao = (^Jbao ~ ^bao) (^BAo)y ( x j h BAo " 
with Xj iB Ao = (^0.2,^0.35)- 



C. Cosmic Microwave Background (CMB) radiation 

By using the CMB data, we can derive the distance to the decoupling epoch z*(p± 
1090) [2I, and hence we constrain the model describing the high- 2 epoch. Since the ex- 
pansion history of the universe from the decoupling era to the present time influences on the 
positions of acoustic peaks in the CMB anisotropics, those are shifted provided that there 
exists dark energy. 

We define the angle for the location of the CMB acoustic peaks as 

0a = , (2.28) 



d^ {^z^ 
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with c ffi = di/ (1 + z) being the comoving angular diameter distance. The acoustic scale 



Ia 
ft 



43 
45 



44 



representing with the CMB multipole corresponding to 9 a and the shift parameter 



46[ are defined as 



Ia(z*j 



7l 

T A 

3tt 
"I 



'1 + 



-xDa{z* 



nfh 2 



in 



R,(a) 



U(z*) 



\ n?h 2 
4^ " 



AQ^h 2 



n 



(0) 



(o) 

K 



Hq sinh 




cq I 



1 + R s (a*) 



Rs(aeq) 



-1 



ft 



1 + z*) D A {z* 



(2.29) 
(2.30) 
(2.31) 

(2.32) 



with z* being the redshift of the decoupling epoch 44J, given by 



1048 



1 + 0.00124 (fiJV 

0.238 



-0.738 



1 + 91 {^h 2 ) 



2\92 



0.0783 ( Qf ] h 2 



9i 



92 



0.560 



l + 39.5[Q^h 2 



0.763 ' 



h 2i.i {n { ° ] h 2 



1.81 ' 



(2.33) 
(2.34) 
(2.35) 



In Eq. (I2.30p . Pb and p 7 are the energy density of baryons and photons, respectively. Fur- 
thermore, a* and a eq are the scale factors at the decoupling epoch and the radiation-matter 
equality time. 

It is seen from Eq. (12.301) that since the change of the cosmic expansion history from the 
decoupling era to the present time influences on the CMB shift parameter, the multipole Ia 
is shifted. The relation between all peaks and troughs of the observed CMB anisotropics is 
represented by [47] 

l n = l A {n- <f> n ) , (2.36) 

where n denotes peak numbers and <p n is the shift of multipoles. For instance, for the first 
peak n = 1 and for the first trough n = 1.5. Moreover, the WMAP 5-year results js} present 
the limit on the shift parameter of ft = 1.710 ± 0.019 (68 % CL). In the flat universe, we 
have 1Z(z*) = y o£?i?o (1 + z*) Da(z*)- Hence, the smaller is, namely, the larger 
is, the smaller ft is. For an estimation executed by the WMAP 5-year data analysis Q of 
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Tl = 1.710, n^h 2 = 0.02265, n^h 2 = 0.1369 and tt ( ° ] h 2 = 2.469 x 10" 5 , with Eq. ( |2"30]) 
we obtain I a = 299. By plugging this value and 0i = 0.265 47J into Eq. (I2.36p . we find the 
first acoustic peak li = 220, which is consistent with the CMB anisotropies observation. 

The data from Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) obser- 
vations [l| on CMB can be used. The \ 2 °f the CMB data is expressed as Xcmb = 

(<CMB " *S&») (McU), (*&UB " *3W *,:,CMB S (/,(,,), K(z,), z.) ^ 

being the inverse covariance matrix. It follows from the WMAP7 data analysis [4j that 
l A (z*) = 302.09, K(z*) = 1.725, = 1091.3, and 

f 2.305 29.698 -1.333 ^ 

29.698 6825.27 -113.180 • (2.37) 
-1.333 -113.180 3.414 



■^CMB 



V 



/ 



As a result, we find the xLai of a11 the observational data xLai = Xsn + Xbao + Xcmb- lt 
is known that there exists a fitting procedure called the Markov chain Monte Carlo (MCMC) 



approach, e.g., CosmoMC 
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Finally, we mention an issue of the origin of the cosmological constant A, which is one of 
the most possible candidates of dark energy because its presence is consistent with a number 
of observations. In fact, however, the origin is not well understood yet 
energy density p v is given by 

d 3 k Vk 2 + m 2 



The vacuum 



Pv 



(27T) d 

4nk 2 dk k 
(2tt) 3 2 



2 

K_ 
IQir 2 



(2.38) 
(2.39) 



In deriving Eq. ( 12. 38ft . we have used the zero-point energy of a field with its mass m, 
momentum k and frequency u, E = to/ 2 = \/k 2 + m 2 /2, and k c (<^ m) is a cut-off scale. 
The vacuum energy density is p v ~ 10 74 GeV 4 , provided that k c = Mp\. On the other hand, 
the current observed value of the energy density of dark energy is pde — 10 -47 GeV 4 . Thus, 
the discrepancy between the theoretical estimation and observed value of the current energy 
density of dark energy is as large as 10 121 . This is one of the most difficult problems in 
terms of the cosmological constant. There also exists another problem why the "present" 
value of the energy density of vacuum in our universe is extremely small compared with its 
theoretical prediction. 
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III. COSMOLOGY OF DARK FLUID UNIVERSE 



In this section, we present a description of dark fluid universe [28|-|30|. We concentrate 
on the case in which there is only single fluid which corresponds to dark energy in general 
relativity 

A. Basic equations 

We represent the equation of state (EoS) of dark energy as 

w D e = — = -1 , (3.1) 

P P 

with 

P = -P- f( P ) , (3.2) 

where f(p) can be an arbitrary function of p. In what follows, for the time being, we will 
investigate the evolution of the universe at the dark energy dominated stage, so that we can 
regard p = pde and P = Pde- For simplicity, the "DE" subscription will be omitted as long 
as there is no need to mention. We also remark that f(p) characterizes a deviation from the 
ACDM model. 

We assume the flat FLRW metric ds 2 = —dt 2 + a 2 (t) Yli=i 23 (dx 1 ) 2 , which is equivalent 
to Eq. f)2.3p with K = 0. In this background, the Einstein equations are given by 

# 2 = yP, (3.3) 
H = -^(P + P). (3.4) 

Equations (13.3j) and (13.41) correspond to Eqs. f)4.2p and ( 14. 3 p with = p, P^ = P, p M = 0, 
and Pu = (i.e., the case in which the scalar field <fi is responsible for dark energy and 
there is no any other matter) shown in Sec. II A, respectively. Furthermore, the continuity 
equation of the fluid is given by 

p + 3H (p + P) = 0. (3.5) 
From Eq. (13.51) . we find that the scale factor is described as 

1 r dp \ 



a = a c exp ( - J j , (3.6) 
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where a c is a constant. In addition, by combining Eqs. (I3.3P and ( 13. 5 j) with Eq. ( 12. 3p . we 
obtain 

Eqs. (|33) and (S3]) lead to 

- = (P + 3P) = ^ (2p + 3f(p)) . (3.8) 
a o o 



B. EoS of dark energy in various cosmological models 

Provided that there exists only single fluid of dark energy, (i) for the ACDM model, wbe = 
— 1 (f(p) = 0), (ii) for a quintessence model, —1 < wde < —1/3 (—2/3 < f(p)/p < 0), (hi) 
for a phantom model, u>de < — 1 (f(p)/p > 0). If our universe lies beyond u>de = — 1 region, 
then its future can be really dark. In other words, in finite-time phantom/quintessence 
universe may enter a future singularity. 



1. Finite-time future singularities and energy conditions 



In the FLRW background (I2.3p . the effective EoS for the universe is given by 12] 



where p c s = 3H 2 /k 2 and P c s = — (2H + 3H 2 ^j / k 2 correspond to the total energy density 
and pressure of the universe, respectively. When the energy density of dark energy becomes 
perfectly dominant over that of matter, one obtains u>de ~ ^eff- Fo r H < (> 0), w c r > 
— 1 (< —1), representing the non-phantom [i.e., quintessence] (phantom) phase, whereas 
w c s = —1 for H = 0, corresponding to the cosmological constant. It is not clear from the 
very beginning how our universe evolves, and if it ends up in a singularity. This should 
always be checked. 

The finite-time future sing ularities can be classified into the following four types 
* Type I ("Big Rip' : 



26 



49| ) singularity: In the limit of t — > t B , all the scale factor, 
the effective energy density and pressure of the universe diverge as a — > oo, p cS — >• oo and 
|P c fj| — > oo. This also includes the case that p c g and P c r asymptotically approach finite 



values at t = t s . 
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* Type II ("sudden" |50|,[5l|) singularity: In the limit oft—> t B , only the effective pressure 
of the universe becomes infinity as a — > a s , p c s —> p s and |P e ff| — > oo. 

* Type III singularity: In the limit of t — > t s , the effective energy density as well as the 
pressure of the universe diverge as a — > a s , p c fr _ > oo and |P c g| — > oo. 

* Type IV: In the limit of t — > t s , all the scale factor, the effective energy density and 
pressure of the universe do not diverge as a — > a s , p c fr — > and |P e ff| — > 0. However, higher 
derivatives of H become infinity. This also includes the case that p e ff and/or |P e ff| become 



finite values at t = t„ 



52]. 



Here, t s , a s (^ 0) and p s are constants. In Ref. 53], the finite-time future singularities 



in F(R) gravity have first been observed. Furthermore, the finite-time 
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uture singularities 



55| . In particular, 



in various modified gravity theories have also been studied in Refs. 
it has recently been demonstrated that the finite-time future singularities can occur in the 
framework of non-local gravity 



in Ref. 



571 ] as well as in f(T) gravity in Ref. [58[. Also, 



various studies on the finite-time future singularities have recently been executed, e.g., in 
Ref. [59]. 

Moreover, there exist four energy conditions. 



(a) The null energy condition (NEC): 



p + P > 0. 



(3.10) 



(b) The dominant energy condition (DEC): 



p>0, p±P>0 



(3.11) 



(c) The strong energy condition (SEC): 



p + 3P>0, p + P>0 



(3.12) 



(d) The weak energy condition (WEC): 



p>0, p + P>0 



(3.13) 



If wde = P/p < — 1, the Type I (Big Rip) singularity appears within a finite time. In such 
a case, all energy conditions are violated. On the other hand, even when the strong energy 
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condition in Eq. A3. 12[) is satisfied, the Type II (sudden) singularity can appear [50|. We 
remark that if the Type I (Big Rip) singularity exists, all the four conditions are broken, 
whereas if there exist the Type II, III or IV singularity, only a part of four conditions is 
broken. Hence, it is expected that the origin of the finite-time future singularities is somehow 
related with the violation of all or a part of energy conditions. 
We mention that 'W singularity has been studied in Refs. 



60 



62[ and parallel-propagated 



(p.p.) curvature singularities 63| have earlier been investigated. For the 'W singularity, 
when t — > t s , a — > a s , p c s — > 0, |P e g| — > 0, and the EoS for the universe becomes infinity. 



2. Example of a fluid with behavior very similar to the ACDM model 



To begin with, we investigate a fluid which behaves very similar to the ACDM model 



f(p) = P q - 1 



(3.14) 



where q is a constant (\q\ <C 1). In this case, from Eq. (13.11) we see that |u>de — 1| «C and 
therefore this fluid model can be regarded as almost the ACDM model. For this fluid, by 
substituting Eq. (13.141) into Eqs. (I3.6P and (13. 7p . we obtain 



a ~ a c 



p-l 



1/(69) 



'-1 



(3.15) 
(3.16) 



where in deriving the approximate equalities we have used \q\ 1. 



3. Generalized Chaplygin gas (GCG) model 



Next, we examine the generalized Chaplygin gas (GCG) model proposed in Refs. 22L |23|. 
This model is a proposal to explain the origin of dark energy as well as dark matter through 
a single fluid. 



where A(> 0) is a positive constant and u is a constant. For u 



original Chaplygin gas model 



22| 



(3.17) 

1, Eq. fEUTD is for the 
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By combining Eq. (I3.17P and the continuity equation (13. 5 j) . we find 



P 



A + 



B 



l/(l+u) 



(3.18) 



a 3(l+u) 

where B is an integration constant. It follows from Eq. (13. 18ft that the asymptotic behaviors 
in the early universe a C 1 as well as the late universe o> 1 are described as 



p ~ Ba 3 , for a <C 1 , 
p ~ A 1/{1+u) , for a > 1 . 



(3.19) 
(3.20) 



In the early universe, the energy density evolves p oc a~ 3 and hence this fluid behaves as non- 
relativistic matter, i.e., dark matter, whereas in the late universe, p approaches a constant 
value of w 4. 1 /( 1 +") and therefore it can correspond to dark energy. As a result, the generalized 
Chaplygin gas model can account for the origin of both dark matter and dark energy. 

ose to the ACDM model 



It is known that the generalized Chaplygin gas model is very c 
by analyzing the gauge- invariant matter perturbations 5u in Refs. 6J, |65j, which satisfy [65 1 

2 



5 M + (2 + 3c" - 6w BE ) H5 



(1 - 3c s 2 - 3w 2 DE + 8w BE ) H* 



c s k 



5 M = 0, (3.21) 



where k is a comoving wavenumber and wbe — P/p with Eq. (13 . 1 7[) is the EoS of the 
generalized Chaplygin gas. In addition, c s is the sound speed, defined by 

dP 



dp 



-uw DE 



(3.22) 



When z 3> 1, both i«de and c s approach to zero. Thus, at the matter-dominated stage 
c s 1, and after it c s grows. If u > 0(< 0), c 2 > 0(< 0) because »de < 0. 

It follows from Eq. (I3.2ip that 5m grows due to the gravitational instability, provided 
that 



Kl<2 



3 f aH\ 



(3.23) 



On the other hand, if \c%\ > (3/2) (aH/k) 2 , the rapid growth or damped oscillation of 5m 
occurs, depending on the sign of c 2 . Around the present time when | Wr>F, I = 0(1) and thus 



| ix. j . The relation (I3.23P imposes the following constraint on u 65j: 



\u\ < 10" 



(3.24) 



Since the case of u = is equivalent to the ACDM model, the viable generalized Chaplygin 
gas model would be very close to the ACDM model. 
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4- Coupled dark energy with dark matter 



Related to the generalized Chaplygin gas model, in which dark energy and dark matter 
are interpolated each other, we explain a coupled dark energy with dark matter and its 



cosmological consequences (for recent discussions, see |66l. |67|). 



The equations corresponding to the conservation law are written as 

Pde + 3H (1 + wde) Pde = -Qpde , (3.25) 
p m + 3H (p m + P m ) = +Qp DE , (3.26) 

where u>de = -Pde/pde- Here, Q describes a coupling between dark energy and dark matter, 
which is assumed to be a constant. We note that the subscription "m" represents quantities 
of cold dark matter (CDM), and therefore P m = 0. It is important to emphasize that the 
continuity equation for the total energy density and pressure of dark energy and dark matter, 
which is the summation of Eqs. (13.251) and (13.261) . is satisfied. 
The solutions for Eqs. ( I3.25P and (13.261) are given by 



Pde = PBE(Q)a-^ 1+WDE) e-^ , (3.27) 
Pm = Qa- 3 J dtp DE[0) a- 3w ™(t)e- Qi , (3.28) 



where Pde(o) is an integration constant. 

In the flat FLRW background ( I2.3p . from Eqs. ( 13. 3 h and ( 13. 3ft the Einstein equations are 
expressed as 

3 

\ (2H + 3H 2 ) = P DE (3.30) 



2 H 2 = p DE + p mj (3.29) 



wbePbe = w DE p DEm a-^ 1+w ^e- Qt , (3.31) 



where in deriving the second equality in Eq. (I3.3ip we have used Eq. (I3.27p . We can obtain 
an exact de Sitter solution of Eq. (13 . 3 1 [) as 

a = a e m , (3.32) 

with 

H = — . (3.33) 

3(1 + u;de) 1 ; 
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Here, ao is a constant, and obeys 



3 

72 



Q 



w D EpDE(o)a 3.34 

3 (1 + wbe) 

If the universe is in the phantom phase, i.e., u>de < —1, H in Eq. ( 13. 33ft is positive. Hence, 
Eq. (I3.34p has a real solution, so that ao can be a real number. 

On the other hand, when only dark energy exists and therefore there is no direct coupling 
of dark energy with dark matter, i.e., p m = and Q = 0, from Eqs. (I3.25P and (I3.29P we 
find /?de + 3H (1 + wde) Pde = and (3/k 2 ) H 2 = p DE , respectively. These equations can 
have an expression of H as 

2/[3(l + ^ E )] 

U-t K ' 

which leads to a Big Rip singularity at t — t s . Thus, one cosmo logical consequence of a 

coupling between dark energy and dark matter is that a de Sitter solution could be realized, 

and not a Big Rip singularity. 

Another cosmological consequence of a coupling between dark energy and dark matter is 

to present a solution for the so-called coincidence problem, i.e., the reason why the current 

energy density of dark matter is almost the same order of that of dark energy. H in Eq. ( 13. 33ft 

may be taken as the present value of the Hubble parameter, which is given by H p = 2.1a x 



10" 42 GeV 



with h = 0.7 



In this case, from Eqs. ( I3.27P and (I3.32|) we see that the 



energy density of dark energy becomes constant 

Pde = Pde(o)«o 3(1+u,de) . (3.36) 
By substituting Eq. (13. 36ft into Eq. (I3.28j) . we acquire 

Pm = Pm(0)GT 3 - (1 + W DE ) PDE(0)«0 3(1+U,DE) > ( 3 - 37 ) 

where p m (o) is an integration constant. By plugging Eq. (I3.37P into Eq. (I3.29P and using 
Eq. (I3.34p . we find that p m (o) = 0. Hence, the energy density of dark matter is also constant 

as 

p m = - (1 + w DE ) p D E(o)ao 3(1+ ™ de) = - (1 + w DE ) pde , (3.38) 

where the second equality follows from Eq. (I3.36p . Thus, provided that the de Sitter solution 
in Eq. (I3.32p is an attractor one, by taking 

^ = _(i + WDE ) * (3.39) 
Pde >j 
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from which we have 

^de~-^, (3.40) 
the coincidence problem could be resolved. Of course, this is rather qualitative presentation. 

C. Phantom scenarios 

According to the present observations, there is the possibility that the EoS of dark energy 
wde would be less than —1. This is called the "phantom phase". It is known that if the 
phantom phase is described by a scalar field theory with the negative kinetic energy, which 
is the phantom model, the phantom field rolls up the potential due to the negative kinetic 
energy. For a potential unbounded from above, the energy density becomes infinity and 
eventually a Big Rip singularity appears. 



1. Phantom phase 

In order to illustrate the phantom phase, we present a model in which the universe evolves 
from the non-phantom (quintessence) phase (u»de > — 1) to the phantom phase (u>de < — 1), 



namely, crossing of the phantom divide line of wde = — 1 occurs 69]. 
The scale factor is expressed as 

a = «.(^)", (3.41) 

where a c is a constant, n(> 0) is a positive constant, t s is the time when a finite-time future 
singularity (a Big Rip singularity) appears, and we consider the period < t < t s . In this 
model, for t <C t s , a(t) behaves as a « t n and hence u>de = — 1 + 2/ (3n) > — 1, whereas for 
f ~*s, »de = -1 -2/ (3n) < -1. 

It follows from Eq. ( 13.41D that the Hubble parameter H is given by 

By combining Eq. (13.31) with Eq. f !3.42p . we find 
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H as well as p becomes minimum at t — t s /2 and those values are 



# min = — , (3.44) 



2 



48n 



(3.45) 



Moreover, p is written as 




--a/^, (3-46) 
3n 2 nt s V 3 

where we have removed t by using the relation between p and t in Eq. (I3.43p . Here, the 
plus (minus) sign denotes the expression of p for t > t s /2 (0 < t < t s /2), i.e., the phantom 
(non-phantom) phase. This implies that the energy density of dark energy increases in the 
phantom phase. Substituting Eqs. (13 .2p and (I3.46P into the continuity equation of dark 
energy (I3.5p . we obtain 



As a consequence, it is necessary for the EoS to be doubled valued, so that the transition 
from the non-phantom phase to the phantom one can occur. Furthermore, we see that 
/(Pmin) = 0. This means that at the phantom crossing point u>de = —1, both H and p 
have those minima, which can also be understood from the definition of wde in Eq. (13.11) . 
In addition, when a Big Rip singularity appears at t — t 8 , from Eq. (I3.47P we find that f(p) 
evolves as f(p) — 2p/ (3n), and by using this relation and Eqs. (13.21) as well as (13.11) we see 
that P = —p — 2pj (3n) and u>de = — 1 — 2/ (3n), which is a constant. On the other hand, 
in the opposite limit of t — > 0, we also have a constant EoS as u>de = — 1 + 2/ (3n). 

When the crossing of the phantom divide occurs, /(/?)= 0. In order to realize this 
situation, the integration part of J dp/f(p) on the right-hand side (r.h.s.) of Eq. (13. 6p 
should be finite. Thus, f(p) need to behave as 

f(p) ~ fc (p - PcY , < a < 1 , (3.48) 

where we have used the condition f(p c ) = 0. In general, it is necessary for f(p) to be 
multi- valued around p = p c because < s < 1. 
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2. Coupled phantom scenario 



Next, we explore another phantom scenario in which dark energy is coupled with dark 
matter. We suppose that dark energy can be regarded as a fluid satisfying Eq. (I3.25p . Since 
the coupling Q can be described by a function of a, H, p m , p m , p DE and pde, we consider 
the case in which the ratio of p m to pde defined as r = p m /pDE is a constant, that is, Q is 
represented by scaling solutions. In this case, p m does not close to zero asymptotically. By 
using Eqs. (I3.25|) and (I3.26p . we obtain 



Q [ — + — ) - 3H (w m - mj D e) 
Pm Pde , 



(3.49) 



where w m = P m / p m . By taking f = in Eq. ( 13.490 . we acquire Q in the presence of scaling 
solutions as 

Q = 3H (w m - de) PDEPm • (3.50) 
Pde + Pm 

As a model which can be treated analytically, we examine Q given by 

Q = 5H\ (3.51) 

with 5 being a constant. Here, cold dark matter is regarded as a dust, and hence P m = 0. 
Moreover, we provide that /(pde) = Pde in Eq. (13. 2p . i.e., 

Pde = -2p DE • (3.52) 
By combining Eq. (I3.29P with Eqs. (I3.25P and (13. 26 p . a solution is given by 

H=l(\ + T^], (3-53) 



3 \t t s -t, 
4 fl 1 \ 1 , 

^ E= 3^G+v^)i^i' (3,55) 

with 

*- s w- (3 ' 56) 

We note that the same solution can be derived if Q is described as 

9Hp DE p m 

2 (PDE + Pm) 
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The expression of H in Eq. (I3.53P is the same as that in Eq. C I3.42[) with n = 2/3, and 
thus a Big Rip singularity appears at t — t s . Incidentally, in this case t s is not determined. 
We also mention that Eqs. (I3.54p and (I3.55P lead to r = (t s — t) /t, which implies that this 
ratio is not a constant but changes dynamically and eventually becomes zero, and hence the 
solutions in Eqs. fl3.53p - fl3.55p do not correspond to scaling solutions. This is reasonable 
because around a Big Rip singularity, the energy density of dark energy becomes dominant 
over that of dark matter completely. 



3. Little Rip scenario 



We study a Little Rip scenario |70l477l| . which corresponds to a mild phantom scenario. 
The Little Rip scenario has been proposed to avoid the finite-time future singularities, in 
particular, a Big Rip singularity within fluid dark energy. In this scenario, the energy density 
of dark energy increases in time with w-de being less than —1 and then wbe asymptotically 
approaches wde = —1- However, its evolution eventually leads to a dissolution of bound 
structures at some time in the future. This process is called the "Little Rip". 

A sufficient condition in order to avoid a Big Rip singularity is that a(t) should be a 
nonsingular function for all t. We suppose a(t) is expressed as 

a(t) = e /w , (3.58) 

where fit) is a nonsingular function. It follows from Eq. (13. 3p that p = 3 (a/a) 2 = 3f 2 . 
The condition for p to be an increasing function of a is that dp/ da = (6/a) // > 0. This 
condition can be met provided 

/>0. (3.59) 

Thus, in all Little Rip scenarios, a(t) is described by Eq. (13.31) with / satisfying Eq. (13.591) . 

We derive the expression for p as an increasing function of the scale factor a and examine 
the upper and lower bounds on the growth rate of p(a) which can be used to judge whether 



a Big Rig singularity appears. By defining N = In a, Eq. (13.31) is rewritten to 70J 




dN . (3.60) 



p(N) 

The condition for the appearance of a finite-time future (Big Rig) singularity to be avoided 
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is that it takes a Big Rig singularity infinite time to appear. This means that 



,00 x 
JNc \/~p{ 



dN -> 00 



(3.61) 



where N c = lna c . For the case P = —p — Ap 1 ! 2 with A being a constant, we have 

2 



P_ 

Pi 



3A 1 a 
m — 



(3.62) 



2^pl \a c 

where p\ is a constant and p — p 1 and a = a c at a fixed time t\. In this case, A > 
is required, so that there can be exist the phantom phase as wbe < — 1- Moreover, the 
expression of p as a function of t is given by 



£_ _ e V3A(t-ti) 

Pi 

Furthermore, by using wde = — 1 — Ap -1 / 2 and Eq. (I3.62p we acquire 

-1 -1 



(3.63) 



wbe 



2 V a J A 



(3.64) 



which can also be derived from the relation (a/p) (dp/da) = — 3 (1 + wde)- In addition, by 
eliminating p from Eqs. f)3.62p and (I3.63P we find 

2v/pT 



a 

— = exp 



3A 



(3.65) 



In this expression, by comparing Eq. (I3.65P with Eq. (I3.58p . we obtain 



2 v / 3A/2(t-ti)-l 



+ In a c 



(3.66) 



In this case, / = (y / piv4/2) e( v ^ j4//2 )^ - * 1 ^ > because of A > 0, and hence the condition 
in Eq. ( I3.59P is satisfied. We mention that the single model parameter A satisfies the 
observational bounds estimated by using the Supernova Cosmology Project 40(. The best 
fit value is given by A = 3.46 x 10~ 3 Gyr" 1 and the range for 95 % Confidence Level fit is 
-2.74 x 1(T 3 Gyr^ 1 < A < 9.67 x KT 3 Gyr^ 1 |70|. 

From Eqs. f)3.64p and f)3.65p . we see that when t ~ ti, a ~ a c and w DE < —1, i.e., the 
universe is in the phantom phase. As the universe evolves, a increases, and when t <C t±, a 
becomes very large and thus wbe asymptotically becomes close to —1. However, it takes a 
as well as p infinite time to diverge due to Eq. f)3.6ip . a Big Rip singularity cannot appear 
at a finite time in the future. This means that a Big Rip singularity can be avoided. 
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As further recent observations on Little Rip cosmology, in Ref. 7l| it has been demon- 
strated that a Little Rip scenario can be realized by viscous fluid. On the other hand, in 
Ref. 72| a new interpretation of a Little Rip scenario by means of an inertial force has been 



presented. It has been shown that a coupling of dark energy with dark matter can eliminate 
a little rip singularity and an asymptotic de Sitter space-time can appear. Moreover, a scalar 
field theory with realizing a little rip scenario has been reconstructed. 

We investigate the inertial force -Fi ncrt on a particle with mass m in the context of the 
expanding universe. When the distance between two points is /, by using the scale factor a, 
the relative acceleration between those is described as la/ a. If a particle with mass m exists 
at each of the points, an inertial force on the other mass would be measured by an observer 
at one of the masses. Thus, F inert is expressed as [to|, \l\ 



inert 



ml- =ml[H + H' 

a 



(3.67) 



where in deriving the second equality in Eq. ( 13. 67ft we have used Eqs. ( 13. 3ft and ( 13. 4ft . If 
-Finert > and -F incrt > F h , where F h is a constant force and bounds the two particles, the 
two particles becomes free, i.e., unbound, so that the bound structure can be dissociated. 
If the Hubble parameter is expressed as 



H 



(3.68) 



(t s - ty 

where h s (> 0) is a positive constant and q(> 1) is larger than or equal to unity. In this 
case, a Big Rip singularity appears in the limit of t — > t s . By substituting Eq. ( I3.68P into 
Eq. (j3TBTj) . we obtain [58] 

q h s 



mlh* 



oo , when t — > t s . 



(3.69) 



i(t s -ty +1 (t s -tr 2 r 

Since in the limit of t — > t s , H and H diverge, -Fi nert becomes infinity. The important point 
is that for a Big Rip singularity, -Fi ncrt diverges in the "finite" future time. 

On the other hand, a representation of the Hubble parameter realizing Little Rip scenario 
is described by [72] 

H = H LR exp(£t) . (3.70) 
Here, Ht,^(> 0) and £(> 0) are positive constants. By plugging Eq. ( I3.70p into Eq. (I3.67p . 



we find 



finert = mlH LR [£ + H LR exp (£t)] exp (£t) — > oo , when t ->■ oo 



(3.71) 
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In the limit of t — » oo, the divergence of H and H leads to the consequence of -Fi nert — > oo. 
Thus, a "Rip" phenomenon happens in both cases of a Big Rip singularity and Little Rip 
scenario. It is remarkable to note that for Little Rip cosmology, it takes Fi ner t the "infinite" 
future time to become infinity. This feature is different from that in case of a Big Rip 
singularity, which is the finite time future singularity. 

In Little Rip cosmology, as a demonstration we estimate the value of F^ PT t at the present 



time t n « Hq \ where H n = 2.1h x 10" 42 GeV 
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with h = 0.7 



681 ] is the present 



value of the Hubble parameter. The bound force F^ s for the Earth-Sun system is given 



by F b ES = GM e M /r 2 _ 



4.37 x 10 16 GeV 2 , where M ffi = 3.357 x 10 51 GeV 



Earth mass and M c . 



© 



1.116 x 10 57 GeV 



7.5812 x 10 26 GeV _1 
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is the mass of Sun, and r, 



©-0 



is the 
1AU = 



is the Astronomical unit corresponding to the distance between 



Earth and Sun. We take £ = H , m = M e and I 



-©■ 



In this case, from Eq. (I3.7ip 



we acquire F incrt = 2.545 x 10 78 e# 2 [(H LR /H ) + e (H LR /H ) 2 ] , where e = 2.71828. In 
order for the current value of F incrt to be larger than or equal to F^ s , H LR should be 
H LR > 4.82 x 10~ 30 GeV. 



4- Pseudo-Rip cosmology 

In addition, as an intermediate cosmology between the ACDM model, namely, the cosmo- 
logical constant, and the Little Rip scenario, more recently the Pseudo-Rip model has been 
proposed in Ref. 78j. In this case, in the limit of t — > oo the Hubble parameter tends to a 
constant asymptotically. In other words, the Pseudo-Rip cosmology is a phantom scenario 
and has a feature of asymptotically de Sitter universe. 



A description realizing Pseudo-Rip cosmology is given by 58 



H(t) = # PR tanh (J-j . (3.72) 

Here, Hp R (> 0) is a positive constant. It follows from Eq. (13. 72ft that we find 

a = a PR cosh T^-J , (3.73) 

with cxpr(> 0) being a positive constant. In the limit of t — > oo, H(t) — > Hp R < oo and 
H(t) increases monotonically in time. Hence, the universe evolves in the phantom phase 
and eventually goes to de Sitter space-time asymptotically. By combining Eqs. (13.721) and 
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(13. ip . we obtain 

3t -MpRsmh (t/to) 

This implies that since H(t) = H PR / [t cosh 2 (t/t )] > 0, the universe always evolves within 
the phantom phase as wde < — 1, and eventually, when t — > oo, w^e — > — 1, similarly to 
that in the Little Rip scenario. By comparing the observational present value of u>de as 
w DE = -1.10 ± 0.14 (68% CL) |4| with w DE in Eq. (|3~74l at t = t ~ # _1 > we obtain the 
constraint # PR > 2.96 x 10~ 42 GeV. 

Furthermore, by using Eq. (13.671) and (13.721) . we acquire js8 | 



inert 



mlH 



PR 



1 



iocosh (t/to] 



+ ifp R tanh z ( ^ 
In 



1 inert ,00 ^ {J<J 



when t 



with 



inert ,00 



mlHl K 



00 , 
(3.75) 

(3.76) 



This means that in Pseudo-Rip cosmology, -Fmert becomes finite asymptotically in the limit 
of t — > 00. This originates from the fact that when t — >■ oo, H — > i^pR and H — > 0. To 
realize the Pseudo-Rip scenario, the relation F^. t ^ > F^ s should be satisfied, so that the 
ES system could be disintegrated much before the universe asymptotically goes to de Sitter 
space-time. For m = M e and / = r e _ , we find H PR > ^/GM Q /r^_ Q = 1.31 x 10~ 31 GeV. 



This constraint is much stronger than that obtained from the present value of wde as i^pR > 
2.96 x 10~ 42 GeV shown above. In Appendix A, it is examined how strong the inertial force 
can constrain the EoS of dark energy. 



D. Finite-time future singularities 

In the ACDM model, f(p) = in Eq. (13.11) and hence wde = — 1- In a quintessence 
model, the type II, III and IV singularities can occur. On the other hand, in a phantom 
model, the singularities of Type I and type II can appear. We explicitly demonstrate that 
the EoS in Eq. (13.11) can lead to all the four types of the finite-time future singularities. 
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1. Type I and III singularities 



We examine the case in which f(p) is given by 

f(p) = Ap a , (3.77) 



where A and a are constants. By using Eq. (13. 6p . we obtain 

,1-a 



a = a c exp 



P 



3(1 -a) A 



(3.78) 



For a > 1, when p — > oo, a approaches a finite value, whereas, for a < 1, when p — > oo, if 
A > (A < 0), a diverges (vanishes). We note that in this sub-subsection (i.e., Sec. Ill D 
1), we concentrate on the limit of p — > oo 1 . We also mention the case in the opposite limit 
of p -> in Sec. Ill E 2. 

Moreover, from Eq. (13 .7p we have 

111 P r 1 / 

t = t s + , for a = - . (3.80) 

For a > 1/2, p can diverge within the finite future or past t = t s . Meanwhile, for a < 1/2, 
it takes infinite time for p to diverge and therefore p — > oo in the infinite future or past. 

It follows from Eq. ( \3.2\\ with Eq. (13. 77ft that P = —p — Ap a . When p diverges, P also 
becomes infinity. By substituting Eq. (I3.77P into Eq. (13. ip . we find 

w DE = - = -l- Ap^ 1 . (3.81) 
P 

For a > 1, when p — > oo, if A > (A < 0), u>de — ► +oo (u>de — > — oo). On the other hand, 
for a < 1, when p — )• oo, if A > (A < 0), wde — > — 1 + (wde — > — 1 — 0), i.e., wde 
approaches —1. 

The above considerations are summarized as follows. 

(a) For a > 1, the Type III singularity can exist. If A > (A < 0), wde — > +oo 
(w DE -)• -oo). 



1 It has also been examined in Ref. that for a < 0, when p — > 0, there can appear the Type II singularity. 
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(b) For 1/2 < a < 1, if A > 0, there can exist the Type I (Big Rip) singularity. If A < 0, 
p — > oo, a — > 0. Such a singularity in the past (future) may be called a Big Bang (Big 
Crunch) singularity When p -> oo, if A > (A < 0), u> DE ->• — 1 + (u> DE -> -1 -0). 

(c) For < a < 1/2, there does not exist any "finite-time" future singularity 



f(p) = - — :vr> ( 3 - 82 ) 



Type II singularity 

Next, we investigate the following form of f(p): 

C 

(pc - py 

where C, p c and 7(> 0) are constants. We note that for p c = 0, (— 1)~ 7 C = A, 7 = —a, 
the expression in Eq. (I3.82p is equivalent to that in Eq. (13.77]) . We concentrate on the 
case in which p < p c . Since P = —p — C / (p c — p) 7 , when p — > p c , P diverges. Hence, 
R = 2k 2 (p — 3P) also becomes infinite. From Eq. (13. 6p . we obtain 

(Pc-P) 7+1 ' 



a = a c exp 



(3.83) 



3C (7 + I) 

Thus, when p — > p c , a is finite as a — > a c and hence a is also finite because H = a/a oc ^fp 
given by Eq. (13.31) . whereas a diverges. Furthermore, by using Eq. (13. 7p . we find 

t ~ t ^ Pc - P ^' + (3 84) 

S kCv^(7 + 1)' ^ 



where t s is an integration constant. Therefore, t — > t s when p — > p c . By combining Eq. ( 13.8^ 
into Eq. (13. ip . we have 



wbe = -1 - • (3.85) 



C_ 

P (Pc - py 

When p — >• p c , wve — > —00 (wde — > 00) for C > (C < 0). As a consequence, there can 
exist the Type II (sudden) singularity if f(p) is given by Eq. (13. 82ft . It is interesting to 
remark that for C < 0, when p becomes around p c , the strong energy condition in Eq. (I3.12p 
is satisfied. This means that the Type II (sudden) singularity can appear in the quintessence 
era. 

Next, we demonstrate a quintessence model as well as a phantom one in which the Type 
II singularity occurs Q|. We consider the case that f(p) is given by 

/(P)=i-^T-, (3-86) 

1 - P/Ps 
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with ((> 0) being a positive constant. In the limit of t — > t s , p approaches p s (> 0). It 
follows from Eq. (13.11) that if f(p) > with p > 0, which is realized for p < Ps with p 
being the present energy density of the universe, %e < — 1 ; an d hence the universe is in the 
phantom phase and p increases until the pressure of the universe diverges. While, for p > 0, 
if —2/3 < f(p)/p < 0, — 1 < u>de < — 1/3 and therefore the universe is in the non-phantom 
(quintessence) phase. For p > Ps, since f(p) < 0, this case can correspond to a quintessence 
model. The energy density p becomes small in time and finally a Big crunch happens at 
p = p s . We estimate how long it takes the Type II singularity to appear from the present 
time to- The remaining time is described by 

_2_ r s ° dx ' 

V3 4 c 2 



to 

x 



1^ 

2 s 



(3.87) 
(3.88) 



By using Eq. (13.61) . we acquire 

P_ _ 

Ps 
E 



1± +6£m(l + *) 

Po 

Ps 

e 



(3.89) 
(3.90) 

(3.91) 



Here, in Eq. (13.891) the "+ (— )" sign denotes the case that 9 > 1 (6 < 1) describing a 
quintessence (phantom) model. Moreover, the EoS for dark energy at the present time is 
expressed as 

WDE(O) = -1 



9(1-6)- 

In addition, from Eqs. (12. 9p . the luminosity distance di is described by 



d L = — (1 + z) 
Ho 



dz' 
Wz 7 



e{z) = dn^ii + zf + n^z 



(0) 



1± V( 1 -©) +6^111 (l + z) 



(3.92) 

(3.93) 

(3.94) 
(3.95) 



Furthermore, the remaining time from the present time to the appearance of the finite-time 
future singularity is given by 



ts ~ t(] 



Hn 



du' 



^ (1 + u>) JO® (1 + zf + flW e (fiO 



)(0) 



(3.96) 
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to to u 



where u = ao/a — 1. This variable changes from u = at t 
exp [— (1 — O) 2 / (6£)] — 1 at the time when P diverges, i.e., in the limit of t — > t s . We 
note that u is equivalent to the redshift z, and hence for u = z, g(u') in Eq. (I3.96P be- 
comes equal to g(z) in Eq. (I3.95p . As numerical estimations, for the phantom phase with 
wbe(o) — —1-10 and = 0.95, t s — t = l.lGyr, whereas for the quintessence phase with 
wde(o) = —0.98 and 6 = 1.05, t s — t = 5.79Gyr and td—t = 5.78Gyr, where is the time 
when the cosmic deceleration starts, i.e., a becomes zero from its positive value 
Hq 1 = 13.6Gyr has been used. 



79|]. Here, 



The deceleration parameter q^ec and the jerk parameter j are defined by 

1 d 2 a 1 [ld{H 2 ] 



80) 



J = 



aH 2 dt 2 
1 d 3 a 1 
aH 3 dt 3 ~ 2H 2 



1 

H 2 



2 dN 
d 2 (H 2 ) 



d{H 2 



2H' 



(3.97) 
(3.98) 



dN 2 dN 

where iV = — In (1 + z) is the number of e-folds and iV = at the present time t = t . The 
values of gdec and j at the present time t = to are written as 



Qdccit — to) = qd, 



ec(0) 

9 

-- (wbe 



9£ 



26(1-9) 

i) ngi 4 



n {0) + ■ 

DE > 
' 2^DE 



-fi (0) - 1 
, DE 1 



j(t = t ) 



Jo 
9 



9^ 



26(1 - 6) 



(w DE + 1) 



1 - 



1 - 
6 



1-6) 
(w DE + 1) 



i4 DE 



n m + i 

"de ^ L 



(3.99) 
(3.100) 
(3.101) 

(3.102) 



1 - 6 

If ifDE(o) ~ —1) except the limit of 6 — > 1, jo is not so different from its value for the 
ACDM model, where j is unity. Moreover, in the range —1.05 < u>de(o) < —0.95, we find 
4tc(o) - 0.16 < g dcc(0 ) < q d ^ c{0) + 0.16 Q, where g^ (0) = -1 is the value of g dec(0 ) for the 
ACDM model. Thus, it can be considered that the model in Eq. (I3.86P describing both the 
quintessence and phantom phases, where the Big Crunch and the Type II singularity even- 
tually happens, respectively, fits the latest supernova data from the Supernova Cosmology 
project 40|] well. Similar results have been obtained also in Refs. 8ll. I82I]. 

We remark that for the case describing the quintessence phase, the dissolution of the 
bound structure before the appearance of the Big Crunch cannot be realized. The inertial 
force is given by 

ml f 1\ 

-Fincrt = — ( WBE + 3 j P ' (3.103) 
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For the quintessence phase, the maximum value of F inert in Eq. (I3.103P is F^^. = mlp/3. 
The inertial force becomes small in time because the energy density of the universe in the 
quintessence phase decreases. 



3. Type IV singularity 



We further explore the case in which 



f(p) 



ABp 



a+/3 



(3.104) 

Ap a + BpP' v 1 

where A, B, a and /3 are constants. For a > f3, we see that f(p) — > Ap a when p — )■ 0, while 
f(p) — > Bp 13 when p — >■ oo. 

For a j£ 1 and /3 7^ 1, Eq. yields 



a = a c exp 



-a+l 



-/3+1 



[a -I) A (/3-l)5 



In addition, for a = 2/3 — 1, by using Eq. (13. 7p . we acquire 



r = -V3nA(t - t s 



- (4/3-3) /2 



2A 



-(2/3-l)/2 



(3.105) 



(3.106) 



4/3 - 3" (2/3-1)5' 
where t s is an integration constant. Equation (I3.106P is available for (3 ^ 1, /3 ^ 3/4, and 
/3 ^ 1/2. Furthermore, from Eq. (13. ip with Eq. f)3.104p . we have 

ABp^P- 1 



U>DE = -1 ~ 



Ap a 



For < /3 < 1/2, when p 0, 



-p - Bp? . 



(3.107) 



(3.108) 



Equation (I3.106P implies that when p — > 0, t — > t s . Hence, from Eq. (I3.109P we see that in 
the limit t -> t s , p and P ->■ 0. For a = 2/3 - 1, from Eq. f!3.104p we find 



hil- 

a r 



3(/3-l) 



P 



1-0 



2A B 



3(/3- 1)5 



(2/3 -1)5 
2A 



T 1 



with 



q=l 



1 



(3.109) 



(3.110) 



2/3 - 1 

Since q > 2, a — > a c in the limit t — > i s . Moreover, 5 and H are also finite. However, 
unless q is an integer, higher derivatives of H, d n H/dt n where n > — 1/ (2/3 — 1), diverges. 
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As a result, there can exist the Type IV singularity. In this case, i.e., for a = 2(3 — 1 with 
< (3 < 1/2, when p — > 0, from Eq. (13. 107[) we see that u>de — > oo (wde — > — oo) for B < 
(B > 0). 



It has been examined in Ref. 
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that if f(p) is expressed as Eq. (13.1041) . there can also 
exist the Type I, II and III singularities. For 3/4 < /3 < 1 with A > 0, the Type I singularity 
can appear. For A/B < with any value of (3 or for (3 < irrespective the sign of A/B, the 
Type II singularity can exist. For /? > 1, the Type III singularity can occur. 



E. Asymptotically de Sitter phantom universe 

We study an example of a fluid realizing asymptotically de Sitter phantom universe. 
We present an important model constructed in Ref. [l^] in which the observational data 
consistent with the ACDM model are satisfied, but it develops the dissolution of the bound 
structure. 

For convenience, we introduce a new variable x = yfp. By using x, Eqs. (13. 6 j) and (13. 7ft 
are rewritten to 

(2 r x'dx'\ 

a=a ' eKp \3Lm)' (3 - m) 

' ^ (3.112) 



V3kJ Xc f(x' 

where x c is the value of x at t = 0. If at x = Xf < oo, the integration in Eq. (13. 1 12j) becomes 
infinity. Hence, it takes the energy density infinite time to arrive at pf = x\. In other words, 
the cosmic expansion becomes exponential behavior asymptotically. In addition, the energy 
density approaches a constant value, i.e., the cosmological constant, whereas the EoS w for 
the universe always evolves within the phantom phase w < — 1. We examine a model in 
which f(x) is given by 

X \ 



f(x) = (3^ 



Xf/ 



(3.113) 



We suppose that the energy density of dark energy varies from zero to pf = x\. It follows 
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from Eq. (13.1111) that p is expressed as a function of the redshift z 

p{z) = * [1 - (1 + zf (1 - A)] 4/3 , 



7 = T P f 



3/3 ^-3/4 

1 

3/4 



poX 

pf) 



(3.114) 
(3.115) 

(3.116) 



where po is the energy density of dark energy at the present time to- Furtheremore, the 
current EoS for dark energy is given by 

27 1 - A 



^DE(O) 



(3.117) 



3 A 

By using a parmeter v varing from zero at t — to to unity in the limit of t — > 00, we can 
obtain the parametric description of the scale factor as 



a(v) 



t(v) =t + 



(1 _ 03/2)2/(37) 



In (v 2 



2^3 

7X f 



7T 



1 2V+ 1 

—= arctan = = 

\/3 V3 6\/3 



In (1 - v) 



(3.118) 
(3.119) 



6 V3 V3 6V3 3 

where in the limit of v — > 1 the last term in Eq. (I3.119P becomes dominant. With these 
parametric descriptions, the scale factor can accurately be expressed as 

(3.120) 

The energy density of non-relativistic matter and baryons behaves as p m = p m (o) (1 + - 2 ) 3 ? 
where p m (o) is the current value of the energy density of non-relativistic matter and baryons. 
Thus, it follows from Eq. (12.91) that the luminosity distance is written as 



/ 2 \ 2/(37) 


Xf (t - 1 ) 


UJ exp 


[ Vs \ 



d L = —(l + z) 
-"0 



h(z) = A" 4 / 3 [1 - (1 



^L 0) (1 + z? + ^lh(z') 
zf (1 - A)] 4/3 , 



-1/2 



dz' 



(3.121) 
(3.122) 



where we have neglected the contribution from radiation because it is much smaller than 
those from non-relativistic matter and baryons and dark energy. In case of the ACDM model, 



= Q\ and h(z) = 1 in Eq. (13. 121ft . where = Pa/ 'Pcrit- Provided that 7 ~ 0, where /3 
is very small, and A w 0, we find h(z) ~ 1. Thus, in this case it is impossible to distinguish 
this model with the ACDM model. From Eq. (13.1171) . we also have u>de(o) ~ — 1. On the 
other hand, if we take A = 0.5 and 7 = 0.075, u>de(o) = —1-05. According to the numerical 



(o) 
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the difference between the distance modulus /x in this model and that 



(ACDM) 



IS 



analysis in Ref, 

in the ACDM model is estimated as 5p = 5 log \di/ 'd^ CDM ^ < 0.016, where d 
the luminosity distance in the ACDM model, for f2m' ) = 0.28 and fij^ = Qa = 0.72. Since 
errors of the modulus of SNe are ~ 0.15, which are larger than the above estimation of S, 
this model can fit the observational data as well as the ACDM model. 
In this model, the Friedmann eqution (13. 3 j) is given by 



k 2 k 2 r 

H 2 = —(p + p m ) = — \p{ [l 



^(l-A)] 4/3 + Pm(0) e-^ 



3 vr rill/ 3 

By substituting Eq. (13.1231) into Eqs. (I3.97P and (I3.98p . we acquire 



<7dec(0) 

Jo 



2 r (A -1 - 1) n A + ^L 0) 



9 



7 2 (A- 1 -^ 



;1 



2 7 (A- 1 -!) 



(3.123) 

(3.124) 
(3.125) 



For A = 0.5, 7 = 0.075, VL { ^ = 0.28 and fi££ = fi A = 0.72, we obtain g dcc(0) = -0.54 and 
jo = 1-11. In the ACDM model, where A = 0, we have q'dec(o) — —0.58 and jo = 1- As a 
consequence, the values of gdec(o) an d jo in this model are near to those in the ACDM model. 

In addition, we explore whether the dissolution of the bound structure can occur. By 
using Eqs. (I3.3p . H3.4[) and (I3.67P F incrt is also represented as 



)(0) 



inert 



-m/— (pde(oJ + 3P D E(a)) = ml— [ 2p DE (a) + - 



da 



The dimensionless discription for F inert is given by 

1 



inert 



2pde(o) H : a 



Pde(o) V ^ a 
where Pde(o) is the current energy density of dark energy. 

The parametric discription of the dimensionless inertial force is given by 



(3.126) 



(3.127) 



F incrt (n) = 2A- 4 / 3 [1 - (1 + uf (1 - A)] 1/3 



and t — to is written by 
t-t 



27 
3 



1 



u 



A) 



(3.128) 



H 



du' 



u 



(1 + u'f + n^h(u' 



.A3 



(0) 



(3.129) 



Here, we have introduced the new variable u = ao/a — 1, which varies from at the present 
time t = t to —1 in the limit of t — > 00. It follows from Eq. (13. 128|) that in the limit of 
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t — > oo, we find F h 



u 



— y 2 A 4//3 . For the Earth-Sun system, the necessary value of F{ 



inert 



74j , which leads to 



for the disintegration of the bound structure is estimated as -Fi nert > 10 23 
the condition for the dissolution A < A min = 10~ 17 . Moreover, from Eq. ( 13.117)) we see that 
for A = A min = 10~ 17 , 7 can change from (in which u>de(o) — — 1) to 3.6 x 10 -18 (where 
w de(o) — —1-24, which is the observational lowest constraint [l]). As a consequence, if 
A < A m i n = 10 -17 and < 7 < 10 -18 , this model can be compatible with the observational 
data of SNe. Accordingly, in this model the dissolution of the bound structure can be 
realized, although it satisfies the observational data which is consistent with the ACDM 
model. 



F. Inhomogeneous (imperfect) dark fluid universe 



In this subsection, we explain inhomogeneous (imperfect) dark fluid universe by following 



Refs. 
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83| . We investigate so-called inhomogeneous EoS of dark energy, which the pressure 



has the dependence not only on the energy density but also on the Hubble parameter H. 



This idea comes from, e.g., a time dependent bulk viscosity in the ideal fluid 
modification of gravity. For a recent study of imperfect fluids, see 
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1. Inhomogeneous EoS 



An inhomogeneous expression of the pressure is described by 

P = -p + f(p)+G(H), 



(3.130) 



where G(H) is a function of H. We note that generally speaking, G is a function of the 
Hubble parameter H, its derivatives and the scale factor a. However, for simplicity we 
consider mainly the case that G depends on only H. By substituting Eq. (13.130)) into 
Eq. ( 13.51) . the continuity equation of the inhomogeneous fluid is represented by 



= p + 3H(f(p) + G(H)) . 



(3.131) 



By plugging H = k^J p/3, which follows from Eq. (1 3 . 3 j) for the expanding universe (H > 
0), into Eq. ( 13 .131)) . we obtain 



(3.132) 
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On the other hand, by combining Eq. (13 . 1 30[) with Eq. (13. 4p . we have a representation of 
G(H) in terms of / as 

G{H) = -/( — — ) + —H . (3.133) 

K K z 

2. Influences on the structure of the finite-time future singularities 

First, as a concrete example, we examine the case in which the relation between P and 
p in Eq. ( 13. 2 p is changed as 

P = w h p + w 1 H 2 , (3.134) 
w h + ^)p, (3.135) 

where is the homogeneous EoS, W\ is a constant, and in deriving Eq. ( I3.135P we have 
used Eq. (13.31) . Thus, for the inhomogeneous fluid, the inhomogeneous EoS is shifted 
from the homogeneous EoS u>h(= w) as 

K 2 Wi 

w — y w ih = w h H — . (3.136) 

From Eq. ( I3.136P , we see that if Wih > — 1 , there does not appear a Big Rip singularity, even 
for Wh < —1. In other words, provided that if W\ —1, the universe can evolve to the 
phantom phase < — 1, even though in the beginning, the universe is in the non-phantom 
(quintessence) phase with w h > —1. 

Furthermore, in order to demonstrate how the inhomogeneous modification in EoS in- 
fluences on the structure of the finite-time future singularities, as another example, we 
investigate the case that 

Zih(p) = f(p) + G(H) = -Ap a - BH 2 ? (3.137) 

= -Ap a -Bp 13 , (3.138) 

where B and j3 are constants, and by using Eq. (13.31) B is defined as 

B = Bl-j . (3.139) 

For > a, when p is large, the inhomogeneous part, i.e., the second term, of Eq. ( 13.138P 
becomes dominant as 

f ih -Bp? . (3.140) 
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On the other hand, for /3 < a, when p — > 0, the inhomogeneous part of Eq. (I3.138P becomes 
also dominant as in Eq. (I3.140p . 

We explicitly describe the cases in which the structure of the finite-time future singular- 
ities are changed due to the presence of the inhomogeneous term. First, we examine the 
limit that p diverges. For 1/2 < a < 1 with A > 0, if there only exists the homogeneous 
term without the in homogeneous term G(H) in Eq. (I3.130p . i.e., the case in Eq. (I3.77P 
in Sec. Ill D 1, there can appear the Type I singularity. However, in the presence of the 
inhomogeneous term as in Eq. ( 13 . 1 38j) . for f3 > 1(> a), in which the situation described 
by Eq. ( 13.1 40 p is realized because (3 > a, in the limit of p — > oo, \P\ also diverges because 
P — >■ —p — Bp 13 . Thus, the Type III singularity appears instead of the Type I singularity. 
For a = 1/2, if /? > 1(> a), the Type III singularity appears, whereas if (a <)l/2 < /? < 1, 
the Type I singularity occurs. For < a < 1/2, if /? > 1(> a) with B(B > 0), the Type 
III singularity appears, while if (a <)l/2 < j3 < 1 with B(B > 0), the Type I singularity 
occurs. 

Second, we consider the opposite limit that p tends to zero. In case of the homogeneous 
EoS in Eq. (13~77D . for < a < 1/2, in the limit of t ->■ t B , a -> a s , p ->■ and \P\ -> 0. The 
substitution of Eq. (13. 79j) into Eq. (13.781) leads to 



The exponent (a — 1) / (a — 1/2) in terms of \t — t s \ in Eq. (I3.14ip is not always an integer. 
Hence, there is the possibility that the higher derivatives of H diverge, even though a 
becomes finite. As a result, the Type IV singularity can appear. For a < 0, in the limit of 
t — > t s , a — > a s , p — 7-0 and \P\ — > oo. Thus, the Type II singularity can occur. 

We explain the case in presence of the inhomogeneous term as in Eq. (13.1381) . For a = 
1/2, if < < l/2(< a), in which the inhomogeneous term becomes dominant over the 
homogeneous one as in Eq. (I3.140p due to the relation /3 < a, the Type IV singularity 
appears, or if f3 < 0, the Type IV singularity occurs. For < a < 1/2, if j3 < 0(< a), the 
Type II singularity appears instead of the Type IV singularity. 




(3.141) 
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3. Implicit inhomogeneous EoS 



Next, in a general case, we suppose that in a proper limit, such as p being large or small, 
F{p) in Eq. (I3.132p is described by 

F{p)~Fp\ (3.142) 

where T and 7 are constants. For 7 ^ 1, we can integrate Eq. (I3.132p as 

?(t-t c )~£ , (3.143) 

1-7 

which is rewritten to 

p~ [(l- 7 )^(t-tc)] 1/(1 " 7) , (3.144) 

where t c is a constant of the integration. It follows from Eq. (13. 3p that the scale factor is 
given by 

a = a c exp { ± - - g ^ [(1 - 7 ) T it - Q] ) . (3.145) 

If 7 = 1, from Eq. (l3~T32l) we find 

p = pe ft , (3.146) 
where p is a constant. By substituting Eq. (I3.146P into Eq. (I3.3p . we obtain 



a = a c exp 



i^/L^/ 2 )* 

T V 3 



(3.147) 



We propose an implicit inhomogeneous EoS by generalizing the expression of F{p) in 
Eq. (|3~T32l) as 

T(P,p,H) = 0. (3.148) 

In order to understand the cosmological consequences of the implicit inhomogeneous EoS, 
we present the following example: 

(P + P?-C c p 2 (l-^ =0, (3.149) 

with C c (> 0) and H c (> 0) are positive constants. 

Plugging Eq. (13. 149[) into Eq. (13 .4p and using Eq. (I3.3p . we acquire 

H 2 = \C C H* (l - §) ■ (3-150) 
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We can integrate Eq. (I3.15ip as 

H= 9C*H c (t-t..)(t + -t)' (3 ' 151) 

with 

t±=t c ±—^—, (3.152) 
6u c n c 

where t c is a constant of integration. By combining Eq. (13.1491) with Eq. (13.1511) and sub- 
stituting Eq. ( 13 . 1 5 1 p into Eq. (I3.3p . we find 

3a 2 



p = - P 



c 



it - t c ) 



(3.153) 



p = ™— . (3.154) 

27C?H?k 2 {t - t^) 2 {t + - t) 2 



Thus, by using Eq. (I3.153p . we have 

w = - = -l-^-(t-t c ). (3.155) 

From Eq. ( 13 . 1 5 1 p . we see that if t- < t < t+, H > because t- < t c < t+. At t = 
t c = (t- + 1 + ) /2, H becomes the minimum of H = H c . On the other hand, in the limit of 
t — > t±, H — > oo. This may be interpreted that at t = i_, there exists a Big Bang singularity, 
whereas at t — t- , a Big Rip singularity appears. It is clearly seen from Eq. (13. 1 5 1 p that when 
£_ < t < t c , w > — 1 (the non-phantom (quintessence) phase, and when t c < t < t + , w < — 1 
(the phantom phase), namely, at t — t c , there can occur the crossing of phantom divide from 
the non-phantom phase to the phantom one. This is realized by an inhomogeneous term in 
the EoS. 

We present another example in which de Sitter universe is asymptotically realized. 



( p '+ py l+ ik {K ~ ' H) ^\~hT^ =0 ' 



16 „ „ x , fh c -H 



where t c , h c and hi are constants and h c > hi > 0. 



H = h c - hie~ t2/t " (3.157) 
p = ^(h c -hie- t2 / t *y , (3.158) 

P = --Ah c - hie-^Y- ^e"^ . (3.159) 
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It follows from Eq. ( 13.1571) that 



2/ll t j.2 1+1 

H = ^-e~ t/t ^. (3.160) 



c 



By using Eqs. O, f lXTBUj) and f[3TT3"Tj) . we obtain Eq. (E3D- For t < 0, H < 0, whereas for 
t > 0, H > 0. For i < 0, io e g > —1, while for t > 0, u> e g < —1- In the limit of t — > ±oo, 
the universe approaches to de Sitter one asymptotically, and hence there appears a Big Rip 
singularity nor a Big Bang singularity. 

We note that in principle, the implicit inhomogeneous EoS can be expressed in more 
general form by including higher derivatives of H as 

F(P t p,H,H,H,...) = 0. (3.161) 

There is a simple example 

2 • 3 (1 + w eff ) 2 
™ = w cS p -H - H . (3.162) 

From Eqs. (13. 3 p and (13 .4ft . we find 



P=^H\ p = WcSp -^H-^H 2 . (3.163) 



It follows from Eqs. (13. 162H and (13. 163H that Eq. (13.1621) is an identity. This implies that 
there exists a solution for any cosmology, provided that Eq. (13. 162ft is satisfied. 
We can also find another example 

P = w cS p-G 1 -^H + G 2 H\ (3.164) 

where G\ and G 2 are constants, and G\ (1 + w c q) > is supposed to be satisfied. For 
G*2 (1 + w c fi) > 0, we have a solution expressing an oscillating universe 

H = h c cos ut (3.165) 
a = a c exp ( — smut ) , (3.166) 



with a coefficient h c and a frequency u, given by 



hc = K \lT7T~r V' (3.167) 

3(1 + w eS ) 



3(1+ Weft) ,„ 1M > 

U= \l Gin 2 ■ (3 - 168) 
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On the other hand, for G 2 (1 + u> c fj) < 0, 



H = h c cosh ut , (3.169) 
a = a c exp ( sinhwt J , (3.170) 



where /i c is given by Eqs. (13. 169H and a frequency a) is described as 



Hence, the above investigations show that a number of models describing cosmology with 
an inhomogeneous EoS can be constructed. 



Finally, we mention that in Ref. [821 ]. cosmological density perturbations around the 
finite-time future singularities with the scale factor being finite have been examined. At the 
present stage, it seems that a number of models with a finite-time future singularity are not 
considered to be distinguishable with the ACDM model by using the observational test. As 
a recent investigation, the cosmological density perturbations in k-essence scenario has been 
investigated in Ref. 871 ] . 



IV. SCALAR FIELD THEORY AS DARK ENERGY OF THE UNIVERSE 

In this section, we explore scalar field theories in general relativity. 



A. Scalar field theories 



The action of scalar field theories in general relativity is given by 

S = J d 4 x^g- f A _ ^(<j>)g^(f>d v( f> - V((f>)\ + J d 4 xC M (g^, *m) , (4.1) 

where u(4>) is a function of the scalar field <fi and V(4>) is the potential of <fi. (i) For u(4>) = 
and V{4>) = A/k 2 with A being a cosmological constant, the action in Eq. ( 14. ip describes 
the ACDM model in Eq. ( 12. ip . (ii) For u(<j)) = +1, this action corresponds to the one 
for a quintessence model with a canonical kinetic term, (iii) For oj(4>) = —1, this action 
expresses a phantom model. We consider the case in which the scalar field is a spatially 
homogeneous one, i.e., it depends only on time t. 
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In the FLRW background (I2.3p . the Einstein equations are given by 

H 2 = y (p* + Pm) , (4.2) 

K 2 

H = -— (p^ + P^ + Pm + Pm) , (4.3) 

where p^ and P^ are the energy density and pressure of the scalar field 0, respectively, given 
by 

P4> = ^(0)0 2 + V{<f>) , (4.4) 
P = l -uj(<P)tf - V(<j>) . (4.5) 

In addition, pm and Pm are the energy density and pressure of matter, respectively. Equa- 
tions ( 14.4p and ( 14. 5ft give the way to express scalar dark energy as fluid dark energy, where 
scalar field equation becomes just the conservation law for such fluid description. 

Here, for clear understanding, by using Eqs. (I4.4p and ( I4.5p . we describe the explicit 
expression of u>de i n the case that there exists only single scalar field of dark energy, i.e., 
dark energy sufficiently dominates over matter, (i) For the ACDM model with oj(4>) = 
and V(<f>) = A/k 2 , wde = — 1- (h) For a quintessence model with oj{<p) = +1> 

if de = , 4.6 

2 + 2V(0) 1 ; 

it follows from which that —1 < u>de < —1/3. (iii) For a phantom model with co(4>) = — 1 

wve = - , 4.7 

02-2^(0) V 1 

which leads to wde < — 1. 

By using Eqs. (j4.2p - (l4.5p . we obtain 

2 



u(<j ) ) ( J ) 2 = -^H-(p M + P M ) , (U 

1 /_„o _--\ 1 



V{(j>) = - (3H 2 + H) - - (p M - P M ) • (4.9) 

If there is no coupling between matter and the scalar field <j), the continuity equations for 
matter and <f) are given by 

Pm + 3#( Pm + Pm) = 0. (4.10) 
P4> + 3H{ P(t> + P^ = 0. (4.11) 
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For a matter with a constant EoS u>m = Pm/pm, from Eq. (I4.10p we find 

Pm = PMca" 3(1+WM) , (4.12) 

where pm c is a constant. 

.t is the interesting case that ou(<j)) and V(<ft) are defined in terms of a single function 7(0) 



as 



W W = --2^ ~ C 1 + ^M)X c e- 3 ( 1+ ^) , (4.13) 

K. (2(2) 

V(0) = 1 (3/ 2 (0) + ^M) - ~ (1 - w u ) X c e~^ + ^ W) . (4.14) 



Here, 7((f>) is defined as 



W)^, (4.15) 



with X(0) and X c being an arbitrary twice differentiate function of </> and an integration 
constant, respectively. Thus, we can acquire the following solutions: 

= t, H = I(t). (4.16) 

For this solution, the equation of motion for <fi is derived from the variation of the action in 
Eq. (14. ip over as 

K0)0 + i^ + 3iM^+^M=O. (4.17) 
From these solutions in Eq. (I4.16p . we have 

f N i/[3(i+«m)] 

a(f) = a c e x W , a c = ^ . (4.18) 



In what follows, we consider the case in which these solutions are satisfied. 

For this case, in the FLRW background (12.3[1 . the effective EoS for the universe is given 
by Eq. dS3D with 

Pcff = P4> + Pm , (4.19) 

P cS = P (j) + P M . (4.20) 

If we define a new scalar field $ as 

$= y #vi^)i, (4.2i) 
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the action in Eq. (14. ip can be rewritten to the form 

" R 1 



d Xyj—g 



2 =f -g^d^d^ - V($) 



+ / d 4 x£ M (s>, *m 



(4.22) 



where the sign in front of the kinetic term depends on that of oj(4>). If the sign of oj(4>) is 
positive (negative), that of the kinetic term is — (+). In the non-phantom phase, the sign of 
the kinetic term is always — , and in the phantom one it is always +. In principle, it follows 
from Eq. (14.211) that <fi can be solved with respect to \ as = Hence, the potential 

V($) is given by V($) = V (<£($)). 



B. Equivalence between fluid descriptions and scalar field theories 

In this subsection, we show the equivalence between fluid descriptions and scalar field 
theories. We first take a fluid and then construct a scalar field theory with the same EoS 
as that in a fluid description. This process leads to constraints on a coefficient function of 
the kinetic term oj(4>) and the potential V((f>) of the scalar field <j) in the action in Eq. (14.11) . 
Through this procedure, we propose a way of expressing a fluid model as an explicit scalar 
field theory. In other words, we can obtain the explicit expressions of uj(4>) and V(0) in the 
corresponding scalar field theory for a fluid model. 

For simplicity, we suppose the dark energy dominated stage, so that we can neglect matter 
and therefore u>de ~ w c s, namely, p c ff ~ p — p<$> in Eq. ( 14. 19ft and p c s ~ p = p^ in Eq. (14.201) . 
In a fluid description, from Eq. ( 13.11) we find 

^eff = -=-l-M. (4.23) 
P P 

While, in a scalar field theory with the solutions in (I4.16p . we obtain 

«> = p + P=-f(p), (4.24) 
V= 1 -(p-P) = p + ^f, (4.25) 

where the second equalities follow from Eq. (13. 2p . Since Eq. ( 13. 3 p presents p = 3H 2 /k 2 , if 
H(= I it)) is given, we acquire the expression of p as a function of t(— 4>) as p = pit) = p(4>). 
Hence, by substituting this relation into Eqs. (14.241) and (14.251) . in principle we have u = ou((p) 
and V = V{4>). Thus, this procedure yields an explicit scalar field theory, which corresponds 
to an original fluid model. On the other hand, as the opposite direction, provided that we 
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have a scalar field theory described by oo(<f)) and V(4>) in the action in Eq. (14. ip . By using 
Eqs. (14.41) and (14. 5p and the solution <p = t and H = I(t) in Eq. (I4.16p . we get the explicit 
expression of w e s in Eq. (I4.23p . Hence, by combining this expression and p in Eq. (14. 4p and 
comparing the resultant expression with the representation of w c s in Eq. (13. ip . we acquire 
f(p) in the fluid description. Consequently, it can be interpreted that the considerations 
on both these directions imply the equivalence between the representation of a scalar field 
theory and the description of a fluid model. 

As concrete examples for a fluid description, we first consider the case of Eq. (I3.14p . By 
substituting Eq. fl3TT4D into Eqs. (j03)) - (05) . we have 

w cS = -1 - - (jfl - 1) . (4.26) 
P 

u = -p q + l, (4.27) 
o q - 1 

V = p+ P -^. (4.28) 

Second, for the case of Eq. (I3.17p . by combining Eq. (13.171) and Eqs. (I4.23j) ( |4.25j) we find 

w cS = --^. (4.29) 
A 



A 


pu+l 


f{ 




p 




2 





pu+l 

A 

p u+l 



(4.30) 
(4.31) 



As the third case described in Eq. fl3T47p . by plugging Eq. (E3ZD into Eqs. fl4T2"3l) - fl4T2"5]) we 
acquire 



2 / An / 3 , 

^ = - lT 5f-W^' (4 ' 32) 



2p An / 3 

"^f-irVsv (4 - 33) 



As the last example, we consider a model of Little Rip scenario given in Eq. (I3.70p . By using 
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Eqs. (13. 3 p and f)3.4j) . we obtain 

u , eff = -l-_^-exp(-et) , (4.35) 

w = -^fff L Rexp(#) , (4.36) 
y= g LR exp(et) [3ffL + (4 37) 

K z 

It is meaningful to remark that if the Hubble parameter H can be represented by t, by using 
the relations ( I4.13P and ( I4.14p and the solutions (I4.16P we can find the explicit expressions 
to = u}(4>) and V = V(4>). This is clearly demonstrated in the following Sec. IV C. 



C. Cosmological models 

In this subsection, we reconstruct scalar field theories corresponding to (i) the ACDM 
model, (ii) quintessence model, (iii) phantom model and (iv) unified scenario of inflation and 
late-time cosmic acceleration. In addition, (v) scalar field models with realizing the crossing 
the phantom divide is also considered including its stability issue. 

1. The ACDM model 

For the ACDM model, we have 

I(<f>) = H c , (4.38) 
where H c is a constant. From Eq. (I4.16p . we find 

H = H C , (4.39) 
a(t) = a c e Hct . (4.40) 

Moreover, by using Eqs. (I4.13P and (I4.14p . we obtain 

w (0) = - (1 + WM ) x ce -3(i+»M)xW ; (441) 

V(<f>) = ^$ ~ \ (1 - wm) X c e~^ 1+W ^ , (4.42) 

with 

!{4>) = H c (j). (4.43) 

In the ACDM model, we acquire 

w cG = -l. (4.44) 
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2. Quintessence model 



As an example of a quintessence model, we investigate the following model: 



(4.45) 



where H c and Hi(> 0) are constants and n(> 1) is a positive (constant) integer larger than 
unity. By using Eq. (I4.16p . we find 



H = H r + — 



a(t) = a c exp 



H r t 



(n - 1) t"" 1 

In addition, from Eqs. ( I4.13P and ( I4.14p we obtain 
2 nH x 



u{<j>) 



V{<t>) 



2 An+1 



K 2 (f) 

1 3 



,2 An+1 



kn-l 



l -{l-w M )T cQ -^ 1+w ™™ 



with 



I{(j>) = H c 



Hi 



in- I) 



kn-l 



Using Eq. (13. 9p . the effective EoS is written as 



2nH 1 t 



n-l 



(4.46) 
(4.47) 

(4.48) 
(4.49) 

(4.50) 
(4.51) 



3(# c ^ + #i) 2 ' 

From Eq. (14.51 [) . we see that w e s > — 1 because Hi > and n > 1. Thus, this model 
corresponds to a quintessence model. 



3. Phantom model 

As an example of a phantom model, we explore the following model: 

= ^ + < 452 > 

where H2 and H3 are constants and t s is the time when a Big Rip singularity appears. We 
examine the range < t < t s . From Eq. (I4.16p . we have 

a(t) = a c (t s - ty H2 e- (H3/ ' } . (4.54) 
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Furthermore, it follows from Eqs. (I4.13P and (I4.14p that 



K z 



Ho 



277, 



V(<f>) 



77 2 (377 2 + 1 



1 + Wm)1c (t B 



■>3(l+w M )H 2 e [3(l+w M )H 3 /<f>] 



(4.55) 



H3 f3H 3 

A3 



677 2 77 3 
(t s -<f>)c 



,3(l+w M )H 2 e [3{l+w M )H 3 /<f>] 



(4.56) 



(4.57) 



(4-^r 

Using Eq. (13.91) . the effective EoS is expressed as 

2[H 2 t 3 -2H 3 (t s -t) 2 ]t 

W c s = — 1 4~ 

3[H 2 t* + H 3 (t s -t)} 2 

It follows from Eq. (14.571) that when the time t approaches t s , the second term on the 
right-hand side (r.h.s.) of Eq. (14.571) becomes negative because H 2 t 3 — 2H 3 (t s — t) 2 > 
and therefore w c g < —1. Hence, there exists the phantom phase in this model and this 
consequence originates from the realization of 77 > 0, i.e., superacceleration. We note that 

1 /3 

when < t <C t s , for < t < (2H 3 t 2 /H 2 ) , w c g > —1 and therefore the non-phantom 
phase exists before the phantom phase appears. 



4- Unified scenario of inflation and late-time cosmic acceleration 



As an example of a unified scenario of inflation and the late-time acceleration of the 
universe, we study the following model: 



7(0) = h 



2 

1 t* 



(4.58) 



9) t 2 + 2 ' 

where h c is a constant, t s corresponds to the time when a Big Rip singularity appears, and 
ti is a time. From Eq. ( I4.16p . we find 



H = Ki 



a{t) = a c 



1 1 

+ 



t 2_ t 2 t 2 i + t 2 
ts + t ^ ' ' e (h£/ti) arctan(t/*i) 

u-t 



Moreover, by using Eqs. (I4.13P and (I4.14p . we obtain 



L>{<f>) 

v{4>) 



8 hi (tj + tg) [0 2 + (g - g) /2] 

(t? + 2 ) 2 (t s 2 - 2 ) 2 
1 /> 2 (t 2 + t s 2 ) 



« 2 (t 2 + 4>r (t s 2 - 2 )^ 



3/> 2 {tj + t 2 ) + 40 ( % 



(4.59) 
(4.60) 



{l + w u )Z c e- 3{1+Wu)m , (4.61) 



(4.62) 
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with 

ft s + 



Using Eq. ( 13. 9ft . the effective EoS is expressed as 



t(t-t+)(t 
2,hl t\ + 1\ 



w cs = -1 - TTTi , 2 , ^ , ( 4 -64) 



where i± = ±a/ (t^ — t\) /2. There exist two phantom phases: £_ < i < and £ > £ + , in 
which w c ff < — 1. On the other hand, there are also two non-phantom phases: — t s < t < t_ 
and < i < t + , in which w c fj > — 1. 

The history of the universe in this model can be interpreted as follows. The universe 
is created at t = —t s because the value of the scale factor a(t) in Eq. ( 14 . 6 j) becomes zero 
a(— t B ) = 0. During — t s < t < t_, there is the first non-phantom phase. The first phantom 
phase in t_ < t < corresponds to the inflationary stage. After inflation, the second 
non-phantom phase in < t < t + becomes the radiation/matter-dominated stages. Then, 
the second phantom phase in t > t + plays a role of the dark energy dominated stage, i.e., 
the late-time accelerated expansion of the universe. Finally, a Big Rip singularity occurs at 
t = t s . As a result, this model can present a unified scenario of inflation in the early universe 
and the late-time cosmic acceleration. Incidentally, phantom inflation has been studied in 
Ref. 



5. Scalar field models with the crossing the phantom divide 

The instability of a single scalar field theory with the crossing the phantom divide was 



examined in Ref. 



90(. In addition, the stability issue in a single scalar field theory as well 



as a two scalar field theory in which the crossing the phantom divide can be realized has 



recently been discussed in Ref. 9_1[. In this subsection, we first examine the stability of a 
single scalar field theory when the crossing the phantom divide occurs. The considerations 
for a two scalar field theory with the crossing the phantom divide are presented in Sec. VI 
A 2. 

We define the new variables Z = <fi and Y = I(<p)/H. With these variables, in the flat 
FLRW background the Friedmann equation ( 14. 2 h with Eq. ( 14 Ah and the equation of motion 
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for (I4.17P are rewritten to 

dZ 

In 



7,^(0) {z 2 - Y) 

dY = Ijrf) (1 - ZY) z 
dN H 2 



3{Z-Y) 



(4.65) 
(4.66) 



where I s< f>(<j>) = dl(cj))/d4> and I = d 2 I((j))/d(j) 2 , N is the number of e-folds and the scale 
factor is expressed as a = e N ~ No with iVo being the current value of N. For the solutions 
f l4.16p . we find Z = 1 and Y = 1. Hence, we examine the small perturbations SZ and 5Y 
around this solution as follows 



Z 
Y 



1 + SZ. 
1 + 6Y 



(4.67) 
(4.68) 



By using the solutions ( I4.16p . we see that these perturbations obey the equation 







' sxA 












\sx x 



(4.69) 



where A4i is a matrix, given by 



H 
HH 



H 

H 2 



H 
H 2 



(4.70) 



Thus, we can obtain the eigenvalues of the matrix .Mi as 

H 



m± 



HH 



H 

H 2 



3 ± 



( H 



HH 



H 

H 2 



4 



(4.71) 



The stability condition for the solutions ( I4.16P is that both of the eigenvalues m + and m„ 
are negative. When the crossing of the phantom divide occurs around H ~ 0, for H > 0, the 
transition from the non-phantom (quintessence) phase (H < 0) to the phantom one (H > 0) 
occurs, whereas for H < 0, the opposite direction transition from the phantom phase to the 
non-phantom one happens. In case of the expanding universe (H > 0), the term H / (^HH^j 
is negative. Thus, around the crossing of the phantom divide, from Eq. f!4.7ip we find 
m + ~ — H / (^HH^j > and m_ ~ 0. As a result, at the crossing time when H — m + 
becomes +oo and therefore the solution in Eq. ( 14. 16ft is unstable when the crossing of the 
phantom divide occurs. In other words, in a single scalar field theory the crossing of the 
phantom divide cannot be realized. 



55 



V. TACHYON SCALAR FIELD THEORY 



19 



In this section, we examine a tachyon scalar by following Ref. (9j. The effective 4- 
dimensional action for the tachyon field which is an unstable mode of D-branes [non- 
Bobomol'nyi-Prasad-Sommerfield (non-BPS) branes] is given by 



S 



d 4 xV(<j>) y/- det + d^d„<f>) , (5.1] 



where is a tachyon scalar field and V(4>) is a potential of 0. From the action in Eq. fl5.ll) . 
the energy-momentum tensor of <f) is derived as 

r<*> = v [ ff a ± - ^W^i + ^H. (5.2) 



In the flat FLRW background Q (with if = 0), p+ = -T o o(0) and = T- W are given by 

V(<f>) , s 

(5.3) 



P^-V^l-^. (5.4) 

We remark that through the use of Eqs. (15. 3p and (15.41) . this model may be written as a 
fluid. By using a/a = H 2 + H and plugging Eqs. fl5T3|) and (JED into Eqs. (I3T3]) and (I3~4l) . 
we find 

« = ^^2£L (i - ^) . (5.5, 

It follows from Eq. ( 15. 5 p that the condition for the accelerated expansion a > is given by 
<fi 2 < 2/3. Furthermore, the EoS for <fi is represented by 

W(t> = ^ = 2 - 1 . (5.6) 

P<t> 

Hence, from the above condition <fi 2 < 2/3 and Eq. (15. 6p we see that the possible range of the 
value of with realizing the cosmic acceleration is —1 < < —1/3, which corresponds 
to the non-phantom (quintessence) phase. 

We derive expressions of V(<p) and <fi in terms of H and H. By using H 2 = 



/c 2 /3) V((f))/Jl - <p 2 , which follows from Eq. (J3l| with Eq. (JE3D, and Eq. (|5T5]) . we have 
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H/H 2 



(3/2) <f) 2 . The combination of these equations leads to 



v(<f>) 




(5.7) 



(5.8) 



We suppose that the scale factor is expressed by a power-law expansion as a(t) oc t p with 
p > 1 being a constant larger than unity. By combining this relation with Eqs. ( 15. 71) and 
(15.8jl . we obtain 



2 

3p 



2 1 




(5.9) 
(5.10) 



where we have taken the integration constant as zero. As a result, we acquire an inverse 
square power-law tachyon potential V(<ft) oc -2 . 



In the open string theory, the form of a tachyon potential V(<j)) is given by 92 j 

Vo 



V{<j>) 



cosh 



(5.11) 



with </>o = \/2 for the non-BPS D-branes in the superstring and 0o = 2 for the bosonic 
string. This form has a ground state in the limit — > oo. Here, Vq and 0o are constants and 
V(<f) = 0o ) = Vo- Moreover, when a tachyon potential appears as the excitation of massive 
scalar fields on the anti D-branes, V(<f>) is given by 93] 



V((f>) = V e m2 ^ /2 , 

where m is the mass of and there exists a minimum of V(<fr) at = 0. 

We mention that a tachyon scalar field can be generalized to so-called k-essence 



(5.12) 



94| , which 



is a scalar field with non-canonical kinetic terms (for its application to inflation, so-called 



k-inflation, see 



95 



96[), and for a unified scenario between inflation and late-time cosmic 



acceleration in the framework of k-essence model, see, e.g., |97|). The action of k-essence is 
described by 



S 



d^xyj—g 



R 



2k 



_-2 



P(0,X) 



d A xC M (fiV' *m) 



(5.13) 
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where P((f), X) is a function of a scalar field and its kinetic term X = — (1/2) g^ v d^,(j)d v (j). 
It is known that the accelerated expansion can be realized by the kinetic energy even without 
the scalar field potential. An effective 4- dimensional Lagrangian describing a tachyon field 



is given by 



98] 



P = -V((f))Vl - 2X. (5.14) 

Finally, we discuss a stability issue. We now consider the action in Eq. ( I5.13p . In this 
case, the energy density and pressure of are given by = 2XP X and P^ = P, respectively, 
where Fx = dP/dX. The EoS for is written as 

W * = 2XP X -P- (5 ' 15) 
If the relation \2XP jX \ <C P is realized, can be close to —1. 

We investigate the stability conditions for k-essence in the ultra-violet (UV) regime. 
In the Minkowski background, we write as <fr(t,x) = <frb(t) + 5<ft(t,x), where (pb(t) is 
the background part and 5<p {t, x) is the perturbed one. By deriving the Lagrangian and 



Hamiltonian for 5<j)(t,x), the second-order Hamiltonian is obtained as 99] 

m = 6 M +6 (^ +6 (^!, (5 , 6) 

with 

^ = (P,x + 2XP xx )>0, (5.17) 

^2 = P,x>0, (5.18) 

& = -i^>0, (5.19) 

where P jXX = d 2 P/dX 2 and P^ = d 2 P/d<f) 2 . The stability conditions are represented by 
the positivity of each three terms on the r.h.s. of Eq. (I5.16p . The propagation speed of is 



defined as 



= 3»£ = 6 (5.20) 

As long as the stability conditions £i > in Eq. (15.171) and £2 > in Eq. (15.181) are satisfied, 
Cg > . In addition, the condition that the propagation speed should be sub-luminal is given 



by [96J 



Pxx>0. (5.21) 
We note tha t the finite-time future singularities in tachyon cosmology have also been exam- 



ined in Ref. 



100] 
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VI. MULTIPLE SCALAR FIELD THEORIES 

In this section, we describe multiple scalar field theories. 



A. Two scalar field theories 



To begin with, in this subsection we investigate two scalar field theories. First, we explain 
the standard type of two scalar field theories. Next, we discuss a new type of two scalar 



field theories with rea 
constructed in Ref. [91] 



izing the crossing of the phantom divide, which has recently been 



27 



101 



1021 ] . The action of 



1. Standard two scalar field theories 

First, we explore the standard type two scalar field theories 
two scalar field theories in general relativity is given by 

S = J dS^Tg f A _ 1 W (0)^<R0 - 1 -a( X )g >lu d,xd,X ~ X )j , (6.1) 

where <r(x) is a functions of the scalar field x an d V(<f>, x) is the potential term of and x- 
Here, we concentrate on the scalar-field part of the action and do not take into account the 
matter part of it. If there does not exist the second scalar field x ; the action in Eq. ( 16. II) is 
the same as the action in Eq. (14. ip without the matter part of it. 

In the FLRW background (12. 3p . the Einstein equations are given by 

,.2 



H 2 = - Pt 



H 



-y(Pt + Pt) 



(6.2) 
(6.3) 



where p t and P t are the total energy density and pressure of the two scalar fields and Xi 
respectively, given by 



p t = \uj(<P)4> 2 + \a(x)x 2 + V((f), x) 



2 



2 

2 + ^(x)x 2 -v(^x) 



(6.4) 
(6.5) 
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Moreover, equations of motion for the scalar fields are given by 

+ 5^3jT^ + 3^(0)0+^^ = 0, (6.6) 
a(M + + SHa( X)X + = . (6 . 7) 

Since it is possible to redefine the scalar fields through a convenient transformation, we 
take = x — t- If a solution H(t) = I{t) is given, where I(t) is a function of t, by plugging 
Eq. (16. 2p into Eqs. ( 16. 6 p and (16. 7p . we obtain 

W(0) = ^^^' = --f-g£— (6 ' 8) 

Here, by taking <p = \ — t into consideration, we can interpret /(</>, x) a s I{t,t) = I(t). 
Thus, we define /(0, x) as 

I(<t>, X ) = -y ^ u;(0)# + y ct(x)^ ■ (6.9) 

We can represent the potential term V(<f), x) by 

nM .^ M+ «S^ + «S^L). (6,6) 



Furthermore, Eq. ( 16. 3 p is rewritten to 

2 dl(t) 



w(t) + <r(t). (6.11) 



k 2 dt 

In addition, we can describe a part of coefficient of the kinetic terms, which is a function of 
a scalar field, Lo(<f) and <r(x), as 

^-1(^4 (6 - i2) 

^(x) = 4^w> ( e - 13 ) 

where #(0) is an arbitrary function of 0. Thus, by su bstit uting Eq. (16. 9p with Eqs. ( 16 . 1 2[) 
and (16.130 into Eq. (16.101) . we acquire (for details, see 102| ) 



V{h X) = ^ (3/ 2 (0, X) + ^ + m - ~9(X)) • (6.14) 

In order to demonstrate an example of two scalar field theories, we examine the following 
model in Eq. (|4.52p : 

/ \ Ho H% 
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In this case, by using Eqs. ( I6.12p and f 16 . 1 5 j) we find that oj(4>) is expressed as 



Ho 



2H* 



(6.16) 



l(t s -t) 2 

and <t(x) is described as in Eqs. ( I6.13p . In addition, it follows from Eqs. (16. 9p . (16.101) and 
(16.15P that I(<f),x) and V{4>, x) are represented as 

#2 



KM 
v(<j>, x ) 



t s - 
i 



§ + / ^<M0) 



3/ 2 (0,x) 



#9 



2i/ 3 



+ <i(<t>) -g(x) 



(6.17) 
(6.18) 



(t s - ty 

We investigate the stability of the system described by the solution in Eq. ( 16.151) . We 
define the following variables: 



Xa 



H 



(6.19) 



By using the variables in Eq. (I6.19p . Eqs. (16. 2p . (I6.6P and ( 16. 7p are rewritten to 

dX^ 1 du{(j))/d(j) 



dN 

dXx 
dN 
dY 



2H u{<f)) 

1 da(x)/dx 
2H a{ X ) 



(X*-1)-3(X+-Y) 
(X 2 X -1)-3(X X -Y) 



K 



[uifiXi {YX^ - 1) + a( X )X x (YX X - 1)] 



(6.20) 
(6.21) 
(6.22) 



dN 2H 2 

where we have used d/dN = (1/H) d/dt. We analyze the perturbations \8X^\ <C 1, \8X X \ <C 
1 and \SY\ < 1 around (X^, X x , Y) = (1, 1, 1) as 



X* = 1 + dX^ , X x = 1 + 5X X , Y = 1 + SY 



(6.23) 



The perturbations satisfy the equation 



_d_ 

dN 

where M is the matrix, defined by 



( sx^ 

5X X 
\ 6Y , 



M 



^ 5X^ 

5X X 



(6.24) 



M = 



( dm((j})/d(j} 
Hu){4>) 





V 



H 



2H 2 





dv{x)/ d X 
Ha{ X ) 

k .2ct(x) 



3 

3 3 

_2w(<W+o-(x) 



2H 2 



(6.25) 
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The characteristic equation for M is given by 



det |M - \E\ = 

/da#V# \ f da(x)/d X 

V HuW J V Ha( X ) 

3k 2 o;(0) / da(x)/dx x \ 3k 2 ct(0) / du(<f>)/d<f) 



3 + A k 



2 ui{ (j>) + (r(x) 
2H 2 



A 



3 + A + 



where A denotes an eigenvalue of M and E is a unit matrix. We impose the following 
conditions 

w(0)^O, a(x)^0, (6.27) 

so that the eigenvalues can be finite without diverging. As a result, if the conditions 
Eq. (16.271) are satisfied, the solution in Eq. (I6.15P does not have infinite instability at the 
crossing of the phantom divide from the non-phantom phase to the phantom one. For an 
illustration, e.g., we take, g(t) = g c /t 3 with g c being a constant and satisfying g c > 2H 3 . 
From Eq. ( 16.171) . we find 



H 2 9c~2H 3 | g 



2x 5 



(6.28) 



Moreover, it follows from Eqs. (16 . 12[) and (16 . 1 3[) that 

H 2 



u{<j>) = - 

9. 



2 

K 2 X 3 



(6.29) 
(6.30) 



Furthermore, by using Eq. ( 16 . 1 8[) we obtain 



1 



v{^x) = -, 



3/ 2 (0,x) 



Ho 



9c - 2H 3 



9c 



(ts - <f>) 

Examples of two scal ar fie ld theories are an oscillating quintom model 
with two scalar fields 



(6.31) 



104j | in the framework of general relativity (see also 



103 1 or a q uintom 



105 



iod |). 



2. New type of two scalar field theories 

A new type of two scalar field theories is given by 



S = J d 4 x^9 ( A + X + u(<P)X 2 -U + V ( X )U 2 - V{<j>, X )j + J d 4 xC M (g, u , * M ) 



(6.32) 
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where X = — (1/2) g^d^d^, U = — (1/2) g^d^x^vX- We suppose that there is no direct 
interaction between the scalar fields and matters, so that the continuity equation of matter 
can be satisfied p' M (N) + 3 (pu(N) + Pm(N)) = 0, where a prime denotes the partial deriva- 
tive with respect to N of d/dN. In the flat FLRW background, for the action in Eq. (I6.32p 
the gravitational field equations are given by 

X + u{<P)X 2 -U + V ( X )U 2 - V(<f>, x) = ~\ (2HH' + 3E 2 ) +p M + ^, (6.33) 

X + 2u(<P)X 2 -U + 2 V ( X )U 2 = --.EE' + p M + ^ . (6.34) 

k z 6 

On the other hand, the Friedmann equation is also represented as E 2 (N) = 
(k 2 /3) (p B (N) + pm(N)), where p s is the energy density of the scalar fields, namely, we have 
expressed the right-hand side by dividing the energy density into the contributions from the 
scalar fields and matter. By using this expression, the gravitational field equations f!6.33j) 
and (I6.34p are rewritten to 

V(<j>, X)-\{X-U)= p s (N) + Mp. , (6.35) 
^)X 2 + V ( X )U 2 + 1 -(X-U) = . (6.36) 

Furthermore, the equation of motion for the scalar fields are given by 

V^cf), x) + 3u^((f))X 2 = —E 2 (1 + 6ooX) f - [EE' (1 + 6ooX) + 3E 2 + QE 2 ooX] </>' , (6.37) 
V x (4>, X ) + 3 V , X ( X )U 2 = E 2 (1 - 6 V U) X " + [EE' (1 - 6r)U) + 3E 2 - 6E 2 V U] X ' , (6.38) 

where the subscription " denotes a partial derivative with respect to 0, e.g., V^(<^, x) = 
dV((fi,x)/d(j) and V x ((j), x) = dV((j),x)/dx- Thus, we acquire the following solutions 

(f) = x = mN, (6.39) 
H 2 (N) = j(J(N)+p M (N)) , (6.40) 
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provided that V(4>, %), w(0) and rj(x) satisfy the equations 

V (rhN, rhN) = J(N) + ^-ffi , (6.41) 

u (fhN) + r,(fhN) = - ^^ , (6.42) 

V^(fhN,mN) + -fh 4 H 4 (N)cj^(mN) = -rh (3H 2 (N) + HH'(N)) 

-3rh 2 H 2 (N)u(rhN) (H 2 (N) + H H' '(N)) , (6.43) 
V iX (rhN, rhN) + |jm 4 # 4 (A07/, x (mA0 = rh (3H 2 (N) + HH'(N)) 

- 3rh 2 H 2 (N)r](rhN) (H 2 (N) + HH'(N)) , (6.44) 

where J(N) is an arbitrary function and rh is a constant. We derive the explicit forms of 
v(x) and V((j),x) m the action in Eq. ( 16.321) so that these should obey Eqs. ( 16.411) - 
(I6.44p . As an example, by using an arbitrary function a and J we can express ou((f)) and 
rj(x) as 

a ^=- «imm % (6 ' 45) 

i(x) = -a(x) ■ (6-«) 

In addition, we define another function J(0, x) as 



^(0,X) = -m 



d(j)' (3rh 2 uo{(j)')H A {(i)' /rh) + 2H 2 ((f) / /rh)) 



+ J d X ' {3m 2 r ] ( X ')H\x'/m) - 2H 2 ( X '/rh)) 



(6.47) 



It follows from Eq. flQZD that J (rhN, rhN) = J(N). With Eq. (IBIETj) . we determine the 
form of V(<f>, x) as 

W,x) = ~m (| d<f>>H\cf>>/m) - J dx'H 2 (x'/m?j + J(cf>, X ) + f + . 

(6.48) 

Thus, through this procedure, it is possible to reconstruct any expanding history of the 
universe by using two arbitrary functions J of <$' /rh or x'/m and a of or x- 

Next, we explore the stability of the solutions ( 16.391) and (I6.40p . There are two ap- 
proaches to study the stability. One is the perturbative analysis. Another is to examine the 
sound speed of the scalar fields. First, we investigate the first order perturbations from the 
solutions, given by 

= 4 >O + 64 > (N), x = Xo + Sx(N), <P = <p + 5x(N), X = Xo + Sy(N) , (6.49) 
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where we have defined Sx(N) = S<f>(N) and 5y(N) = 5x(N). From the Friedmann equation, 
there exists the following constraint between the perturbations 

K 2 Sp s 



5H 



6H 



0(1 + 3u,0 2 )6x+[ + V >4 ) 5<f> + X ("I + 3r/x 2 ) Sy + (^f^ + V X )5 X - (6.50) 



Hence, in this system there are four degrees of freedom (</>, x, <fi, X)- 
Eqs. ( 105D -( jOHD and fhlD and using Eq. (IBToTJ]) . we obtain 



By combining 



where M.\ is a matrix, given by 







/ 5<A 




= .M 2 




5x 




5x 


\*y) 









H- 1 








i 



\ 

r-l 



M? 1 ' M? -Mf } -Mf } 



(6.51) 



(6.52) 
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with 



M 31) 



2HH' - 2H' 2 + m 3 ujH 3 H' + m 2 K 2 H 2 (1 + m 2 uH 2 ) 



H(l + 3m 2 uH 2 ) 



+ 



H 



AHH' + 3H' 2 + HH" 



m 2 K 2 



(1 + m 2 ujH 2 ) (H 2 + HH') 



m 2 K 2 (1 + m 2 uH 2 ] 
H(l + 3m 2 uH 2 ) 



- (H 2 + HH') (-1 + 3rh 2 r]H 2 ) - H 2 



— 2 2 

m k 



(l + m 2 uH 2 ) 



3m 2 H (2uH' + fhujH) 



3m 2 ujH 2 



M 



(34) 



m 2 K 2 1 + m 2 uiH 2 



2 1 + 3m 2 uH 2 



(-1 + 3m 2 T]H 2 ) 



M 41) 



a4 42) 



171 K ( — 1 + m 2 7]H 

H (-1 + 3m 2 r]H 2 ) 



2\ V 



(H 2 + M') (1 + 3m 2 wtf 2 ) 



2M' - 2H' 2 - 3m s r] !X H s H' + m 2 n 2 H 2 (-1 + m 2 ^iJ 



if 2 

r2\ 



(6.53) 
(6.54) 

(6.55) 
(6.56) 
(6.57) 



M 43) 



1 

+ H 

rh 2 K 2 



H (-1 + 3m 2 r]H 2 ) 

— 2 2 

AHH' + 3H' 2 + HH" + (-1 

2 v 

1 + rh 2 rjH 2 



m 



■rjH 2 ) (H 2 + HH') 



-1 + 3m 2 rjH 2 



(l + 3rh 2 ujH 2 ) 



M 



(44) 



m 2 K 2 



(-l+rh 2 r)H 2 ) 



3m 2 H (2r]H' + m?; iX i7) 



, (6.58) 
(6.59) 
(6.60) 



2 \ ■ - i -l + 3m 2 r]H 2 

The condition for the solutions ( I6.39P and f)6.40p to be stable is that the real part 
of all the eigenvalues of the matrix Ai 2 f !6.52j) should be negative. The characteristic 
equation for M 2 is given by det \M 2 - XE\ = A 4 + A^ 3 + A 2 X 2 + A 3 X + A A = 0, 



where A\ 

A 3 = 



-M 



(33) 
2 



(44) 



A 2 



-Mr-^2 



M 31) .M 2 42) 



M 32) a4 41) ) 



H 



^M 2 
H~ l , and Aj = 
The solutions of this equation are given by A = 
Bi ± \/—Bi — 2E ± A^E — £> 3 j /2 — Ai/4, where E satisfies the cubic equa- 
tion S 3 - (B 1 /2)E 2 - B 3 E + BxB 3 /2 - B 2 /8 = 0. Here, B 1 = -3A 2 /8 + A 2 , B 2 = 
A 3 /8 - A ± A 2 /2 + A 3 , and B 3 = -3A 4 /256 + A 2 A 2 /16 - A x A 3 jA + A 4 . The solution of 
this cubic equation is given by E = (-V - VP 2 + Q 3 ) 1/3 + (-V + VV 2 + Q 3 ) 1/3 + Bi/6, 
with P = -5 3 /216 + 51^3/6 - B 2 /16 and Q = -52/36 - B 2 /3. We explore the real so- 
lution of the cubic equation so that the maximum of the real parts of A, which is described 
by ReA max , can be negative. As a result, H is expressed as follows, (i) For Q < and 



31 



V 2 < \Q 



and V > 0, E 



2-v^Qjcos [(l/3)arccos (-V\Q\- 3/2 )] + B ± /6. (ii) For Q < 0, V 2 > \Q 3 



(V + VV 2 + Q 3 ) 1/3 - (V - VV 2 + Q 3 ) 1/ " + Si/6, (iii) For Q < 0, 
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1/3 



V 2 > \Q 3 \ and V < 0, 5 = (-V + VV 2 + Q 3 ) 1/3 + (-P - \JV 2 + Q 3 ) 1 ^ + Si/6, (iv) For 



,V3 



Q > 0, 3 = - (V + V^ 2 + Q 3 ) 1/3 + (-"P + VV 2 + Q*y /0 + Si/6. 
Moreover, ReA max is represented as follows. If 3 2 > S 3 , 



1/3 



ReA r 



i ( V / 2S^+ ^-23 - ^ + 4v^^^) - ^ 



for 3 > — and - 23 - S x + 4/ > 
- 2 



v 2S^-^i<0 
2 v 4 ~ 



for 3 > — - and 
~ 2 



23 - Si + 4/ < 



A S 

22 - Si + 4^3 - S 3 - for 3 < -y and - 23 - Si + 4/ > 

A. B 
— -± for 5 < -y and - 23 - Si + 4/ < 



< 0. 



(6.61) 



On the other hand, if 3 2 < S3, 



Re A 




Si + J(2E + B 1 



-23 - Si + d (23 + Si) z + 16 (3 - S 3 



for 3 < — - 
2 



< 0. 



16 



4 



4 



(6.62) 



Thus, if the conditions described by Eqs. (16.611) and 1 16. 62ft are satisfied, the solutions ( 16.391) 
and (16.40P can be stable. Since these solutions (16. 39ft and (16.401) correspond to those (14. 16ft 



for a single scalar field theory discussed in Sec. IV, through the investigations in Sec. IV 
C 5, the stability of these solutions implies that the crossing of the phantom divide can be 
realized in the new type of two scalar field theories described by the action in Eq. (I6.32p . 

In principle, the stability condition from the perturbative analysis can yield constraints 
on the forms of, e.g., a'(rhN), a(rhN), J{N) and 0J^(fhN). However, it is difficult to derive 
the explicit analytical expressions of such a constraint on a'(rhN) or J(N). On the other 
hand, the stability condition obtained by the sound speed of the scalar fields can present 
the analytical representations of constraints on a(mN). Therefore, we explore the sound 
speed of the scalar fields. Its square has to be positive for the stability of the universe. 
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The sound speed c 2 j, where j = 0,x, of the scalar fields are defined as = P^x I ' P<t>,x = 
(1 + 2cuX) j (1 + 6wX) and c s 2 x = P X) x/p x ,x = (-1 + 2r]X) / (-1 + 6r/X). Hence, the sta- 
bility condition is expressed as < c 2 j(< 1). For the solutions ( I6.39P and (16.40 p . this 
condition can lead to constraints on the function a. as a(0 = ffiN) > J' / (3m 4 // 4 ) or a(<f) = 
rfiN) < J 1 / (3m 4 # 4 ) - 1/ (m 2 H 2 ), and a{<j> = rhN) > or a(<f> = rhN) < -1/ (m 2 H 2 ). 

It is interesting to mention that in Sec. IV C 5, we have shown that the crossing of 
the phantom divide cannot occur in a single scalar field theory represented by the action 
in Eq. (14.11) because of the instability of the solutions (14.161) . whereas in the new type of 
two scalar field theories whose action is given by Eq. (16.321) . the crossing of the phantom 
divide can happen due to the stability of the solutions (I6.39P and (16.401) . This result can 
be a proposal of a clue for the searches on the non-equivalence of dark energy models on a 
theoretical level. 



B. n(> 2) scalar field theories 

We generalize the investigations for two scalar field theories in Sec. VI B. The action of 
n scalar field theories in general relativity is given by 

S = j d 4 xV=9 (j^ - \ "iiWW^i - V(<f>i, ■ ■ ■ , <f> n ) j , (6.63) 

where u>i(4>i) is a function of a scalar field 0j and V(0i, 02, • • • , 4>n) is the potential term of n 
scalar fields <f)\, <p2, ■ • • , 4> n - For n = 2, i.e., two scalar field theories, this action in Eq. (I6.63P 
with (pi = (f>, (f>2 = — w (0) an d ^2(^2) = cr(x) is equivalent to that in Eq. (16. ip . In 

the FLRW background (12. 3p . the Einstein equations are given by Eqs. (16. 2p and (16. 3 p with 
p t and P t being 



Pt 

i=l 



1 n 

- "Mi* + V{<f> u 2 , • • • , 0„) , (6.64) 

n 

P t = J2"i(<t>i)<i>l (6-65) 

i=l 

We can apply the same procedure executed in the case of two scalar field theories for the 
n scalar field ones. From Eq. (16. lip , we analogously find 

/ x 2 dl(t) 

i=l 
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Thus, from Eqs. (I6.12p and (I6.14p . we obtain 

(a\ 2 91{<j>w" ,<t>n) fRa ^ 
^) = -- 2 ^ , (6.67) 

V{K ■■■An)- ^ ( 3J*(fc, ■■■An) + f t > ( 6 - 68 ) 

where I(<j>i, ■ ■ ■ ,4> n ) can be regarded as /(£,••■ ,t) = I(t). Furthermore, we acquire the 
following solutions: 

<j>i = t, H(t) = I(t). (6.69) 
In addition, from Eq. (I6.66P we take 

wi(^i) = - Jj + sfa(^i) + • • " + , (6.70) 

2 

= —9j(0j) , 3 = 2, 3, . . . , n . (6.71) 

Hence, there exist n — 1 arbitrary functions g2, - ■ ■ ,g n which can reproduce the solution 
( I6.69p . and therefore the reconstruction can be executed through the procedure explained 
above. 

Here, we explicitly demonstrate the equivalence of multiple scalar field theories to fluid 
descriptions. By combining the consequences in Eqs. f l6.67p - fl6.69p and the investigations 
in Eqs. f l4.23p -f l4.25p in terms of the equivalence between fluid descriptions and scalar field 
theories, we obtain 

n 

/(p) = -X>(t), (6.72) 

i=l 

1 

Pt = - P t - f( P ) = -v(t) + - J>iW , (6.73) 

i=l 

Pt = V-^ = V(t) + ±j^u i (t), (6.74) 

w Pt_ -2V-(t) + EIU^) f675 , 

where in deriving the second equalities in Eqs. (I6.73P and f)6.74p we have used Eq. (16.721) . 
and Eq. (I6.75P follows from Eqs. (16.731) and (16.741) . Accordingly, these results imply that 
multiple scalar field theories can also be represented by a fluid description, similarly to that 
for a scalar field theory as shown in Sec. IV B. 
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VII. HOLOGRAPHIC DARK ENERGY 



In this section, we study holographic dark energy scenario and its generalization by 



following Ref. 



27] through the analogy with anti de Sitter (AdS)/ conformal field theory 



(CFT) correspondence. In particular, we investigate the case in which the infrared (IR) cut- 
off scale is represented by a combination of the particle and future horizons, the time when 
a Big Rip singularity appears (namely, the life time of the universe), the Hubble parameter 
and the length scale coming from the cosmological constant . 



A. Model of holographic dark energy 



First, we make an overview for a model of holographic dark energy 25]. The energy 
density of holographic dark energy is proposed as 

with Ch being a numerical constant. Here, Lh is the IR cut-off scale with a dimension of 
length. At the dark energy dominated stage, we suppose that the contribution of matter is 
negligible. From Eq. (I3.3p . we have the following Friedmann equation: 3H 2 /k 2 = p^. By 
combining this equation and Eq. (17. ip . we find 

H = ^ , (7.2) 

where Ch(> 0) is assumed to be positive in order to describe the expanding universe. 

It is necessary for us to discuss how to take the IR cut-off Lh because if the IR cut-off is 
identified with the Hubble parameter, the cosmic acceleration cannot be realized. We define 
the particle L p h and future Lfh horizons as 

/"* dt' 

^^Lwr (7 - 3) 

f°° dt' 

= a(t) / — - . (7.4) 
Jt a{t) 

In the flat FLRW space-time in Eq. (12. 3p with K = 0, provided that the IR cut-off scale Lh 
is identified with the particle horizon L p h or the future horizon Lfh, we obtain 

d { C h \ .1 



±-, (7.5) 
dt \ an I a 
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where the plus + (minus — ) sign corresponds to L p h We can solve Eq. (I7.5P as 

a = a c t hc , (7.6) 

Thus, for Lh = L^, an accelerated expansion of the universe can be realized, because the 
power index h c in terms of t in Eq. (17. 7p is larger than unity and hence power-law inflation 
can occur. On the other hand, for = L p h, the universe will shrink due to h c < 0. 

We suppose that the theory is invariant under the change of the time direction as t — > —t. 
In addition, by shifting the origin of time appropriately, we have the following expression 
for a instead of Eq. (17.61) : 

a = a c (t B - t) hc . (7.8) 

Thus, for h c < 0, a Big Rip singularity will appear at t — t s because a diverges at that time. 
If we change the direction of time, the particle horizon becomes like a future one as 

r ,s [ U dt' . , f 00 da 

From Eq. (EU) with Eq. (ESD or Eq. (O]) . we obtain 

«- = -l + 3£. (7-10) 

We remark that when we take as the future horizon in Eq. (I7.4p . we can acquire 
a de Sitter solution 

a = a c e HdS \ (7.11) 

#ds = y, (7.12) 

where I is a constant denoting a length scale and hence H^s is a constant Hubble parameter 
describing de Sitter space. It follows from Eq. (17. ip and = I that ph = 3C^/ (k 2 1 2 ). For 
Ch = 1, by using Eq. (I7.12p we identically find the Friedmann equation (I3.3P as 3H 2 /n 2 = p h , 
whereas for C*h 7^ 1, the de Sitter solution (I7.12p cannot be satisfied. When we choose as 
the particle horizon, there does not exist the de Sitter solution because the particle horizon 
in Eq. (17. 3p varies in time as L p h = (l — e*/ z ) // and not a constant. 
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B. Generalized holographic dark energy 



Next, we explain generalized holographic dark energy |26|, |27j. It has been pointed out 



in Ref. 



1071 ] that if Lh is taken as L p h, the EoS vanishes because L p h behaves as being 



proportional to a 3//2 (t), although the value of the energy density is compatible with the 
observations. Therefore, we examine the generalization of holographic dark energy. In more 
general, Lh could be represented as a function of both L p h and L^,. Provided that the life 
time of the universe is finite, then t s can correspond to an IR cut-off, and Lfh in Eq. (I7.4p 
is not well-defined because of the finiteness of the cosmic time t. Thus, the future horizon 
may be re-defined by 

dtf , , f 00 da 



fh 



a(t) 



<t') a(t) Jo 



(7.13) 



as in (17. 9p . By analogy with AdS/CFT correspondence, we suppose that L h may be described 
by 

Lh = Lh(L p h, Lfh, t s ) , (7-14) 
as long as L p h, Lfh and t s are finite, because there exist a lot of possible choices for the 



IR cut-off 



1081 ] . We note that holographic dark energy fr om t he Ricci scalar curvature 



10j}. 



Moreover, there exist a 



the so-called Ricci dark energy has been explored in Ref. 
number of s tudie s on holographic dark energy in theoretical aspects as well as observational 
constraints H0(. A concrete example is given by 26] 

Lh _ X 

a = 

l/h c 



X 




z-B(p,q) 



L p h + Lfh 

where h c > 0, and B(p, q) is a beta-function defined by 

POO 

B(p,q) = I 

p = n 



(7.15) 
(7.16) 



ff>- 



p+q 



\l + t) 
h c , q = l-h c . 



It follows from Eq. (IT. 151) that the solution is given by 



H = h c - + 



ts-t 

lie 



t B -t 



(7.17) 
(7.18) 

(7.19) 
(7.20) 
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Furthermore, since we have 



r ts dt ( t \ hc 

La + L m = a(t) / - = t s [ B(p, q) , (7.21) 

Jo a V^s — tj 

by substituting Eq. ( 17.2 lj) into Eq. (17.15p . we acquire 

where the second equality follows from Eq. (17.201) . Clearly, Eq. (17.221) is equivalent to 
Eq. (O. 

We mention that for the solution (17.201) . from Eq. (13.91) we find w c q = — 1 + 
2(t s — 2t)/ (3h c t s ). When t — > 0, w c s — > —1 + 2/ (3h c ) > —1, i.e., the universe is in 
the non-phantom phase. At t — t s /2, w c g = —1. After that, in the limit of t — > t s , 
w c ff — > —1 — 2/ (3/i c ) < —1, i.e., the universe is in the phantom phase. Consequently, it can 
occur the crossing of the phantom divide. 

We also investigate the case that there exists matter with its EoS being w m = P m / p m . 
In what follows, with w m we define h c as h c = (2/ 3) / (1 + w m ). We assume the existence of 



an interaction between holographic matters 111] . The equation for matter corresponding 
to the continuity equation is given by 

p m + 3H ( Pm + P m ) = , (7.23) 



where p c is a constant. We also suppose 

h (l ^V' /2 X (7 24) 

In this case, the Friedmann equation (13 .3p becomes 3H 2 /k 2 = Ph + Pm- Thus, we obtain the 
solution (17.201) as well as p m , given by 

) . (7.25) 
t t B -tJ 

By using the Friedmann equation 3H 2 /k 2 = p h + p m , Eqs. (17.201) and (I7.25p . we find that 
the ratio of the energy density of holographic dark energy in Eq. ( 17. II) to that of matter 
becomes constant and it is given by 

P h 3h c f 1 _ k^A (? 2g) 



Pm K 2 p c \ 3k 
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This consequence may be a resolution of coincidence problem between the current energy 
density of dark energy and that of dark matter. 

In Ref . 11 2j , a naive model of such an interaction scenario between dark energy and dark 



matter in (I7.23P realizing a similar result with a constant w e g has been proposed, although 
in the present cas e oi H in Eq. ( I7.19P we have a dynamical w c g. In addition, it has been 



examined in Ref. 



113| that in more general case, the ratio of the energy density of dark 



energy to that of matter is not constant. 

In addition, by extending the relation (I7.14p . we examine more general cases that Lh 
depends on the Hubble parameter H and the length scale which originates from the cosmo- 
logical constant as A = 12// 2 , i.e., L h is expressed as 



Lh — L h (L ph , 1/^, t s , H, I) 



(7.27) 



or 



Lh — L h (L ph , Lfh, t s , if, I) . (7.28) 
The proposal of generalized holographic dark energy in the form in Eqs. (I7.27P or (I7.28P has 



been made in Ref. 



271 ] where instead of H, the scale factor a was used, supposing that such 



a cut-off may depend on the scale factor and its derivatives (i.e., also from H). Hence, this 
is the most general proposal for the IR cut-off which eventually covers all known proposals. 
An example in such an extended class of generalized holographic dark energy is given by 

2" 



Ch 
Lh 



1 



ahL 



rii 



1 + 2 (al 



2«h 



(7.29) 



where c*h(> 0) is a positive dimensionless parameter. Here, we find L h > due to 
Ch > and «h(> 0). This can be seen from the following relations: (a^ — ah — l) 2 — 



2 (a* 



2al 



a h + 1) 
«h («h — l) 2 + 1 > 0. 



[at 



1/2Y 



2a h - 3/4 < and al - 2a\ + a h + 1 

£« h +i e t/i 



L c a h [1 + tj (a h l)] 
t 



j ih 



ah [1 + tj (a h l)] 

with L c being an integration constant. The solution in ( IT.30I) with (IT.3 II) yields 



(7.30) 
(7.31) 



H 



1 + a h [1 + t/ {a h l)Y 
t[l+t/ (a h l)] 
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(7.32) 



We mention the case that «h is negative. From the denominator of Eq. (17.30ft . we see 
that at t — —a^l a Big Rip singularity appears because a diverges. Furthermore, can 
be negative, so that when matter is not included, H = C^/ L h can also be negative and 
therefore the universe will be shrink. 

When t is small, i.e., in the limit of t — > 0, from Eq. (I7.32p we find 

H -+ + i . (7.33) 

This means that in this limit, the energy of the universe is dominated by that of a fluid 
whose Eos is given by w m = — («h — 1) / [3 («h + !)]• While, when t is large, i.e., in the 
opposite limit of t — > oo, from Eq. (I7.32p we see that H becomes a constant as 

H ->■ - . (7.34) 

This implies that the universe asymptotically approaches to de Sitter space. We explore the 
following model: 

2 ~ ' wr x r + (My (7 . 35) 

Also in this case, the solution ( 17.301) with ( 17. 311) is again satisfied. 



0, = H 2 (ag - a h - 1) / Lh\ 

L h a h + l L m I a h +l \a h l J 



C. The Hubble entropy in the holographic principle 

It follows from Eqs. (I3.3p and (17. ip that the Friedmann equation is described by 3H 2 /n 2 = 
3C*h/ (K 2 L^)+p m . We define the energy of matter E m , the Casimir energy Eq and the Hubble 
entropy Sh as 

E m = Pm L 3 h , (7.36) 
Ec = ^ • (7.37) 

^h^^^. (7.38) 



By substituting Eqs. ( I7.36p -( 17.38p into the Friedmann equation sho wn ab ove, we have 114 1 



the Cardy-Verlinde holographic formula of the Friedmann equation H 



SI = {k 2 L{) 2 E C (E C + E m ) . (7.39) 



75 



Suppose that p m can be neglected, Eq. (17.381) is rewritten to 

18nC 2 L 2 18t< 

Su = = -^w ' (7 - 40) 

where in deriving the second and third equalities, we have eliminated H and by using 
Eq. (17. 3p . respectively. It is seen from Eq. (17.401) that Sh depends on time, not a constant 



as in case of de Sitter universe in Ref. 
(I7.19p . from Eq. (I7.40p we acquire 



1151 ]. In the model in Eq. (I7.15P with the solution 



18ttC^ 2 (t s - 1) 



2 



*hie ■ {7A1) 



c s 



Therefore, in this model Sh = at t = and t = t s , and Sh becomes maximum at t — t s /2. 

On the other hand, for a generalized model in Eq. (I7.24p with the interaction represented 
in Eq. (I7.23p . we have 

ft = A. ^V'W («.-«)', (7.42) 



3/i c J kP 1 h^t^ 

In comparison with Eq. (17.4 ip . the form in Eq. (I7.42p is multiplied by the first constant 
factor on the right-hand side. In case of Eq. (I7.29P or the Friedmann equation 3H 2 /k 2 = 
3Cl/ {i^L 2 ) + p m , we have 

Sh= i8,qy|i + «/M)] 2 (743) 

K 2 {l + a h [l+t/ (a h l)] 2 } 
Therefore, in this model Sh = at t = 0, whereas, in the limit of t — > oo, Sh approaches to 
a constant as follows. 

S^ 1 -^. (7.44) 

We remark that Sh in Eq. (I7.43P is positive, even though a c < 0. Moreover, Sh in 
Eq. (17.401) is always positive. However, for the case that Sh is given by Eq. (17.421) . if 
K 2 p c / (3h c ) > 1, Sh can be negative. This implies that entropy of the universe should be 
nega tive, provided that Sh corresponds to the upper bound on the entropy of the universe. In 



Ref. 



1161 ]. negative entropy in the phantom phase has been observed. While, if the phantom 



phase is transient in the late time, the entropy of the universe may remain positive [29]. 

To connect the holographic dark energy scenario with the reconstruction of the corre- 
sponding scalar field theory in Sec. IV, we explore another model. 

^ = fr 3 -^r + 7, (7.45) 

76 



where a(> 0), 0(> 0) and 7(> 0) are positive constants and satisfy 

2B[Jj 



1 > 



3 V a 



(7.46) 



Here, T is defined by 



T 



H - 7 + 4a/ (3fi 2 ) (H + 3f) In (5/8) 



a/ (3/3 2 ) (H + 3 7 ) 2 - f3 - 4a/ (3/3 2 ) In (2/6) 

P _ «(*) p, 
"~o(0) ' 

»+oo 

O = / dtexp 



^t^ + h 2 -^ 

12 2 ' 



(7.47) 
(7.48) 
(7.49) 



Equation (17.471) informs us that if becomes zero only once as a function of T. We also have 

OL. 



7 



■fT c 3 + ^T c . 



(7.50) 



Here, T c (< 0) is a negative constant and it is determined that when T = T c , H(T = T c ) = 0. 

We reconstruct a corresponding scalar field theory in the holographic dark energy sce- 
nario. We take a concrete form of 7(0) in Eq. ( I4.16P as 



a 



1(0) = - (T c + 



^(T c + 0)+ 7 



(7.51) 



where 7 is given by Eq. ( I7.50j) . In this model, uj(<f)) and V(<f>) in the action in Eq. ( 14. ID is 
expressed as 



(0) = -_[a(T c + 0) 2 -^] 



V{4>) 



ft" 



(T c 



- 2a$ (T c 



+ « 7 (T c 



+ (ft + 3^ 2 ) (T c 



2^7(T c + 0) + 3 7 2 ] . 



(7.52) 



(7.53) 



It follows from Eq. (14.161) that we have the solution = t, H — I(t), and hence a is given 



by 



a = a c exp 



^(Tc + t) 4 -|(Tc + t) 2 + 7 (T c + t) 



(7.54) 



The scale factor a has a minimum because H = at t = 0. Therefore, for t < the universe 
will shrink, whereas for t > it will expand. From H = a (T c + t) 2 — j3, we see that H — 
at 



t = t± = -T c ± \l '- > . 

ft 



(7.55) 
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By using Eqs. (13. 9 j) and iJ = i(t) with Eq. (I4.16P and = t, w e g is written as 

2 \a (T c + t) 2 - ^1 , . 

iweff = -1 L , - • 7.56 

a(T c + 0) 3 -3/3(T c + 0) + 37 

Accordingly, it is seen from Eq. (13. 9 p that for t_ < t <t + , w c s > — 1, i.e., the universe is in 
the non-phantom phase, whereas for < t < £_ or t + < t, w c g < — 1, i.e., the universe is in 
the phantom phase, and that at t — 1±, w e g — — 1. In summary, at i = the crossing 

of the phantom divide can occur from the non-phantom (phantom) phase to the phantom 
(non-phantom) one. 

In the model in Eq. (I7.45p . the Hubble entropy in Eq. (I7.38P is given by 

s «*SL ( , 57) 

K 2 [(a/3) (( + T c ) 3 -/3(( + T c )+7] 

Thus, Sh is always positive. It follows from Eqs. (I7.50P and (I7.57P that at t — 0, S'h diverges. 
From Eq. (I7.40p . we see that S'h oc H~ 2 . Thus, for < t < t_ (the phantom phase with 
if > 0), S'h decreases. For t- < t < t + (the non-phantom phase with H < 0), Sh increases. 
For t + < t (the phantom phase with if > 0), Sh again becomes small. Eventually, in the 
limit of t — > oo, Sh — > 0. As a result, in the non-phantom phase Sh grows, whereas in the 
phantom phase, Sh decreases. 



VIII. ACCELERATING COSMOLOGY IN F{R) GRAVITY 



In this section, we study an accelerating cosmology in F(R) gravity. First, we consider re- 
lations between a scalar field theory in the Einstein frame and an F(R) theory in the Jordan 
frame. Furthermore, we show how to obtain the ACDM, phantom-like or quintessence-like 
cosmologies in F(R) gravity by following Refs. 117H120| . We mention that the reconstru- 
crion of F(R, T st ) gr avity has also been investigated, where T st is the trace of the stress-energy 
tensor, e.g., in Ref. [121| . 



A. F(R) gravity and a corresponding scalar field theory 



The action describing F(R) gravity with matter is given by 

-F(R) 



S = d Xy/—g- 



2k? 



d xC M (flV, * 



MJ 



(8.1) 
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where Cm is the Lagrangian of matter. By making a conformal transformation, we move to 
the Einstein frame: 



— O 2 

w _ dF{R) 
dR 



Moreover, R is represented by using R as 



R = e 



1/V3k4 



1 



R + VSU (k6) - -g^du (k6) d v (k 



with 



0(kc 



where the potential V(4>) is represented by 



v{<j>) 



f r r-f 



\2 • 



(8.2) 
(8.3) 
(8.4) 



where a tilde denotes quantities in the Einstein frame. We define a new scalar field as 



5.5) 



(8.6) 



(8.7) 



As a result, we acquire the action in the Einstein frame 122 ] 
S E = J d 4 x^g I £5 - hi^d^d^ - V(<j>) J + J d 4 xC M ((Fr)" 1 (0)^, *m) , (8.8) 



■9) 



2k 2 (F r , 

We note that an important cosmological application of the relations between F(R) gravity 
in the Jordan frame and its corresponding scalar field theory in the Einstein fra me to the 
time variation of the fine struct ure c onstant in non-minimal Maxwell- _F(i?) gravity 123| has 
recently been executed in Ref. 124| by using a novel consequence of a static domain wall 
solution 125| in F(R) gravity. Such non-minimal Maxwell theories with its coupling to a 
scalar field or the scalar curvature break the conformal invariance of the elect rom agnetic 
fields, so that the large-scale magnetic fields from inflation can be generated 1261 ] . This 
would be considered to be a significant cosmological implication of the investigations of the 
present section. 
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B. Reconstruction method of F(R) gravity 



To begin with, we explain the reconstruction method of F(R) gravity [1171 Ill8| . We 
introduce proper functions P{4>) and Q(<fr) of a scalar field and rewrite the action in 
Eq. (18.1 p to the following form 

S = J Avq^ (P(<f))R + Q{<f>)) + J d^xC M (g,*,, *m) • (8.10) 

Since the scalar field does not have the kinetic term, it may be regarded as an auxiliary 
scalar field. It follows from Eq. (I8.10P that the equation of motion of reads 

= i_rm R+ iQM (8U) 

d<p dcp 

Hence, by solving Eq. ( 18. lip in terms of R in principle we find the expression = <f)(R). 
By combining this expression and the action in Eq. (I8.10p . we obtain the representation of 
F(R) as 

F(R) = P((f>(R))R + Q((f>(R)). (8.12) 
Furthermore, the variation of the action in Eq. (I8.10p leads to the gravitational field equation 



- 9(a/ (P(0)i? + Q(0)) - R^Pty) - v,DP(0) + V M V,P(0) + k 2 T^ = . (8.13) 



Here, V M is the covariant derivative operator associated with and □ = p^V^Vj, is the 
covariant d'Alembertian for a scalar field. Moreover, T$P is the energy-momentum tensor 
of matter. We take the flat FLRW space-time ds 2 = —dt 2 + a 2 (t) X^=i 23 {dx 1 ) 2 . In this 
background, the (/i, v) = (0, 0) and (/i, v) = — 1, • • • , 3) components of Eq. (I8.13P 

become 



-6H 2 P(<P(t)) - Q(0(t)) - QH dP ^ t {t)) + 2k 2 P m = , (8.14) 
2 f£mi + 4H^t)l + (± H + 6H ^ P{m + Q{m + 2 K 2 P U = . (8.15) 

Here, we have expressed the sum of the energy density and pressure of matters with a con- 
stant EoS WMi as pm and Pm, respectively, where the subscription "i" denotes a component 
of matters. By using Eqs. (18. 14p and ( 18. 15ft . we eliminate Q{4>), and eventually we have 

— iJ — \-2HP(<j)(t)) + k (pm + Pm) = 0. (8.16) 



SO 



If we appropriately redefine the scalar field 0, it can be taken as = t. In addition, we 
describe the form of a(t) by 

a(t) = a exp (<?(£)) , (8.17) 

with a being a constant and g(t) being a proper function of t. In this case, the Hubble 
parameter is given by H = dg(<f)) / (d<fr). By using this expression, Eq. (18.161) can be rewritten 
to 

d 2 P(<t>) _ a^dP{^)_ + 2 jW) p(fl 
d(p 2 d<f) d<p dip 2 

+ k 2 (! + wmO p Mi a" 3(1+u,M!) exp [-3 (1 + w Mi ) ^(0)] = . (8.18) 

i 

Here, pMi is a constant. In addition, by solving Eq. (18.141) in terms of Q{4>), we acquire 

~dg{4>)~ 



Q(0) = -6 



d~g{<j>)dP{(j)) 
d^ d* 

2K 2 J2p l d~ 3{1+WMl) exp [-3 (1 + u> M *) ~g{4>)\ ■ (8.19) 



5 = / d 4 x^4/ + / d 4 x£ M (g,u, *m) , (8.20) 



It is significant to emphasize that by redefining the auxiliary scalar field <j) as <fi = $(y?) 
with a proper function $ and defining P(<^) = P($(</?)) and Q(<p) = Q($(v?)), we obtain 
the new form of the action 

:HR) 
2k 2 

where 

F{R) = P{ip)R+Q{cp). (8.21) 

Since F(R) = F(R), this action in Eq. (I8.2ip is equivalent to that in Eq. (I8.10p . Furthermore, 
tp is the inverse function of $, and therefore by using <fi = <j)(R) (p can be solved with respect 
to R as <f = f(R) = < &~ 1 (0(-R)). Accordingly, there exist the choices in as a gauge 
symmetry, and hence can be identified with time t as = t. This can be c onsid ered as 



a gauge condition which corresponds to the reparameterization of = 4>(ip) 119| . As a 
result, if we obtain the solution t = t(R), by solving Eq. (I8.18P and (I8.19P and substituting 
these solutions into Eq. (18.121) . the explicit expression of F(R) can be acquired. It should 
be noted that in a naive model of F(R) gravity the crossing of the phantom divide cannot 
be realized because F(R) has to be a double- valued function in order that the crossing of 
the phantom divide can occur. In fact, however, if the action of F(R) gravity is extended 
to the form of P(0)i? + Q(0), the crossing of the phantom divide can happen. We show an 
explicit example to realize the crossing of the phantom divide in Sec. VIII C 3. 
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C. Reconstructed F(R) forms and its cosmologies 



1. The ACDM cosmology 



We demonst rate the reconstruction process of F(R) gravity in which the ACDM cosmol- 



ogy is realized 



1171 ] (f or ot her method of the reconstruction of F(R) gravity reproducing 



the ACDM model, see 



1271 ]) . In the flat FLRW background, if there exist a matter with its 



EoS wm and the cosmological constant A, the Friedmann equation (I3.3P is written as 

- 2 1 



H = — Pmc« 



-3(1+u> m ) 



72 ' 
'A 



(8.22) 



where pMc is a constant and I a is a length scale related to the cosmological constant A. From 
Eq. (18. 22 p . we find that the scale factor can be expressed by Eq. (18.171) with g(t) being 



m 



3 (1 + w M ) 



In 



K 2 l\ 



( 3(1 + WM) (t~t A ) 



(8.23) 



where t& corresponds to an integration constant. We suppose that geometrical dark energy 
originating from F(R) gravity is dominant over matter and therefore matter contribution 
can be neglected. In this case, by using Eq. (I8.23p . Eq. (I8.18P is expressed as 



d 2 P((f)) 
dcf) 2 



lcoth( 3 ' 1 + mM )( f - tA )^ PW 
l A \ 2/ A 

3 (1 + WM) sinh" 2 ( 3(1 _ + " M) (t - t A ) ) P(0) = . 



72 
'A 



2l A 



(8.24) 



This expression can further been rewritten to a Gauss's hypergeometric differential equation 
by replacing the variable <j) with z as 



z(l-z 



d 2 P(z) 



7 



sinh 



dz 2 

I 3 (1 + w M ) 

21a 



a + P + l) z 
(t-t A ) ) , 



dP{z) 
dz 



- aj3P(z) = 



(8.25) 
(8.26) 



with 



7 = 4 + 
a + /3 + l = 6 + 

a(3 = 



1 



3 (1 + w M ) 
1 

3(1 + w M ) 
1 



3 (1 + w M ) ' 



(8.27) 
(8.28) 
(8.29) 
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By using the Gauss's hypergeometric function, a solution of Eq. (I8.25P is described as 

P(z) = P cl P G (a, 0, f, z) + P c2 (1 - zf- & -$ F G (7 - 5, 7 - 7; 2) , (8.30) 
where P c \ and P C 2 are constants and Fq is the Gauss's hypergeometric function, defined by 

r(a)ros)^ r(7 + n) n\ 

with r being the V function. For simplicity in order to obtain the form of F(R), we set 
P C 2 = 0, namely, we take only the first linearly independent solution in Eq. ( I8.30p . It follows 
from Eq. ( I8.19P that we obtain the form of Q as 



= - 3(1 z)Pcl ( 2F G {a, 0, 7; z) + 1 ( - + ^ + - + ^ j _ (g 32) 



In the limit of t = — >■ oo, from Eq. ( 18.261) we see z — > 0. Thus, in this limit, by substituting 
Eqs. (IQjjl and fl8T32|) into Eq. fl8TT2p we find 



F(R) = P(<f>(R))R + Q{<P{R)) -)• Pel f - p J ■ (8.33) 
As a consequence, by comparing Eq. f!8.33j) with Eq. (12. ip . we see that if we take 

Per = ^ , (8.34) 



^ = 2A, (8.35) 

'A 

the general relativity with the cosmological constant is realized. Moreover, by plugging 
Eq. (I8.33P with Eqs. (I8.34p and (I8.35P into the action in Eq. (I8.8p . we acquire the action 
describing the ACDM cosmology which corresponds to the one in Eq. (I2.ip . 

2. Quintessence cosmology 



Next, we reconstruct F(R) gravity in which quintessence-like cosmology is produced 117 ]. 
We investigate the case that g{t) in Eq. (I8.17P is given by 

£(</>) = %) In f£) , (8.36) 

where <p c is a constant and h(4>) is a function varying slowly in 0. Therefore, an adiabatic 
approximation is applied to h(4>), so that the derivative of h{4>) can be neglected, i.e., 
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dh{4>) j d(f) ~ d 2 h(4>)/d(j) 2 ~ 0. Through the same procedure as in Sec. VIII C 1, we derive 
the solutions of P{4>) and Q{4>)- The equation of P(4>) in (18.181) is given by 

d 2 P(<P) _ ~m dP(4) _ 2 H^ P U) 
d<fi 2 (ft d(p (j) 2 

/ i \ -3{l+w M i)h{<f>) 

+ k 2 (1 + ™Mi) PM t a- 3{1+WM > ) (f) = . (8.37) 

We acquire a solution of Eq. ( I8.37P as 

P{<P) = P + 4> n+ + P_0 n - + P i (0)0- 3 ( 1+w Mi)^W+2 ? ( 8>38 ) 

i 

n±(0) = ~ - 1 ± y/h^ + Gh^ + lj , (8.39) 

p/,\ = (i + w M »)PM»q | M,; 0c = (saq) 

A) ~ 6(l + w M i)(4 + 3 M ;MO/i 2 (0)-2(13 + 9M;M J )M0)+4' { ' ' 

where P± are arbitrary constants, and Wmt — 1/3 for radiation and % m = for non- 
relativistic matter. Furthermore, by using Eq. ( 18. 191) . Q is written by 



Q{<t>) = -6h(<f>) (P + CH<t>) + n + (0))0 n+w " 2 + P-{h{<f>) + n_(<P))<p n ~^- 2 ) 

+ J2 {-^{<f>) [~ (2 + 3w M i) K<t>) + 2] Pi(<f>) + p Ml d^ 1+w ^^ l+w ^ h ^] 

i 

x 0~ 3 ( 1+u, Mi)M<M _ (8.41) 



On the other hand, from Eq. (18.361) we find H ~ h(t)/t and hence R = 6 yH + 2H 2 
6 (— + 2h 2 {f)) /t 2 . Here, we provided that in the limit of — >■ 0, h{4>) — > hi, and in the 
opposite limit of <fi — > oo, h{4>) —> hf. As a form of h(4>), we take 



IT 



*(0 = ^r, ( 8 - 42 ) 



where $ is a small constant enough for h(<f)) to be a function varying slowly in <fi. The 
substitution of Eqs. M and (ESQ into Eq. (13321) yields 

F(P) = P($ c (i?)) J R + Q($ c (P)), (8.43) 
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where 

$ C (R) = yj&/ a + &l a , (8.44) 
$ ± = J f - A ± \/ /5 " §H , (8-45) 




01 = 



16 [J? + 3/if (l — 2/if) 0]' 

6/ii (1 - 2/i f ) , (8.46) 



270 3 /? 3 

2 [R + 3hf (l - 2/i f ) 0] [i? + 2 (3/n + 3/i f - 2/ii/if) 0] 



30i? 

4 [i? + 3/if (l - 2/i f ) 0] 2 # + 2(3/^ + 3^-2/^)0 
A = 30^ 0^R ■ (8 - 47) 

At the dark energy dominated stage, w DE ps w eff because the energy density of matter can 
be negligible compared with that of dark energy. In the present model in Eq. (I8.42p . from 
Eq. (I3.9P we have wrw. = —1+2/ (3/if). Thus, for < h { < 1, — 1 < wde < —1/3. As a result, 
this means that in the reconstructed F(R) gravity model, quintessence-like cosmology can be 
realized. In addition, for the reconstructed F(R) gravity model in Eq. (I8.43p . in the late limit 
of 4> — > oo, which can be regarded as the limit of the present time, the asymptotic behavior 
is given by a power-law description as F(R) ~ R s with s = — — 5 + \/h^ + 6/if + 1 \ /4. 
By using this expression, we find that in the action of a scalar field theory in Eq. (18. 8ft . the 
potential V(<j>) in Eq. (18. 9 p is written as 



1 5_i i 5_i ri i -(?-2)/(--i) 



I e (2/3)K0 



(8.48) 



2k, 2 s 2 Ik 2 s 2 

Here, in deriving the second equality, we have used R = [(1/s) e*- 2 / 3 -'^] 1 s 1 \ which follows 
from Eq. (18.51) . This can be interpreted as quintessence potential. 



3. Phantom cosmology 



Further more , we reconstruct F(R) gravity in which the crossing of the phantom divide 
is realized 120] and eventually phantom-like cosmology is produced. 

For matter to be neglected because of the dark energy domination, as an example, we 
examine the case that g(t) in Eq. (18.171) is described by 



-10 In 



-a 



X7+1 
1 to) 



(8.49) 
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where 7(> 0) and C p (> 0) are positive constants and to is the present time. We note that 
since there occurs a Big Rip singularity at <p = t s = t C p 1// ^ 2 ^ +1 ^ W e investigate the period 
< t < t s in order for g{4>) to be a real number. From Eq. (I8.49p . the expression of the 
Hubble parameter H(t) = dg(<ft)/d<j) is given by 

7(0Ao)- ( " +1) + (7 + l)C p (0A o r 



H(t) = 



\t 



7+1 



10 



tsCp 



1/(27+1) 



(0/to)" 7 -C p (0/t o 
7 + (7 + l)(tAs) 27+1 



(t/t 6 



\27+l 



where in deriving Eq. (I8.5ip we have used the relation t 
solution of equation ( 18 . 1 8[) for P(<f>) is derived as 

2 



-1/(27+1) 



and 



P{4>) = exp 

1 ± y/l + 1007(7 + f 



0: 



± 



(8.50) 

(8.51) 
: t. A 

(8.52) 
(8.53) 



with p± being arbitrary constants. 

It follows from Eq. (13. 9p that u>de ~ ^ e a a t the dark energy dominated stage is expressed 

as 

, -7 + 47 (7 + 1) (t/t s ) 27+1 + (7 + 1) (t/t s ) 2(2 ^ +1) 

"^DE — — J- ; ; • (O.04J 



15[7+(7+l)(t/t s )^ +1 ]' 

In the limit of t -> 0, namely, i/i B < 1, from Eqs. (l830|) and ( 1H34|) we find if (t) ~ 107/t and 

wbe ~ — 1 + 1/ (157) (> —1), respectively. This is the non-phantom (quintessence) phase. 

While, in the opposite limit of t — > t B it follows from Eqs. ( I8.50P and (I8.54p that H(t) ~ 

10/ (t a — t), which leads to a(t) ~ a(t s — t)~ 10 , and wbe ~ — 16/15(< —1), respectively. 

Moreover, d(wDE + 1) /dt monotonously decreases in time. Hence, first the universe is in 

the non-phantom phase. As the time passes, when t closes to t s , the universe enters the 

phantom phase. Thus, the crossing of the phantom divide line of u>de = — 1 can occur at 

r ( 1 1/(27+1) 

t = t c = t s [-27 + v/47 2 + 7/ (7 + 1) 

In addition, by using Eq. (I8.52p and substituting it into Eq. ( I8.19p . we acquire the forms 

of P{t) and Q{t) as 

(t/t y 



p( t ) 

Q(t) 



-1 5 



1 - 

-6H 



(t/t s r +i 

' (t/tof 

l-{t/t s ) w 



(8.55) 
(8.56) 
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On the other hand, from R = 6 + 2H 2 j , the relation between R and t is given by 
7 (20 7 - 1) + 447 (7 + 1) (t/t s ) w + (7 + 1) (207 + 21) (t/i s ) 2(2 ^ +1) 1 



60 



i2 



(8.57) 



*» [i-(t/t s r + T 

Accordingly, in principle, if we obtain the relation t = by solving Eq. ( I8.57P reversely, 
the substitution of it into Eqs. (18.551) and (I8.56P and the combination of those with Eq. (I8.12p 
yield the explicit form of F(R). We show analytically solvable cases below. In the limit 



of t -»• (t/t s < 1), from Eq. (15371) we obtain t ~ ^607 (2O7 - 1) / R. By using this 



asymptotic relation, we acquire 



128 



F(R) 



^x/607 (207- 1) 



1 57 



to 



R -5l/2+l 



57 - & - 1 

207-1 



[6O7 (207 - 1)]^ /2 i?-^ /2 



f8.58) 



In the opposite limit of t t s , it follows from Eq. f[g37D that t~t s - 3y/14Q/R. With this 
asymptotic relation, for the large curvature regime as t^R > 1 we find 



F(R) w Fi? 7/2 




(8.59) 
(8.60) 



.3\/l40(27 + 1) Vt . 

For the power-law form of F{R) in Eq. (18.601) . the potential V{4>) in Eq. (18. 9p of the action 
for a scalar field theory in Eq. ( 18. 8ft is described as 



2 (2/3)k0 



7 



-3/5 



(8.61) 
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This can be interpreted as phantom potential. 

Using the reconstruction program with auxiliary scalar fi elds as discussed in this section, 
or without the use of auxiliary scalar fields, following Ref. 1291 ] . one can eventually recon- 
struct any dark energy cosmology studied in this review. For instance, Little rip cosmology 



for modified gravity has been presented in Refs. 



71 



73(]. 



D. Dark energy cosmology in F(R) Hofava-Lifshitz gravity 

As a candidate for a renormalizable gravi tatio nal theory in four dimensions, the Hofava- 
Lifshitz gravity h as be en proposed in Ref. 



mology, see, e.g., 



1301 ] (for a review on the Hofava-Lifshitz cos- 



131] ), although it cannot maintain the Lorentz invariance. In addition, 
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its extension to an F(R) formalism has been executed in Ref. 132 1 . In this subsection, we 
study cosmology for dark energy in F(R) Hofava-Lifshitz gravity 133 1 . 



1. F(R) Hofava-Lifshitz gravity 



by 



he model action describing F(R) Hofava-Lifshitz gravity with the matter action is given 



132) 



S = J dtd 3 xy/g^^N l F(R) + J d 4 xC M , 

KijK ij - XklK 2 + R^ + 2/i HL V M {n< x V v n v - n"V ' v n») 



R 

K = gVKy = g ij ^ (//;;' - vl 3) Nj - vf ] N t 



(8.62) 

(8.63) 
(8.64) 



with 



kl 



G 



ijkl 
AhL = 

E ij = 



9ik 9jl 



E^GijkiE 
1 
2 

AhL 

3A HL - 1 
1 5W[g 



9u gjk 



AHLfl'ij gn > 



(3)l 
kl J 



(8.65) 
(8.66) 

(8.67) 
(8.68) 



'9" 

where, k, I run over 1, 2, 3, N\ is the lapse variable in the Arnowit t-Deser-Misner (ADM) 



decomposition in (3 + 1) space-time, N l is the shift 3- vector [134J, [135], Cm is the Lagrangian 
of matter which is a perfect fluid, K~ij is the extrinsic curvature, R^ is the spatial scalar 
curvature, n M is a unit vector perpendicular to a constant time hypersurface, R is the di- 
mensionless gravitational coupling, and Ahl and /^hl are constants and these break the full 
diffeomorphism invariance. Furthermore, G^m is the inverse of the generalized De Witt 
metric and it exists only for Ahl ^1/3 because it follows from Eq. ( 18.671) that if Ahl = 1/3, 
Gijki becomes singular. Moreover, E % i is constructed so that t he d etailed balance principle 
restricting the number of free model parameters can be met 1301 ] . In the Hofava-Lifshitz 
gravity, there exists the difference of the scaling properties between the space and time co- 
ordinates as x % — > bx l and t — > b z t, with b being a constant and z being a dynamical critical 
exponent. Here, if z — 3 in (3 + 1) space-time dimensions, the theory is renormalizable, 
whereas for z — 1, it is the general relativity. Such scaling properties give the theory only 



SS 



the foliation preserving diffeomorphisms, given by 5x l = ((x l ,t) and 5t = £(£), where ( and 
£ are functions of x l as well as t and t, respectively. We also mention that in the Hofava- 
Lifshitz gravity 1301 ] . Ni is supposed to depend only on time in order for the projectability 
condition to be satisfied and it is set to be unity (N\ = 1) by applying the foliation preserving 
diffeom orph isms. For z = 2 and z = 3, the expression of W^ly ] in Eq. (I8.68P are presented 



in Ref. 



1361 ] . In the spatially flat FLRW space-time ds 2 
3(1-3Ahl + 6/i H l)# 2 , 6/i 



-N?dt 2 + j: t=l „Jd 



x 



R 



11L 



d ( H 



N 2 



By varying the action in Eq. (I8.62p with respect to N\ and g\, , we acquire 



Ni dt \Ni 

(3) 



we have 
(8.69) 



drx 



— ll , — . 



- 6 (1 - 3A HL + 3// H l) ^ 2 - 6//hl (H-HRf"(R) )-k 2 Pm = 0, (8.70) 



(8.71) 



F(R) - 2 (1 - 3A HL + 3/i HL ) I H + 3H ) F (R) — 2(1— 3A HL ) RF (R) 



+ 2/ihl ( R 2 F"'(R) + RF"{R) ) + k 2 Pj 



M 







where the prime denotes the derivative with respect to .R, and pu and Pm are the energy 
density and pressure of a perfect fluid, respectively. In deriving Eq. ( 18. 70 p . we have used the 
projectability condition for JVj to have only the time dependence, and in obtaining Eq. f)8.7ip . 
we have taken N\ = 1. For Ahl = Hkl = 1, these resultant equations are reduced to those 
in ordinary F(R) gravity. 

By using the continuity equation pu + 3H (pm + -Pm) = in terms of a perfect fluid and 
executing the integration of Eq. f)8.70p . we find 



F(R)-6 (1 - 3A HL + 3/i HL ) H 2 + fi HL H F (R)+Q/i HL HRF (R)-k 2 Pm 



liL 



0, (8.72) 



with C*hl being a constant of integration, which has to be chosen to so that the constraint 
equation (18.701) can be satisfied. 

We note that for F(R) = R, in the flat FLRW background the gravitational field equations 
are written as H 2 = {k 2 / [3 (3A H l — 1)]} Pm and H = — {k 2 / [2 (3A H l — 1)]} (pm + -Pm) with 
Ahl > 1/3 due to the consistency, and for Ahl — > 1 these equations become the ordinary 
Einstein equations in general relativity. 



2. Reconstruction of F(R) form 

We further analyze Eq. (18. 72 p . It follows from the continuity equation of a perfect fluid 
with its constant EoS w M = Pm/ Pm that p M = pM C a" 3(1+WM) e" 3(1+u,M)Ar with p Mc being 
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a constant, where N = In (a/a c ) with a c being a constant is the number of e-folds. We 
replace the cosmic time t as a variable with N, so that Eq. (I8.69j) can be rewritten to 
R = 3(1- 3A HL + 6/xhl) G(N) + 3/i HL (dG(N)/dN), where we have defined G(JV) = H 2 
and used it in order to analyze Eq. (I8.69P easier. Since this equation can be solved as 
N = N(R) and H can be represented as H = H(N), Eq. (I8.69P can be described as an 
equation of F(R) in terms of R. Accordingly, Eq. ( I8.69P can be rewritten to 



F(R) - 6 



18/iHL 



1 - 3A HL + 3/xhl) G + 



/j HL dG(N) 
2 dN 



dF(R) 
dR 



,^dG „d 2 G{N) 
3A H l + 6 Mhl )G— + llm £- m - 



d 2 F(R) 
dR 2 



K 2 p Mc a^ l+w ^e-^ 1+w ^ N 



0. 



(8.73) 



First, we reconstruct the form of F(R) with realizing the ACDM cosmology, in which the 
Friedmann equation ( 13. 3 j) can be described by G(N) = H 2 = H 2 + (k 2 /3) Pm c <\ 3e_3Ar with 
H c being a constant, where we have used wm = because a perfect fluid is considered to a 
non-relativistic matter. For general relativity, H 2 = A/3 as seen in Eq. (]2.4p . Moreover, we 
have e" 3Ar = [R-3{1- 3A HL + 6/i HL ) #c] / {^PM C a c ^ 3 [1 + 3 (/i HL - A HL )]}- With these 
equations and % = 0, Eq. (I8.73P is represented as 

(1 - 3A HL + 3/xhl) F(R) -2(l- 3A HL + ^hl) R + 9/i H l (1 - 3A HL ) H l^J^ 



6/lHL (-R - 9/iHL^c) - 3H 2 (1 - 3A HL + 6 /iHL 

J R-3(l-3A H L + 6/i H L)^ c 2 = 0, 



rf 2 F(£) 
dR 2 



(8.74) 



The homogeneous part of Eq. (18.741) is rewritten to 

T (l- T )f£^L + [ 7 _ (a + ^+1) r] - c^F(£) = 



r = 



dT 2 
R - 9/i HL # ( 



dr 



3H 2 [1 + 3 (/i HL - A HL )] 



(8.75) 
(8.76) 



where 



1 

a + (3 + 1 

a/3 



1 

"2 ' 

1 - 3 (A HL + /xhl/2) 

3/iHL 
1+3 (/i H L - A H l) 



6/i 



HL 



(8.77) 
(8.78) 

(8.79) 
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We describe the complete solution of Eq. (I8.76P with the Gauss's hypergeometric function 
Fq as 

F{R) = F cl F G (a, 0, 7; r) + F c2 r 1 ^F G (a - 7 + 1, $ - 7 + 1, -7 + 2; r) 

R-2A, (8.80) 



3A HL - 1 



3 [1 - 3 (A HL - 3^ HL )] g c 2 
A "" 2[1-3(A hl -^hl)] ' (8 ' 81) 

with F cl and F c2 being constants. Thus, a class of the reconstructed F(R) theories in 

Eq. (I8.8ip can represent the ACDM cosmology. We remark that for /ihl = Ahl — 1/3, 

R becomes a constant, so that the solution in Eq. f l8.81D can be expressed as F(R) = 

F(R)/ (3A HL - 1) - 2A, where A = (3/2) (3A HL - 1) # 2 - 

Next, we reconstruct an F(R) form describing the phantom cosmology. Since we examine 

the dark energy dominated stage, for simplicity, non-relativistic matter contributions are 

neglected. At the dark energy dominated stage, by using the continuity equation in terms 

of dark energy, the Hubble parameter can be represented as H = H p h/ (t s — t) with = 

— 1/3 (1 + wjde), where at t — t s , a Big Rip singularity appears. In what follows, From this 

expression, we have G(N) = H 2 (N) = iJ p he 2Ar ^ ph . In this case, the relation of N to R is 

described by e 2N ^ Hph = Rj [H ph (A ph H ph + 6/xhl)] with A ph being a constant. By combining 

these relations and Eq. (I8.73p . we acquire the Euler equation 

,,, _ {A ph H ph + 6/^hl) (A ph H ph + 3/i H l) r Q Q q\ 

6yU HL {A ph H ph + 12/xhl) 

W - (A>h# P h + 6^hl) 2 , s 

Wph2 = 12^^ + 12^) " ( } 

A solution of Eq. (I8.82p is given by 

F(R) = F ph+ RF*> h+ + F ph _R^- , (8.85) 



W phl ± V /(W ph i-l) 2 -4W ph2 



V P h± = — - , (8.86) 

where -F p h+ and -F p h_ are constants. Hence, a reconstructed form of F(R) in Eq. (18. 85ft can 
describe the phantom cosmology. As a result, it is considered that by using the procedure 
explained above, in principle, an F(R) form with representing any cosmology could be 
reconstructed. 
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E. F(R,T) gravity 



In Ref. 



34j . the formulations of a novel modified gravitational theory, the so-called 



F(R, T) gravity with T being the trace of the stress-energy tensor, which can explain the 
late-time cosmic acceleration, have been investigated. In this subsection, we review this 
latest theory. 



1. Formulations 



The action of F(R, T) gravity is given by 34] 

-F(R,T) 



S 



d Xy 



16tt 



d X\J—gL 



M 



(8.87) 



where T = g^T^ is the trace of the stress-energy tensor of matter, defined as [137] T^ v = 
— (2/ \J—g) 8 (^—gLu) /SgQ^u, Lm is the Lagrangian density of matter, and F(R, T) is an 
arbitrary function of R and T. Here and in this subsection, we use the unit of G = c = 1. 
As past related studies, a theory whose Lagrangian density is described by an arbitrary 
funct ion of R and the Lagr angi an density of matter as F(R,Lm) has been explored in 



Ref. 



138j . Moreover, in Ref. 139| a theory in which the cosmological constant is written by 



a function of the trace of the stress-energy tensor as A(T) has been investigated. 
From the action in Eq. (18.871) . the gravitational field equation is given by 

F R (R, T) - l -F{R, T)g, u + (g^D - V,V U ) F R {R, T) = 8vr (-F T (R, T)) % u - F T (R, T)6 



(8.88) 

with 6^ = g a P (STap/Sg^), which follows from the relation 5 (g af3 T a p/ '5g^) = % v + Q^, 
and F R (R, T) = dF(R, T)/dR, F T (R,T) = dF(R,T)/dT. The contraction of Eq. (ETggjl 
yields F R (R, T)R + 3DF R (R, T) - 2F(R, T) = (8tt - F T (R, T)) T - F T (R, T)Q with = 
g^Qfiu- Combining Eq. (I8.88P and the contracted equation and eliminating the \3F R (R, T) 
term from these equations, we find 



F R (R,T) (r^ - l -Rg,)j + ~F(R,T)g^ 



l7r-F T (R,T))[%„--Tg^ 



F T (R,T)[Q 



V»V V F R {R,T). (8.89) 



On the other hand, the covariant divergence of Eq. (18.871) as well as the energy-momentum 
conservation law V [F R (R, T) - (1/2) F(R, T)g^ u + (jg^D - V M V„) F R (R, T)\ = 0, which 
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corresponds to the divergence of the left-hand side of Eq. (18.871) . we acquire the divergence 
of % v as 

VM7 "^ = sf-fyfrT) [{%V + Q(lu) W M FT ^ T) + V " ^ ] ' (8 ' 90) 
In addition, from T^ u = g^ u Lyi — 2 (dL M / 'dg^ v ) we have 



ST aj 3 _ ( 5g a p 1 \ 1 d 2 L 



Using the relation 5g a p/5g tl1 ' = —ga P g^u^ v with 6f%, = 5g p<T /5g^ u , which follows from 
g ap g pp = 8P, we obtain 

8^ = -2%, + g„L M - 2g a ?g^^ ■ (8.92) 

Provided that matter is regarded as a perfect fluid, T^ u is expressed as T hLV = 
(Pm + Pvi) u^u v — Pug^u, where being the four velocity satisfying g^u^Uv = —1, and 
Pm and Pm are the energy density and pressure of the perfect fluid, respectively, we acquire 
Qfiu — -~2T ilv — P\ig P u- 

2. Example 

As an example, we examine the case that F(R,T) = R + 2iq(T) with Pi(T) being an 
arbitrary function of T and matter is a perfect fluid. In this case, Eq. f!8.88j) becomes 

G>„ = (svr + % v + (Vm^S + FiCT)) g, v , (8.93) 

where G^ v = R^ v — (1/2) Rg^ v is the Einstein tensor. In the flat FLRW background, for the 
matter to be a dust, i.e., Pm = 0, and Pi(T) = AT with A being a constant, the gravitational 
field equations are given by 

a \ ( otx 



, +A p M , (8.94) 
a I V 3 



2 

a la' 



2-+ - =Ap M . (8.95) 
a \a J 

From these equations, we have 

The solution of this equation is given by H = p/t with p = [2 (87r + 3A)] / [3 (8n + 2A)]. 
Thus, we obtain a = t p with p > 1, and consequently the accelerated expansion of the 
universe can be realized. 



93 



IX. f(T) GRAVITY 



In this section, we explore f(T) gravity 2 . It is known that as a candidate of an al- 
ternative gravitational theory to general relativity, there exists "teleparallelism" in which 
the Weitzenbock connection is used 3l|. In this theory, there is only torsion T and the 



curvature R defined by the Levi-Civita connection does not exist. Recently, to account 
for the l ate t ime accelerated expansion of the universe as well as inflation in the early 
universe 14(| . by extending the teleparallel Lagrangian density described by the torsion 



scalar T to a function of T as f(T) 32|, |33j, various studies in f(T) gravity have been ex- 



ecuted. This concept has the same origin as the idea of f(R) gravity. In order to examine 
whether f(T) gravity can be worthy of being an alternative theory of gravitation to gen- 



eral relativity, recent 
in the literature 58, 



y a number o f asp ects of f(T) gravity have widely been investi gated 



141 



147 



non-trivial conformal frames 
gularities 



147 



148l | . For examp le, the lo cal Lorentz invariance 1431 ] . 



58 



144f | , thermodynamics 1451 . Il46| , and finite-time future sin- 



1481 ] . In this review, we concentrate on the issues on the finite-time future 



singularities in f(T) gravity and review the results in Ref. 58]. 



A. Basic formalism and fundamental equations 

We use orthonormal tetrad components ca{x^) in the teleparallelism. At each point 
of the manifold, an index A runs over 0, 1,2,3 for the tangent space. The relation to the 
metric g^ v is given by = rjAB^t&vi where \i and v are coordinate indices on the manifold 
and run over 0, 1,2,3. Moreover, forms the tangent vector of the manifold. We define 
the torsion T p ,,„ and contorsion K pu „ tensors as 

flu p 

n^ e A(a/-^) (9.i) 

and 

K ,u p = _ 1 ^ _ T u^ _ ^ ^ ) (g 2) 

respectively. In general relativity, the Lagrangian density is described by the Ricci scalar 
R, whereas the teleparallel Lagrangian density is represented by the torsion scalar T. With 

2 For clarity, we use the notation U F(R)" gravity and "/(T)" gravity throughout this review. 
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the torsion tensor as well as the contorsion tensor, we first define the following quantity 

1 



S 



[IV 



( K^ v A- r)^ T au — r) u T a>1 ) 

\ p ' u p a p a) 



(9.3) 



By using this quantity in Eq. (19. 3p . the torsion scalar T is described as 



T = SrT> 



(9.4) 



The modified teleparallel action of f(T) gravity with the matter Lagrangian £ M is expressed 



as 



33] 



d x\e\ 



f(T) 
2k 2 



M 



(9.5) 



Here, |e| = det (e^) = y/—g. By varying the action in Eq. ( 4.51) with respect to the vierbein 



vector field e^, we acquire the gravitational field equation 



32] 



\d» {eSjD f ~ e X A T^ x S p »»f + S/ V d, (T) f + \e\f = ^7^/, (9.6) 

with T^ M ) p v being the energy-momentum tensor of all perfect fluids of ordinary matter such 
as radiation and non-relativistic matter. 

We assume the flat FLRW space-time with the metric ds 2 = dt 2 — a 2 (t) $^t=i 2 3 (dx 1 ) 23 . 
Therefore, we have g^ iu = diag(l, —a 2 , —a 2 , —a 2 ) and the tetrad components e p = (1, a, a, a). 
From these relations, we find T = —6H 2 . In this flat FLRW universe, we can write the 
gravitational field equations in the same forms as those in general relativity 

(9.7) 

Pde) , (9.8) 

(9.9) 

(9.10) 
(9.11) 
(9.12) 

Here, F = df/dT and F' = dF/dT. Moreover, we express the energy density and pressure 
of all perfect fluids of generic matter as pm and Pm, respectively. These perfect fluids satisfy 



H 2 = 


K 2 

y (pM + Pde) , 




K 2 


H = 


\Pm + Pm + Pde 






Pde = 


2V 1 ' 


-Pde = 


~ (4J 2 + Ji) , 


Ji = 


-T-/ + 2TF, 


J 2 = 


(1 — F — 2TF') H . 



In this section, the metric signature of (+, — , — , — ) is adopted. 
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the continuity equation pM + 3H (pu + Pvi) = 0. In addition, for the representations of 
Pde in Eq. (I9.9p and Pde in Eq. (19.101) . the standard continuity equation can be met as 
p DE + 3H (p DE + P DE ) = 0. 



B. Reconstruction of f(T) gravity with realizing the finite-time future singularities 

1. Finite-time future singularities in f(T) gravity 



We provided that an expression of the Hubble parameter [54J realizing the finite-time 
future singularities and the resultant scale factor obtained from the form of the Hubble 
parameter are given by 



(t 8 ~ty 
(t s -ty 



H~ u \ )q , forg>0, (9.13) 
H ~ H s + ^ ks , for q < -1 , -1 < q < 0, (9.14) 



a ~ a s exp 



q-1 



(t s - t) 



for 0< q < 1, 1 < q, (9.15) 



a ~ a& (+ for g = 1 ' ^ 9 ' 16 ^ 

with /i s (> 0), H s (> 0) and a s (> 0) being positive constants, q(^ 0, —1) a non-zero constant, 
and t s the time when the finite-time future singularity appears. We only deal with the period 
< t < t s because H should be a real number. For the expression of H in Eqs. (I9.13P and 
(I9.14p . the finite-time future singularities occur in the case of each range of the value of 
q. It follows from p DE w p cff = 3H 2 /k 2 in Eq. P DE w P cff = - + 3# 2 ) /k 2 in 

Eq. (19.1 Op . and the relation T = —QH 2 that in the limit of t — > t s , if PT — > oo, both pde and 
Pde diverge; if H becomes finite but H does infinite, Pde diverges, although p DE does not. 
We note that J\ in Eq. (19. lip only depends on T, i.e., H and J 2 in Eq. (19.121) is proportional 
to H. In Table HJ we summarize the conditions to produce the finite-time future singularities 
in the limit of t — > t s . 

2. Reconstruction of an f(T) gravity model 

Next, we reconstruct an f(T) gravity model with realizing the finite-time future singu- 
larities in the limit of t — > t s . From Eqs. (19.91) and (19.101) . the EoS of dark energy is given 
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TABLE I: Conditions to produce the finite-time future singularities in the limit of t — > t s . 



Type q& 0, -I) a H H Pde ^de 

I q > 1 o^oo ff^oo ij^oo J 17 LQ J 17 LQ 



or J 2 ^ 



III 


< q < 1 


a - 


■> a s 


iJ - 


■> oo 


H - 


■> oo 




II 


-1< g < 


a - 


■> a s 


F - 


*H S 


H - 


■> oo 




IV 


q<-l 


a - 


■> a s 


F - 


*H S 


H - 


■> 





(<7 integer) (Higher 



derivatives 
of i/ diverge.) 



by 



DE 



4 (1 - F - 2TF 1 ) H + (-T - / + 2TF) 



Pde ~T-f + 2TF 

This expression can be rewritten to a fluid description explained in Sec. Ill A as 



Pde = -Pde + J(H, H) , 
J = -X \2(1-F- 2TF) H 



(9.17) 

(9.18) 
(9.19) 



Here, J corresponds to —/(/?) in Eq. ( 13. 2p . The comparison of Eq. (19 . 1 8[) with P c g = 
-Pes - 2H/k 2 leads to H + (k 2 /2) J(H, H) = 0. By combining Eq. and this equation, 

we obtain H (F + 2TF') = 0. This yields the condition F + 2TF = because H ^ for H 
in Eqs. f l9TT3j) and (j9TTi]l . 

On the other hand, Eqs. (19. 7p and (19. 8p can be reduced to 



- / + 2TF = 
F + 2TF' = 



(9.20) 
(9.21) 



where we have also used Eqs. fl9T9|) - (l9~T2|) and H ^ for H in Eqs. (I9~T3|) and (l9~T4p . We 
see that Eq. (19.211) is equivalent to the above condition and Eq. (19.201) corresponds to a 
consistency condition. For a power-law model given by 



f(T) = AT 
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(9.22) 



with A(y^ 0) and 0) being non-zero constants, which with A = 1 and a = 1 corresponds 
to general relativity, Eq. ( EOT]) becomes F + 2TF' = A (-6)"" 1 (2a - 1) F 2 ^" 1 ) = 0, and 
Eq. ( E2D} reads -f + 2TF = A (-6)° (2a - 1) # 2q = 0. In the limit of t -> t s , both of these 
equations have to be satisfied asymptotically. If a = 1/2, these equations are always met. 
Hence, for q > 0, a < 0, while for g < 0, a = 1/2. We remark that these are not sufficient 
but necessary conditions to realize the finite-time future singularities. In fact, if a < 0, 
the Type I singularity occurs faster than the Type III and finally the Type I singularity 
happens because the speed of the divergence depends on the absolute value of q. From the 
same reason, if a = 1/2, the Type IV singularity appears faster than the Type II singularity 
and eventually the Type IV singularity occurs. We also note that for an exponential model 
f{T) = C e exp (A e T) with C e (^ 0) and A c (^ 0) being non-zero constants and an logarithmic 
model f(T) = -Diln(7iT) with D\(^ 0) being a non-zero constant and 7i(> 0) being a 
positive constant, both Eqs. ( 19 . 1 3[) and (I9.14p cannot be simultaneously met, and therefore 
these models cannot produce the finite-time future singularities. 



For "u>" singularity, the scale factor is given by 61] 




n(t) = a. | 1 - ^ | — I^ll-^ •f | <! 



1 - 2/ (3a) t 
n-2/(3a) T s 

VW na , L _ L f^MMr'V' (L) 2 " 3 " . (9 . 23) 



n — 2/ (3<t) \ 3a { n — 1 J J \t 



Here, a and n are arbitrary constants. When t — > t s , both H and H becomes zero, whereas 
the effective EoS for the universe w c r = (1/3) (2c/dec — 1) — >■ oo with gdec = —da/a 2 being the 
deceleration parameter. Thus, in the above limit of t — > t s , we find F + 2TF' = because 
H{t -> t s ) = 0. For a power-law model in Eq. (I9~22|) with A^0 and a > 1, Eq. f l9~2T|) can 
asymptotically be met due to H(t — > t s ) = 0, and hence 'W singularity can occur. 

3. Removing the finite-time future singularities 

We examine the possibility to remove the finite-time future singularities by taking a 
power-law correction term f c (T), given by 



f c (T) = BT> b . (9.24) 
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TABLE II: Necessary conditions for the appearance of the finite-time future singularities on a 
power-law f(T) model in Eq. f|9.22|) and those for the removal of the finite-time future singularities 
on a power-law correction term f c (T) = BT@ in Eq. (|9,24p . 



Type 


?(^0, -1) 


Final appearance 


f(T) = AT a 


/ C (T) = BT? 


I 


?>1 


Occur 


a < 


13 > 1 


III 


< q < 1 




a < 


P > 1 


II 


-1< q < 




a = 1/2 


/3/1/2 


IV 


q < — 1 (q ^ integer) 


Occur 


a = 1/2 


/3/1/2 



Here, 0) and /3(^ 0) are non-zero constants. By plugging the total form of 

f(T) = AT a + BT 13 into Eqs. (ICTl) and (ICTj) . we obtain -/ + 2TF = A (2a — 1) T a + 
B (2/3 - 1) T' 3 ^ and —F - 2TF' = -Aa (2a - 1) T^ 1 - B/3 (2/3 - 1) T 13 ' 1 ^ 0. From 
the investigations in Sec. IX B 2, the latter inequality is satisfied when the condition that 
for q > 0, /3 > 0, whereas for g<0,/3^1/2, is met. Accordingly, if /3 > 1, in the limit of 
t — > t s both ( 19. 7p and ( 19. 8 p cannot be met. Thus, the finite-time future singularities in f(T) 
gravity can be removed by a power-low correction as T 13 , where /3 > 1. It is remarkable to 
mention that a T 2 term can cure all the four typ es of the finite-time future singularities in 
f(T) gravity, similar to that in F(R) gravity 12[. In Table Hi] we show necessary conditions 
for a power-law f(T) model in Eq. (I9.22|) to prodece the finite-time future singularities and 
those appearance, and necessary conditions for a power-law correction term f c (T) = BT 13 
in Eq. (I9.24p to remove the finite-time future singularities. 

We also remark that In terms of the 'W singularity, if a power-law correction term in 
Eq. ( I9.24p with B ^ and /3 < is taken, the gravitational field equations ( 19. 7p and (19. 8p 
cannot be met asymptotically. As a consequence, the power-law correction term can cure 
the "w" singularity. 
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TABLE III: Forms of H and f(T) with realizing (i) inflation, (ii) the ACDM model, (iii) Little Rip 
scenario and (iv) Pseudo-Rip cosmology. 



Cosmology a H f(T) 

(i) Power-law inflation a = a inf i hinf H = h inf /t , f{T) = AT a , 

[when t — > 0] ai n f > h 1Ti t > 1 a < or a = 1/2 



(ii) ACDM model a = a A exp (H A t) , H = VA/3 f(T) = T - 2 A , 
or exponential inflation oa > = constant A > 

A > 

(iii) Little Rip scenario a = olr H = H^exp (£t) , f(T) = AT a , 
[when t — > oo] x exp \(H LR /^) exp (£<)] , i/ LR > and ^>0 a < or a = 1/2 

olr > 

(iv) Pseudo-Rip cosmology a = apRCOsh (t/to) , H = ffpRtanh (t/to) j f{T) = As/T 

a P R > .ffpR > 



C. Reconstructed f(T) models performing various cosmologies 

Furthermore, we describe the reconstructed f(T) models in which the following various 
cosmologies are realized: (i) inflation, (ii) the ACDM model, (iii) Little Rip scenario and 
(iv) Pseudo-Rip cosmology. We present expressions of a, H and f(T) realizing the above 
cosmologies in Table HTT1 Here, a in f(> 0) and olr(> 0) are positive constants and h in {(> 1) 
is a constant larger than unity. We note that the form of /(T) and the conditions for it are 
derived so that the gravitational field equations ( 19. 71) and ( 19. 81) . 

As another quantity to show the deviation of a dark energy model from the ACDM model, 
in addition to the EoS u>de of dark energy in Eq. ( 13. 9ft . the deceleration parameter q^ec hi 
Eq. (I3.97p . and the jerk parameter j in Eq. (I3.98p . the snark parameter s is used, which 



defined as 



80| 



(9.25) 



" 3 (g dcc - 1/2) 

For the ACDM model, u>de = — 1, <?dcc = — 1, j = 1 and s — 0. Hence, by examining the 
deviations of (iode, Qdec, j, s) from (—1, —1, 1,0) for the ACDM model and using these four 
parameters as a tool of observational tests, we can distinguish a dark energy model from the 
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TABLE IV: Expressions of u>de> Idee, 3 and s at the present time t = to for the ACDM model, 
Little Rip scenario and Pseudo-Rip cosmology [581 ] . 



Model ^de(o) <?dcc(o) Jo so 

ACDM -1-1 1 

model 

Little Rip -l-(2/3)x, -1 - x 1 + X (x + 3) - [2* (* + 3)] 

scenario x = #o/ {H LR e) x [3 (2* + 3)] _1 

< 0.36 , 

e = 2.71828 

Pseudo-Rip -1 - [25/ (3s 2 )] , -1 + (5 2 t 2 - l) / (5 2 t 2 ) , 1 + (l - 5 3 t 2 ) / (<5 3 t 2 ) , [2/ (35)] 
cosmology x (5 3 t 2 — l) 

x (<5 2 t 2 + 2)" 1 

5 = H /H PR s 2 = sinh 2 1 = 1.38 t 2 = tanh 2 1 = 0.580 

< 0.497 



ACDM model. In Table HVl we display the expressions of u>de, <7dcc, j and s at the 
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present 
. Here, 



time to for the ACDM model, Little Rip scenario and Pseudo-Rip cosmology 
s = s(t = t ). 

Finally, we mention another feature of f(T) gravity. It has been discuss ed t hat in the 
star collapse, the time- dependent matter ins tabi lity found in F(R) gravity 1491 ] . which is 



related to the well-studied matter instability 



in the relativistic star formation process 
gravity 



50| leading to the appearance of a singularity 



1511 ] . can also happen in the framework of f(T) 



D. Thermodynamics in f(T) gravity 



In this section, to explore whether f(T) gravity is worthy of an alternative gravitational 
theory to general relativity, we investigate thermodynamics in f(T) gravity. In particular, 



the second law of thermodynamics around t 
by applying the procedure proposed in Refs. 



re finite- time future singularities is st udiec 



146 



1521 ] . Black hole thermodynamics 153 1 
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suggested the fundamental relation of gravitation to thermodynamics (for recent reviews, see, 



e.g. 



154j ) . With the proportionality of the entropy to the horizon area, in general rela tivity 



the Einstein equation was obtained from the Clausius relation in thermod ynam ics 



This consideration has been extended to more general gravitational theories 



156 



155]. 



157]. 



1. First law of thermodynamics 



It is known in 146l . Il52l . |X58| ] that when the continuity equation of dark component is met 
as /?DE + 3-ff (pde + -Pde) = 0, we can have an equilibrium description of thermodynamics. In 
the flat FLRW universe with the metric ds 2 = h a/ sdx a dx 13 + f 2 dQ 2 , where r = a (t)r, x° = t 
and x 1 = r with the two-dimensional metric h a p = diag(l, — a 2 (t)), dQ 2 is the metric of 
two-dimensional sphere with unit radius. The radius f ^ of the apparent horizon is described 
by ta = 1/H. The relation h al3 d a fdpf = leads to the dynamical apparent horizon. The 
time derivative of ta = 1/H yields — df^/f^ = HHdt. Combining the Friedmann equation 
(19. 7p with this equation presents [1/ (4tiG)} dtA = f\H (p t + P t ) dt with p t = p DE + p M and 
P t = Pde + Pm being the total energy density and pressure of the universe, respectively. 
The Bekenstein-Hawking horizon entropy in general rela tivity is written by S = A/ (4G). 
Here, A = 47rf A is the area of the apparent horizon 1531 ] . Thus, with the horizon entropy 
as well as the above relation, we obtain [1/ (2^7^)] dS = 47rr\H (p t + P t ) dt. The Hawking 
temperature Th = |ft S g|/ (27r) is considered to be the associated temperature of the apparent 
horizon, and the surface gravity K sg is given by 159 ] 



1 



K 



--d a 



-hh a0 d, 



ta 



[G r A (l-3w t )p t . (9.26) 



2V-h " V "J f A \ 2HrAj 3F 

Here, h is the determinant of the metric h a p and w t = Pt/ Pt is the EoS for the total of 
energy and matter in the universe. From Eq. (I9.26p . we see that if w t < 1/3, n sg < 0. The 
substitution of Eq. ( 19. 26ft into Th = \n S g\/ (27r) yields 

1 



Tn 



r A 



(9.27) 



2nr A \~ 2Hf A 

By combining the above relation of the horizon entropy S with Eq. (19.271) . we find 
TudS = Anr A H (p t + P t ) dt — 2i\f 2 A (p t + P t ) dr~A- Moreover, the Misner-Sharp energy 16o| 
is expressed by E = ta/ (2G) = V pt, where V = 47rf^/3 is the volume inside the apparent 
horizon. From the second equality, we see that E corresponds to the total intrinsic energy. 
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With this eq uatio n, we obtain dE 



work density 



161] is defined by W 



Anr\H (p t + P t ) dt + 4irr A p t dr a- In addition, the 
- (1/2) {T^Kp + T^^hap) = (1/2) (p t - P t ) 



with T^ a ? and T( DE ) a/3 being the energy-momentum tensor of matter and that of dark 
components, respectively. We plug the work density W into the relation on dE derived 
above, so that we can represent the first law of equilibrium thermodynamics as 



T K dS = -dE + WdV . 



(9.28) 



Thus, we acquire an equilibrium description of thermodynamics. It follows from the 
gravitational equations (19. 7p and ( 19. 8 j) as well as the above relation on dS that S = 



M 2 m\( Pt + P t ) = (&k/ G) ( T/T 



-(2tt/G) H/(3H 



> 0. Accordingly, for the 



expanding universe (H > 0), if the null energy condition p t + P t > in Eq. ( 13 . 1 0[) is 
satisfied, namely, H < 0, S always becomes large. 



2. Second law of thermodynamics 

We then explore the second law of thermodynamics. The Gibbs equation of all the matter 
and energy fluid is expressed as T K dS t = d (p t V)+P t dV = Vdp t + (p t + P t ) dV . We can write 
the second law of thermodynamics as dS sum /dt = dS/dt + dS t /dt > 0. Here, S sum = S + S t 
with St being the entropy of total energy inside th e ho rizon. If the temperature of the 
universe is the same as that of the apparent horizon 1621 ]. this can be represented as 

(9.29) 



2 

dS 

sum 

6n ( T \ 1 




dt G \ T ] AHT + T 



where we have used V = 4irr A /3, Eqs. (I9.8|) . (I9.27P and the relation on S shown in 
the last part of Sec. IX D 1. Consequently, the condition that Y = — (^AHT + T 



12H (2H 2 + H) > 



146| is met, the second law of thermodynamics can be verified. When 
t t s , for all of the four types of the finite-time future singularities in Table [H the relation 
2H 2 + H > is always satisfied. Since H > for the expanding universe, the second law of 
thermodynamics described by Eq. (I9.29p can be met around the finite-time future singular- 
ities including in the phantom phase (H > 0), although at the exact time of the appearance 
of singularity of t = t s this classical description of thermodynamics could not be applicable. 
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X. TESTING DARK ENERGY AND ALTERNATIVE GRAVITY BY COSMOG- 
RAPHY: GENERALITIES 



Next, we move to the comparison of the theoretical studies on dark energy and modified 
gravity with the observational data. In this section, we introduce the idea and concept 
of cosmography to observationally test dark energy and alternative gravitational theory to 
general relativity. 

The observed accelerated expansion of the cosmic fluid can be faced in several equivalent 
ways. In other words, both dark energy models and modified gravity theories seem to be 
in agreement with data. As a consequence, unless higher precision probes of the expansion 
rate and the growth of structure will be available, these two rival approaches could not 
be discriminated. This confusion about the theoretical background suggests that a more 
conservative approach to the problem of the cosmic acceleration, relying on as less model 
dependent quantitie s as possible, is welcome. A possible solution could be to come back 
to the cosmography 1631 ] rather than finding out solutions of the Friedmann equations and 
testing them. Being only related to the derivatives of the scale factor, the cosmographic 
parameters make it possible to fit the data on the distance - redshift relation without any a 
priori assumption on the underlying cosmological model: in this case, the only assumption 
is that the metric is the FLRW one (and hence not relying on the solution of cosmological 
equations). Almost eighty years after Hubble's discovery of the expansion of the universe, 
we can now extend, in principle, cosmography well beyond the search for the value of the 
only Hubble constant. The SNela Hubble diagram extends up to z = 1.7 thus invoking the 
need for, at least, a fifth order Taylor expansion of the scale factor in order to give a reliable 
approximation of the distance- redshift relation. As a consequence, it could be, in principle, 
possible to estimate up to five cosmographic parameters, although the still too small data 
set available does not allow to get a precise and realistic determination of all of them. 

Once these quantities have been determined, one could use them to put constraints on 
the models. In a sense, we can revert to the usual approach, consisting with deriving the 
cosmographic parameters as a sort of byproduct of an assumed theory. Here, we follow 
the other way of expressing the quantities characterizing the model as a function of the 
cosmographic parameters. Such a program is pa rticularly suited for the study of alternative 



theories as F(R) or f(T) gravity 



147 



164, 



165| and any equivalent description of dynamics 
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by effective scalar fields. As it is well known, the mathematical difficulties in analyzing 
the solution of field equations make it quite problematic to find out analytical expressions 
for the scale factor and hence predict the values of the cosmographic parameters. A key 
role in F(R) gravity and f(T) gravity is played by the choice of the function. Under quite 
general hypotheses, it is possible to derive relations between cosmographic parameters and 
the present time values of F^ n \R) = d n F/dR n or f {n) (T) = d n f/dT n , with n = 0, . . . , 3, 
whatever F(R), f(T) or their equivalent scalar-field descriptions are. 

Once the cosmographic parameters are determined, the method allows to investigate the 
cosmography of alternative theories matching with observational data. 

It is worth stressing that the definition of the cosmographic parameters only relies on the 
assumption of the FLRW metric. As such, it is however difficult to state a priori to what 
extent the fifth order expansion provides an accurate enough description of the quantities 
of interest. Actually, the number of cosmographic parameters to be used depends on the 
problem one is interested in. 

To illustrate the method, one can be concerned only with the SNela Hubble diagram so 
that one has to check that the distance modulus Hcp{z) obtained using the fifth order expan- 
sion of the scale factor is the same (within the errors) as the one /^de(^) of the underlying 
physical model. Being such a model of course unknown, one can adopt a phenomenological 
parameterization for the dark energy EoS and look at the percentage deviation A/^/^de as 
a function of the EoS parameters. Note that one can always use a phenomenologica l da rk 
energy model to get a reliable estimate of the scale factor evolution (see, for example, 166| ). 

Here, we will carry out such an approach using the so called Chevallier-Polarski-Linder 
(CPL) model [3ja, |36j, introduced below, and verified that A/^//^ DE is an increasing function 
of z (as expected), but still remains smaller than 2% up to z ~ 2 over a wide range of 
the CPL parameter space. On the other hand, halting the Taylor expansion to a lower 
order may introduce significant deviation for z > 1 that can potentially bias the analysis 
if the measurement errors are as small as those predicted by future observational surveys. 
However, the fifth order expansion is both sufficient to get an accurate distance modulus 
ov er the red shift range probed by SNela and necessary to avoid dangerous biases. As shown 
in 
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1681 ] . the method can highly be improved by adopting BAO and Gamma Ray Bursts 
(GRBs) as cosmic indicators. 

As stated above, the key rule in cosmography is the Taylor series expansion of the scale 
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factor with respect to the cosmic time. To this aim, it is convenient to introduce the following 
functions: 

H{t) = +±%, (10.1) 
a at 

, . 1 d 2 a 1 

^ s -SdFfl5- (10 ' 2) 
^ = < 10 ' 3 > 

. . 1 d A a 1 ,„ „ lS 

(10 ' 5) 

which are usually referred to as the Hubble, deceleration, jerk, snap, and /erA; parameters, 
respectively. (Here, for clear understanding, we have again defined the Hubble, deceleration 
(as q, which is the same as gdec in Eq. ( 13.971) ). jerk and snap parameters.) It is then a matter 
of algebra to demonstrate the following useful relations : 

H = -H 2 (l + q), (10.6) 

H = H%j + 3q + 2) , (10.7) 

d 3 H/dt 3 = H 4 [s- 4j - 3q(q + 4) - 6] , (10.8) 

d 4 H/dt 4 = H 5 [I - 5s + 10(q + 2)j + 30(g + 2)q + 24] . (10.9) 

Equations (I10.6l) -( fl0.9l) make it possible to relate the derivative of the Hubble parameter to 
the other cosmographic parameters. The distance - redshift relation m ay then b e obtained, 



starting from the Taylor expansion of a(t) along the lines described in 169l - ll71 |. 



By these definitions, the series expansion to the 5th order in time of the scale factor is 
a{t) = a(t ) jl + H (t - t ) - |tf 2 (t - t ) 2 + §# 3 (* - *o) 3 + 

+ |^-^o) 4 + |^ 5 (t-to) 5 +0[(t-t ) 6 ]} , (10.10) 



from which, we find 



^ = 1 + H (t- to) - 2 (t - to) 2 + §if 3 (t - t ) 3 + 

+ ^Ho(t ~ toT + - t ) 5 + 0[(t - t ) 6 ] . (10.11) 
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It is easy to see that Eq. ( 110. lip is the inverse of redshift z, being the redshift defined by 

a(t ) 



l + z 



a(t) 



The physical distance travelled by a photon that is emitted at time and absorbed at the 
current epoch t is 

D = c J dt = c(t — t*) , 

where c is the speed of light. Assuming t* = to — ® and inserting it into Eq. ( 110. lip , we 
have 

a(t ) 



l + z 



a(t - £ 



1 



(10.12) 



The inverse of this expression becomes 

i + ,=i + > + ( 1+ f)(^)V + u +? „ 



Jo\ 



6 / V c 



L> 3 + 



+ 1 + 2qn + -On H 

1 yo 4 y ° 6 2 12 



1 — 1 



10.13) 



Then, we reverse the series z(D) — > D(z) to have the physical distance D expressed as a 
function of the redshift z 



z(D) = Z' D 
+ Z 



i / -"o 



D 



? 2 







F 3 



H n D\ _ fH n D\ l _, fH D\ 3 , ^ 4 ^# £ x 



(10.14) 



with 



Z 1 



z z 

D 



J7 4 



1 + 



'7o 



i + % + , 

6 

1 + 2g + \ql 



Jo _ £o 
3 24 ' 
gojo Jo 
6 2 



12 



(10.15) 
(10.16) 

(10.17) 

(10.18) 

(10.19) 
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From this, we obtain 
D(z) 



cz 



{V° z +V]z + V\ z 2 + V\ z 3 + V\ z 4 + 0(z 5 ) } 



with 




v° z = 




v\ = 


-(*+?)■ 


v\ = 


l + .o + f- 


v\ = 






1 + 2q + 3q 2 



5 . s 



12 



21 



So 3 



7 



1 



2?o - pojo ~ ^oJo ~ |<foSo ~ Jo + y| 



fo 
6 



120 



In standard applications, other quantities can result in become useful 
• the luminosity distance: 



d L = — m ( a (*o)r 0; 



the angular-diameter distance: 



where tq(D) is given by 



a(t - 

= 77 / (a(*o)?"o) 

a (to J 



>i -f o(t) 



K = 0; 



(10.20) 



(10.21) 
(10.22) 

'10.23) 

(10.24) 

(10.25) 



10.26) 



10.27) 



10.2* 



If we consider the expansion for short distances, namely, if we insert the series expansion of 
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a(t) into r (D), we find 

"*° c dt f to c dt 



ro(D) 



r° c dt r° c dt r TT . . / 

/ -7tt= / \l + H (t -t) + (l 

Jt -R a(t) J to _n a I V 



go 
2 



1 + qo 



Jo 
6 



H 2 (t -t) 2 + 
H*(to - tf- 



l + 2 qo+ 6 - q 2 + q -f + 3 j-^ + l )H^t -tf + 0[(t -tr] 



D 



l#n£> 



«o 



2 c 



2 + go 
6 



#0^ 



+ 



6 + 6g + jo 
24 



+ 
+ 



HnD 



+ 



24 + 36g + 6q 2 + 8j - s 
120 



12 + 24g + 9q 2 + 2g j + 6j ~ *o + 12/ c 



72 



+ 







(10.29) 



To convert from physical distance travelled to r coordinate, we have to consider that the 
Taylor series expansion of sin-sinh functions is 



r (D) 



to 



*° c dt 

D 



a(t) 



k 



t 



to- 



c dt 
d a(t) 



+ 



to 



to- 



c dt 
_d a{t) 



-l 5 N 



so that Eg. (HO. lip with the spatial curvature K term becomes 



r (D) = °\k° d + ^ HoD 



a 



D~ 



+ n 



D 



HnD 



+ 7Z 



D 



HnD 



HnD 







HnD 



HnD 



+ 



with 



(10.30) 



(10.31) 



D 



D 



V 2 
IX. D 



1 , 
1 
2 
1 

6 
1 

24 
1 

120 
1 

144 



2 + g - 



Kc 2 

P/2 2 
- n u 



Kc 2 

6 + 6g + Jo-6— 

ii a OJ 

9 5Kc 2 (7+ 2q o y 
24 + 36g + 6g 2 + 8 Jo - s \ m> 

24 + 48g + 18g 2 + 4g jo + 12j - 2s + 24/ c 



3Xc 2 (15 + 10g +Jo) 



n 2 H 2 



(10.32) 
(10.33) 

(10.34) 
(10.35) 
(10.36) 
(10.37) 
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Using these definitions for luminosity distance, we acquire 

cz 



d L (z) 



He 



{V° L + V\ z + V\ z 2 + V\ z 3 + V\ z 4 + 0{z 5 )} 



(10.36 



with 



"2 (-! + «!) 

If 2 Kc 2 

-g l-go-3g 0+Jo + M 



a 2 
o"o 



2 - 2q - 15q 2 - 15q 3 + 5j + 10g jo + + 



1 

24 
1 

120 

97 . , in . 2 „ , 5A'c 2 (l + 8g + 9g 2 -2j ) 
27j + 10j -llso-to 



2Ac 2 (l + 3g ) 



[-6 + 6q + 81g 2 + 165g 3 + 105g 4 - 110g j - 105g 2 j - 15g s + 



a H 



While, for the angular diameter distance we find 



(10.39) 
(10.40) 

(10.41) 
(10.42) 
(10.43) 
(10.44) 



d A (z) 



cz 



{V° A + V A z + V 2 A z 2 + V 3 A z 3 + V\ z A + 0(z 5 )} 



(10.45) 



with 



r>2 



T) 4 



■ (3 + g ) , 

11 + 7q + 3ql - j 



Kc 2 

F72 2 
- n a J 



1 

24 



50 + 46g + 39g 2 + 15g 3 - 13j - 10g Jo - s 



2Kc 2 (5 + 3q ) 



H 2 n 2 
- n o u o 



120 



[274 + 326g + 411g 2 + 315g 3 + 105g 4 - 210g Jo - W5q 2 oJo - 15q s - 



137j'o + 10j 2 - 21s - k 



5Ac 2 (17 + 20go + 9g 2 -2j ) 

n 2 H 2 
a n Q 



(10.46) 
(10.47) 

(10.48) 
(10.49) 
(10.50) 
(10.51) 



We define Qq — 1 + lir^; which can be considered a purely cosmographic parameter, or 



f2o — 1 ~ = + ^r U; + &>x' ) where corresponds to the current fractional densities 



)(0) 



,(o) 



)(0) 



)(0) 
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of dark energy, if we explore the dynamics of the universe. With these parameters, we obtain 
V°r = 1, (10.52) 



1 



= ~2 (-3 + *) , (10.53) 

V\ y = ~ (12 - 5g + 3g 2 - j - fi„) , (10.54) 

2>L = ^ [ 52 " 20 <?o + 21g 2 " 15g 3 - 7 j + 10q o j o + s - 2O (l + 3g )] , (10.55) 

T>l y = ^ [359 - 184g + 186g 2 - 135g 3 + 105g 4 + 90g Jo - 105g 2 Jo - 15q s + (10.56) 

- 57j'o + 10j 2 + 9s -l - 5fi (17 - 6g + 9g 2 - 2j )] , (10.57) 



and 



V% y = l, (10.58) 



1 

2 
1 



T>\ v = —(! + «,) , (10.59) 



Z>a,„ = " g (-9o - 3g 2 + jo + Ho) , (10.60) 

T>% y = ~ (-2q + 3g 2 + 15g 3 - j - 10q o j - s + 2fi ) , (10.61) 

V% y = (1 - 6g + 9g 2 - 15g 3 - 105g 4 + 10q o3o + 105g 2 j + 15^0+ (10.62) 

- 3j - 10j 2 + s + l + 5fi ) • (10.63) 

Previous relations have been derived for any value of the curvature parameter. To il- 
lustrate the method, however, we can assume a spatially flat universe, using the simplified 
versions for K = 0. Now, since we are going to use supernovae data, it will be useful to 
give as well the Taylor series of the expansion of the luminosity distance at it enters the 
modulus distance, which is the quantity about which those observational data inform. The 
final expression for the modulus distance based on the Hubble free luminosity distance, 
fi(z) = 51og 10 d L (z), is 

H(z) = :— - • (log z + M x z + M 2 z 2 + M 3 z 3 + M 4 z 4 ) , (10.64) 
log 10 v ' 
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with 



M 1 
M 2 
M 3 



"2 (-! + «>) 



21 



(7-10g -9g 2 + 4j ) 



24 



(5 - 9g - 16g 2 - 10g 3 + 7j + 8g Jo + s ) , 
(-469 + 1004g + 2654g 2 + 3300^ + 1575g, 

24/ ) . 



2880 

2620g Jo - 1800g 2 j - 300g s - 324s 



200j 2 - 1148j'o+ 



(10.65) 
(10.66) 
(10.67) 

10.68) 



XI. AN EXAMPLE: TESTING F(R) GRAVITY BY COSMOGRAPHY 



The cosmographic approach can be used to deal with F(R) gravity 164J. However, similar 
considerations perfectly hold also for f(T) gravity 1471 ] or any scalar-tensor gravity model. 
In order to construct the cosmographic apparatus, we describe the Friedmann equation (13. 3p 
in the FLRW space-time in F(R) gravity as 

Pm 



F'(R) 



Pc 



11.1) 



where the prime denotes the derivative with respect to R, the gravitational coupling is taken 



as k 2 = 1 and p curv is the energy density of an effective curvature fluid: 



- (F(R) - RF'(R)) - 3HRF"(R) 



11.2) 



Assuming there is no interaction between the matter and the curvature terms (we are in the 
Jordan frame), the matter continuity equation gives the usual scaling pm = Pm(^ = ^o) a_3 — 
3HqQm(i~ 3 , with Qyj* the matter density parameter at the present time. The continuity 
equation for p curv then reads 

3H$n$RF"(R) 



+ 3if (1 + w cuvv )p 

curv 



with 



W c 



-1 



(F'(R)Y 

RF"(R) + R RF"'(R) - HF"(R) 



11.3) 



11.4) 



[F(R) - RF'(R)} /2 - 3HRF"(R) 
the barotropic factor of the curvature fluid. It is worth noticing that the curvature fluid 

quantities p curv and u> CU rv only depend on the form of F(R) and its derivatives up to the 

third order. As a consequence, considering only those current values (which may naively be 
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obtained by replacing R with Rq everywhere), two F(R) theories sharing the same values of 
F(Ro), F'(Ro), F"(Rq), F'"(Rq) will be degenerate from this point of view. One can argue 
that this is not strictly true because different F(R) theories will lead to different expansion 
rates H(t) and hence different current values of R and its derivatives. However, it is likely 
that two F(R) functions that exactly match with each other up to the third order derivative 
today will give rise to the same H(t) at least for t ~ to, so that (Rq, Rq, Rq) will be almost the 
same. Combining Eq. (111.3)) with Eq. (jll.ip . one finally gets the following master equation 
for the Hubble parameter 

zY L 

2F'(R) 

Expressing the scalar curvature R as function of the Hubble parameter as 



3H%d®aT 3 + RF"(R) + +R (RF'"(R) - HF"(R)\j . (11.5) 



R = -6[H + 2H 2 j (11.6) 

and inserting the resultant expression into Eq. flll.5p . one ends with a fourth order nonlinear 
differential equation for the scale factor a(t) that cannot easily be solved even for the simplest 
cases (for instance, F(R) oc R n ). Moreover, although technically feasible, a numerical 
solution of Eq. (II 1.5ft is plagued by the large uncertainties on the boundary conditions (i.e., 
the current values of the scale factor and its derivatives up to the third order) that have to 
be set to find out the scale factor. 

Motivated by these difficulties, we now approach the problem from a different viewpoint. 
Rather than choosing a parameterized expression for F(R) and then numerically solving 
Eq. ( jl 1.5ft for given values of the boundary conditions, we try to relate the current values 
of its derivatives to the cosmographic parameters (qo,jo, sq, Iq) so that constraining them in 
a model independent way can give us a hint for what kind of F(R) theory is able to fit the 
observed Hubble diagram. Note that a simila r ana lysis, but in the context of the energy 



conditions in F(R), has yet been presented in [172[. However, in that work, an expression 
for F(R) is given and then the snap parameter is computed in order for it to be compared to 
the observed one. On the contrary, our analysis does not depend on any assumed functional 
expression for F(R). 

As a preliminary step, it is worth considering again the constraint equation All .6ft . Dif- 
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ferentiating with respect to t, we easily get the following relations : 



R = -Q \ H + AHHJ , (11.7) 
R = -6 (d 3 H/dt 3 + AHH + 4i7 2 ) , (11.8) 
d 3 R/dt 3 = -6 (d 4 H/dt 4 + AHd 3 H/dt 3 + 12HH\ . (11.9) 

Evaluating these at the present time and using Eqs. fll0.6p - fll0.9p . we finally obtain 

R = -QH 2 (1 - q ) , (11.10) 

Ro = -QH 3 (j - q - 2) , (11.11) 

Ro = -6if 4 (s + q 2 Q + 8q + 6) , (11.12) 

d 3 i?oM 3 = -6i7 5 [l -s Q + 2(q + 4) Jo - 6(3g + 8)g - 24] , (11.13) 

which will turn out to be useful in the following. 

We come back to the expansion rate and master equations flll.ip and fill. 51) . Since they 
have to hold along the full evolutionary history of the universe, they are naively satisfied 
also at the present time. Accordingly, we may evaluate them in t = to and thus we easily 
obtain 

H ° = nRo) + ' (1L14) 



• _ 3i7 2 < RlF"\Ro) + [Ro - H Ro ) F"(R ) 

- H °-2nRV) + 2fv^) • (1L15) 

Using Eqs. flTDTgp — flTUrgp and flll.iup - flll.13p . we can rearrange Eqs. (I11.14p and fj!1.15j) as 
two relations among the Hubble constant Hq and the cosmographic parameters (qo,jo, so), 
on one hand, and the present day values of F(R) and its derivatives up to third order. 
However, two further relations are needed in order to close the system and determine the 
four unknown quantities F(R ), F'(Rq), F"(Rq), F"'(Ro). A first one may be easily obtained 
by noting that, inserting back the physical units, the rate expansion equation reads 

H 2 = J^y(p m + Pc Urv FW), (H.16) 

which clearly shows that, in F(R) gravity, the Newton's gravitational constant G (restored 
for the moment) is replaced by an effective (time dependent) G e s = G/F'(R). On the other 
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hand, it is reasonable to assume that the value of G e s at the present time is the same as 
that of the Newton's one, so that we can acquire the simple constraint 

G eS (z = 0) = G^ F'(Ro) = 1 . (11.17) 

In order to find the fourth relation we need to close the system, we first differentiate both 
sides of Eq. (I11.5P with respect to t. We thus obtain 



H 



R 2 F"'(R) + (A- HR) F"(R) + 3H^a~ 3 R 3 F^ (R) + (3RR - HR 2 ^j F"'(R) 



2^RF"(R)^j \F'(R)) 2 2F'(R) 
(d 3 R/dt 3 -HR + Hit} F"{R) - 9H^Ha- 3 



11.18) 



2F'(R) 

with FW(R) = d 4 F/dR 4 . We now suppose that F(R) may be well approximated by its 
third order Taylor expansion in terms of (R — R ), i.e., we set 

F{R) = F(Ro) + F'( J R )( J R - R Q ) + ^"{^{R - R ) 2 + \f"\Rq){R - R ) 3 . (11.19) 

2 o 

In such an approximation, we find F^(R) = d n F/R n = for n > 4, so that naively 
F^(Rq) = 0. Evaluating Eq. (111.181) at the present time, we acquire 

RqF"'(Rq) + (ji — HoRo^ F"(Rq) + 3ifo^M (^R R — HR^j F"'(Rq) 
2(r F"(R )Y 1 (F'(R )) 2 2F '™ 



(d 3 R /dt 3 - H R + HqRo) F"(R ) - 9H 3 n$ 

- • (11-20) 

2F'(R ) y ' 

Now, we can schematically proceed as follows. We evaluate Eqs. fll0.6j) (I10.9j) at z = and 
plug these relations into the left-hand sides of Eqs. ( 111.14}) . (111.151) and ( 11 1.20f) . Then, we 
insert Eqs. (Ill.10p - flll.13p into the right-hand sides of these same equations, so that only 
the quantities related to the cosmographic parameters (qo,jo, So, lo) and the F(R) term can 
enter both sides of these relations. Finally, we solve them under the constraint (111.1 7j) with 
respect to the current values of F(R) and its derivatives up to the third order. After some 
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algebra, we eventually end up with the desired result 

F(R ) Vo(qo,jo,So,lo)Sl M + Qo(go, jo, s , l ) 



F'(Ro) = 1 , 

F"(Rq) V2(qo, jo, s )tt M + Q 2 (go, jo, so) 
(QH^)' 1 ~ TZ(q ,j ,s ,l ) 

F'"(Ro) V 3 (q ,jo, s , k)^M + 63(90, Jo, «o, k) 



where we have defined 



(jo -9o - 2)n(q ,jo,s J ) 



V = (jo - go - 2)/ - (3s + 7 j + 6q 2 + 41g + 22)s 



- [(3g + 16)j + 20g 2 + 64g + 12] j - (H + 25g 3 + 96g 2 + 72q + 20) 
Qo = (g 2 - io?o + 2g )/ + [3g s + (4g + 6) j + 6gjj + 44g 2 + 22g - 12] s - 

+ [2j 2 + (3g 2 + 10g - 6) jo + 17g 3 + 52g 2 + 54g + 36] j + 

+ 3q 5 + 28g 1 + 118gg + 72g 2 - 76g - 64 , 
V 2 = 9s + 6j + 9g 2 + 66g + 42 , 

Q 2 = - {6(g + l)s + [2j - 2(1 - g )] Jo + 6g 3 + 50g 2 + 74g + 32} , 

V 3 = 3l + 3s - 9(g + 4) jo - (45g 2 + 78g + 12) , 

Q 3 = - {2(1 + q )l + 2(j + q )s - (2j + 4g 2 + 12g + 6)j + 

- (30g 3 + 84g 2 + 78g + 24)} , 
n = (jo - go - 2)/ - (3s - 2j + 6g 2 + 50g + 40)s + 

+ [(3g + 10) j + ll^o + 4g + 18] Jo - (3g£ + 34g 3 + 246g + 104) . 



(11.21) 
(11.22) 
(11.23) 

(11.24) 



(11.25) 



(11.26) 
(11.27) 
(11.28) 
(11.29) 

(11.30) 

(11.31) 



Equations f lll.2ip - flll.3ip make it possible to estimate the current values of F(R) and its first 
three derivatives as function of the Hubble constant Hq and the cosmographic parameters 
(go, Jo, s o, ^0) provided a value for the matter density parameter £1$ is given. This is a 
somewhat problematic point. Indeed, while the cosmographic parameters may be estimated 
in a model independent way, the fiducial value for fijjj is usually the outcome of fitting a 
given dataset in the framework of an assumed dark energy scenario. However, it is worth 



noting that all the different models converge on the concordance value Q 



(o) 

M 



0.25 which 



is also in agreement with astrophysical (model independent) estimates from the gas mass 
fraction in galaxy clusters. On the other hand, it has been propos ed that F (R) theories may 



avoid the need for dark matter in galaxies and galaxy clusters 



173 



- Il77| . In such a case, 
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the total matter content of the universe is essentially equal to the baryonic one. According 
to the primordial elements abundance and the standard Big bang nucl eosynthesis (BBN) 
scenario, we therefore find ~ u b /h 2 with u b = VL^h 2 ~ 0.0214 [178| and h the Hubble 



constant in units of lOOkm/s/Mpc. Setting h = 0.72 in agreement with the results of the 



HST Key project 68|, we hence obtain ilw = 0.041 for a universe consisting of baryons only. 



M 

Thus, in the following we will consider both cases when numerical estimates are needed. 

It is worth noticing that H only plays the role of a scaling parameter giving the correct 
physical dimensions to F(R) and its derivatives. As such, it is not surprising that we need 
four cosmographic parameters, namely (go, jo, so, ^o), to fix the four F(R) related quantities 
F{Rq), F'{Rq), F"(Ro), F"'(R ). It is also worth stressing that Eqs. fflX2TT) -f TTl~Ml) are 
linear in the F(R) quantities, so that (go, Jo, so, lo) can uniquely determine the former ones. 
On the contrary, inverting them to acquire the cosmographic parameters as a function of the 
F(R) ones, we do not obtain linear relations. Indeed, the field equations in F(R) theories 
are nonlinear fourth order differential equations in terms of the scale factor a(t), so that 
fixing the derivatives of F(R) up to the third order can make it possible to find out a class 
of solutions, not a single one. Each one of these solutions will be characterized by a different 
set of cosmographic parameters. This explains why the inversion of Eqs. flll.2ip -f lll.3ip 
does not give a unique result for (go, Jo, s o, M- 

As a final comment, we again investigate the underlying assumptions leading to the above 
derived relations. While Eqs. ( 111. 141) and (111.151) are exact relations deriving from a rigorous 
application of the field equations, Eq. flll.20j) heavily relies on having approximated F(R) 
with its third order Taylor expansion ( 111.191) . If this assumption fails, the system should not 
be closed because a fifth unknown parameter enters the game, namely F^(Rq). Actually, 
replacing F(R) with its Taylor expansion is not possible for all the class of F(R) theories. 
As such, the above results only hold in those cases where such an expansion is possible. 
Moreover, by truncating the expansion to the third order, we are implicitly assuming that 
higher order terms are negligible over the redshift range probed by the data. That is to say, 
we are assuming that 

F {n) {R ){R-R ) n «J2 Y^(R-R ) m forn>4 (11.32) 



ml 

m=0 



over the redshift range probed by the data. Checking the validity of this assumption is not 
possible without explicitly solving the field equations, but we can guess an order of magnitude 
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estimate considering that, for all the viable models, the background dynamics should not 
differ too much from the ACDM one at least up to z ~ 2. Using the expression of H(z) for 
the ACDM model, it is easily to see that R/Ro is a quickly increasing function of the redshift, 
so that, in order Eq. (jll.32p should hold, we have to assume that F^(Rq) « F"'(Rq) for 
n > 4. This condition is easier to check for many analytical F(R) models. 

Once such a relation is verified, we have to still worry about Eq. (111.171) relying on the 
assumption that the cosmological gravitational constant is exactly the same as the local 
one. Although reasonable, this requirement is not absolutely demonstrated. Actually, the 
numerical value usually adopted for the Newton's constant G is obtained from laboratory 
experiments in settings that can hardly be considered homogenous and isotropic. Similarly, 
the space-time metric in such conditions has nothing to do with the cosmological one, so 
that strictly speaking, matching the two values of G should be an extrapolation. Although 
commonly accepted and quite reasonable, the condition Gi oca i = G cosmo could (at least, 
in principle) be violated, so that Eq. ( 111.171) could be reconsidered. Indeed, as we will 
see, the condition F'(Rq) = 1 may not be verified for some popular F(R) models which 
has recently been proposed in the literature. However, it is reasonable to assume that 
G e ff(z = 0) = G(l + e) with e « 1. When this be the case, we should repeat the 
derivation of Eqs. (111.211) — (111. 24[) by using the condition F'(R ) = (1 + e)" 1 . By executing 
the Taylor expansion of the results in terms of e to the first order and comparing with the 
above derived equations, we can estimate the error induced by our assumption e — 0. The 
resultant expressions are too lengthy to be reported and depend on the values of the matter 
density parameter the cosmographic parameters {qo, jo, soJo) and £ in a complicated 
way. Nevertheless, we have numerically checked that the error induced on F(R Q ), F"(R ), 
F"'(Rq) are much lower than 10% for the valu e of e as high as an unrealistic e ~ 0.1. 



However, results are reliable also for these cases 



164]. 



A. The CPL model 

A determination of F(R) and its derivatives in terms of the cosmographic parameters need 
for an estimate of these latter from the data in a model independent way. Unfortunately, 
even in the nowadays era of precision cosmology, such a program is still too ambitious to 
give useful constraints on the F(R) derivatives, as we will see later. On the other hand, the 
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cosmographic parameters may also be expressed in terms of the dark energy density and EoS 
parameters so that we can work out what are the current values of F(R) and its derivatives 
giving the same (qo,jo, Sq, /o) of the given dark energy model. To this aim, it is convenient to 
adopt a parameterized expression for the dark energy EoS in order to reduce the dependence 
of the results on an y und erlying theoretical scenario. Following the prescription of the Dark 
Energy Task 



by setting 



35 



orce 



179J, we will use the CPL parameterization for the EoS of dark energy 



w = w + w a (l - a) = w + w a z(l + z) 1 



11.33) 



so that, in a flat universe filled by dust matter and dark energy, the dimensionless Hubble 
parameter E(z) = H/H in Eq. (12.131) reads 



E\z) 



nS?(i + *) 8 



3(1+W0+W a ) 



exp 



3WgZ 

'l + z 



11.34) 



with 



1 — fijjj because of the assumption that the universe is flat. Here and in the 



following, we omit the inferior "DE" of Wo and w a . In order to determine the cosmographic 
parameters for such a model, we avoid integrating H(z) to get a(t) by noting that d/dt = 
— (1 + z)H(z)d/dz. We can use such a relation to evaluate (H, H,d 3 H/dt 3 ,d 4 H/dt A ) and 
then solve Eqs. ( 110.61) — (110.911 . evaluated at z — 0, with respect to the parameters of interest. 
Some algebra finally gives 



go 

Jo 
so 

lo 



l+l(l-^)w , 



) [3wo(l + w ) + w a ] 



n 



(0) 



= 1 + 2^ 



7 — Q 



■U' n 



w + 



Q {0) 



)(0) 



3 - 



3 

o • 



35 

y 

2 



(0) 
M 



213 



7-n 



n 



w a 



(0) 
M 



489 + 9 82 - 21ft 



W a 



67-21fiS ) + - ^23- lift 



)(0) 
M 



W n 



Wr 



+ 2 1 (1 _ n 



47 - 24ft 



™; ! + y ( 1 - ftM 



3 - 2ft 



(11.35) 
(11.36) 

(11.37) 

+ 

(11.38) 



Inserting Eqs. (111.350 — (jll.38p into Eqs. flll.2ip — ( 111.310 . we acquire lengthy expressions 
(which we do not report here) giving the current values of F(R) and its first three deriva- 
tives as a function of (ft^ ,Wo,w a ). It is worth noting that the F(R) model obtained is 
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not dynamically equivalent to the starting CPL one. Indeed, the two models have the same 
cosmographic parameters only today. For instance, the scale factor is the same between the 
two theories only over the time period during which the fifth order Taylor expansion is a 
good approximation of the actual a(t). It is also meaningful to stress that such a procedure 
does not select a unique F(R) model, but rather a class of fourth order theories all sharing 
the same third order Taylor expansion of F(R). 

B. The ACDM case 

With these caveats in mind, it is significant to first examine the ACDM model, which is 
described by setting (wo,w a ) = (—1,0) in the above expressions, and hence we have 



When inserted into the expressions for the F(R) quantities, these relations give the remark- 
able result : 



so that we obviously conclude that the only F(R) theory having exactly the same cosmo- 
graphic parameters as the ACDM model is just F(R) oc R, i.e., general relativity (GR). 
It is important to mention that such a result comes out as a consequence of the values of 
(<?cbjo) m the ACDM model. Indeed, should we have left (sq,Iq) undetermined and only 
fixed (go, jo) to the values in (jll.39p . we should have got the same result in (111.401) . Since 
the ACDM model fits a large set of different data well, we do expect that the actual values 
of (qo,jo,So,lo) do not differ too much from those in ACDM model. Therefore, we plug 
Eqs. pi.2ip ~ pi.3ip into the following expressions : 



q = q$x(l + e q ), j =j A x(l + 5 i ), s = s£x{l + e s ) , l = l*x (1 + e { ) , (11.41) 



with (q^jjo, Sq, Iq) given by Eqs. (I11.39|) and (e q , Ej, e s , e{) quantify the deviations from the 
values in the ACDM model allowed by the data. A numerical estimate of these quantities 



< 




(11.39) 



F(Ro) = R + 2A , F"(Ro) = F'"{R Q ) = 



(11.40) 
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may be obtained, e.g., from a Markov chain analysis, but this is outside our aims. Since we 
are here interested in a theoretical examination, we prefer to study an idealized situation 
where the four quantities above all share the same value e « 1. In such a case, we can easily 
investigate how much the corresponding F(R) deviates from GR by exploring the two ratios 
F"(R )/F(R ) and F"'(R )/ F(R ). Inserting the above expressions for the cosmographic 
parameters into the exact (not reported) formulae for F(Rq), F"(Rq) and F'"(Ro), taking 
those ratios and then expanding to first order in e, we finally find 

64-6fii? te[ 0) 



^ 20 = t?i p ^ o - 7 : — — : — — ; : — T2 x 7^ > I 11 - 42 ) 



3 ( 9fig } + 74 ) fij? - 556 



F'"(R ) fi i oifi M - no) fig } + 



fi£ } + 16 



F(i? c 



3 (9fig } + 74) fig } - 556 



ngn +16 243 « ) 



which are dimensionless quantities and hence more suitable to estimate the order of magni- 

)(0) 



tudes of the different terms. Inserting our fiducial values into fig^ we have 



0.15 x e for fig 5 = 0.041, 
V20 ^ { (11-44) 
-0.12 x e for fig 5 = 0.250, 



)(0) 

r/30 ~ ^ (11.45) 

)(0) 
5 M 



Axe for Q^ 1 = 0.041 



-0.18 x e for fi[? = 0.250 



For values of e up to 0.1, the above relations show that the second and third derivatives 
are at most two orders of magnitude smaller than the zeroth order term F(R ). Actually, 
the values of 7730 for a baryon only model (first row) seems to argue in favor of a larger 
importance of the third order term. However, we have numerically checked that the above 
relations approximates very well the exact expressions up to e ~ 0.1 with an accuracy 
depending on the value of fig , being smaller for smaller matter density parameters. Using 
the exact expressions for 7720 and 7730, our conclusion on the negligible effect of the second 
and third order derivatives are significantly strengthened. 

Such a result holds under the hypotheses that the narrower are the constraints on the 
validity of the ACDM model, the smaller are the deviations of the cosmographic parameters 
from the values in the ACDM model. It is possible to show that this indeed the case for the 
CPL parameterization we are considering. On the other hand, we have also assumed that 
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the deviations (s q , Ej, e s , Si) take the same values. Although such hypothesis is somewhat 
ad hoc, we argue that the main results are not affected by giving it away. In fact, although 
different from each other, we can still assume that all of them are very small so that Taylor 
expanding to the first order should lead to additional terms into Eqs. ( I11.42l) - (lll.43p which 
are likely of the same order of magnitude. We may thus conclude that, if the observations 
confirm that the values of the cosmographic parameters agree within ~ 10% with those 
predicted by the ACDM model, we must recognize that the deviations of F(R) from the GR 
case, F(R) oc R, should be vanishingly small. 

It should be emphasized however, that such a conclusion only holds for those F(R) 
models satisfying the constraint (jll.32p . It is indeed possible to work out a model having 
F{R ) oc #0, F"{Rq) = F"'{R ) = 0, but F^(Ro) ^ for some n. For such a (somewhat ad 
hoc) model, Eq. fl 1 1 .32j) is clearly not satisfied so that the cosmographic parameters have to 
be evaluated from the solution of the field equations. Accordingly, the conclusion above does 
not hold, so that one cannot exclude that the resultant values of (go, jo, so, fo) are within 
10% of those in the ACDM model. 



C. The constant EoS model 



Let us now take into account the condition w = — 1, but still retains w a = 0, thus 
obtaining the so called quiescence models. In such a case, some problems arise because both 
the terms (jo — qo~ 2) and 1Z may vanish for some combinations of the two model parameters 
(fi^ \ w Q ). For instance, we find that j 



Qo 



n 



(o) 

M 



W 2 



1 

3 



1 + 



2 = for w = (wi, w 2 ) with 
1 



i-n 



■off) 



n 



(0) 
M 



4-fi 



(11.46) 



11.47) 



On the other hand, the equation TZ(Q^,wq) = may have different real roots for w de- 
pending on the adopted value of . Denoting collectively with w nu u the values of wo 



that, for a given fi^ , make (jo — qo ~ 2)7£(^m\ Wq) taking the null value, we individuate 



(0) 



a set of quiescence models whose cosmographic parameters give rise to divergent values of 
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FIG. 1: The dimensionless ratio between the current values of F"{R) and F{R) as a function of 
the constant EoS wq of the corresponding quiescence model. Short dashed and solid lines refer to 
models with = 0.041 and 0.250, respectively. 

F(Rq), F"(Rq) and F"'(Rq). For such models, F(R) is clearly not denned, so that we have 
to exclude these cases from further consideration. We only note that it is still possible to 
work out an F(R) theory reproducing the same background dynamics of such models, but 
a different route has to be used. 

Since both go and Jo now deviate from those in the ACDM model it is not surprising 
that both F"(Rq) and F"'(Rq) take finite non null values. However, it is more interesting 
to study the two quantities 7720 and 7730 denned above to investigate the deviations of F(R) 
from the GR case. These are plotted in Figs. [T] and [2] for the two fiducial values of . 
Note that the range of wo in these plots have been chosen in order to avoid divergences, but 
the lessons we will draw also hold for the other wo values. 

As a general comment, it is clear that, even in this case, F"(Rq) and F"'(Rq) are from 
two to three orders of magnitude smaller than the zeroth order term F(Rq). Such a result 
could yet be guessed from the previous discussion for the ACDM case. Actually, relaxing the 
hypothesis wo = — 1 is the same as allowing the values of the cosmographic parameters to 
deviate from those in the ACDM model. Although a direct mapping between the two cases 
cannot be established, it is nonetheless evident that such a relation can be argued and hence 
make the outcome of the above plots not fully surprising. It is nevertheless worth noting 
that, while in the ACDM case, 7720 and 7730 always have opposite signs, this is not the case for 
quiescence models with w > — 1. Indeed, depending on the value of Qyj* , we can have F(R) 
theories with both 7/20 and r/30 positive. Moreover, the lower is, the higher the ratios 7/20 
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FIG. 2: The dimensionless ratio between the current values of F"'(R) and F(R) as a function of 
the constant EoS wq of the corresponding quiescence model. Legend is the same as in Fig. [U 

and rj 30 are for a given value of Wq. This can be explained qualitatively noticing that, for a 
lower the density parameter of the curvature fluid (playing the role of an effective dark 
energy) must be larger and thus claim for higher values of the second and third derivatives. 

D. The general case 

Finally, we study evolving dark energy models with w a ^ 0. Needless to say, varying 
three parameters allows to get a wide range of models that cannot be discussed in detail. 
Therefore, we only concentrate on evolving dark energy models with wq = — 1 in agreement 
with some most recent analysis. The results on 7720 and 7730 are plotted in Figs. [3] and H] where 
these quantities are displayed as functions of w a . Note that we are exploring models with 
positive w a so that w(z) can tend to wo + w a > wo for z — > 00 and the EoS for dark energy 
can eventually approach the dust value w — 0. Actually, this is also the range favored by 
the data. We have, however, excluded values where 7720 or rj 30 diverge. Considering how 
they are defined, it is clear that these two quantities diverge when F(Rq) = 0, so that the 
values of (w ,w a ) making (7720,^30) diverge may be found by solving 



where Vo(w ,w a ) and Qo(w ,w a ) are obtained by inserting Eqs.f lll.35|) - flll.38p into the 
definitions f lll.25p -f lll.26p . For such CPL models, there is no F(R) model having the same 
cosmographic parameters and, at the same time, satisfying all the criteria needed for the 




(11.48) 
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FIG. 3: The dimensionless ratio between the present day values of F"(R) and F(R) as function of 
the w a parameter for models with wq = — 1. Legend is the same as in Fig. [TJ 

validity of our procedure. In fact, if F(R ) = 0, the condition ( 111. 32[) is likely to be violated 
so that higher than third order must be included in the Taylor expansion of F(R) thus 
invalidating the derivation of Eqs. ( II 1 .2Tj) - (ll 1 ,24p . 

Under these caveats, Figs. [3] and H] demonstrate that allowing the EoS for dark energy to 
evolve does not change significantly our conclusions. Indeed, the second and third deriva- 
tives, although being not null, are nevertheless negligible with respect to the zeroth order 
term, and therefore the consequence is in favor of a GR-like F(R) with only very small 
corrections. Such a result is, however, not fully unexpected. From Eqs. (111.351) and (111. 36ft . 
we see that, having set wq = —1, the parameter q is the same as that in the ACDM model, 
while jo reads j$ + (3/2) ^1 — Q^^j w a- As we have stressed above, the Einstein-Hilbert 
Lagrangian F(R) = R + 2A is recovered when (qo,jo) = (SchioO whatever the values of 
(sq, Iq) are. Introducing a w a ^ makes (so, lo) differ from those in the ACDM model, but 
the first two cosmographic parameters are only mildly affected. Such deviations are then 
partially washed out by the complicated way they enter in the determination of the values 
of F(R) at the present time and its first three derivatives. 

XII. THEORETICAL CONSTRAINTS ON THE MODEL PARAMETERS 

In the preceding section, we have worked out a method to estimate F(R ), F"(R ) and 
F"\Rq) resorting to a model independent parameterization of the EoS for dark energy. 
However, in the ideal case, the cosmographic parameters are directly estimated from the data 
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FIG. 4: The dimensionless ratio between the present day values of F"'(R) and F(R) as function 
of the w a parameter for models with wq = — 1. Legend is the same as in Fig. [TJ 

so that Eqs. fll 1.21D — fll 1.31H can be used to infer the values of the quantities related to F(R). 
These latter can then be used to put constraints on the parameters entering an assumed 
fourth order theory assigned by an F(R) function characterized by a set of parameters 
p = (pi, . . . ,p n ) provided that the hypotheses underlying the derivation of Eqs. (Ill .21j) — 
f lll.3ip are indeed satisfied. We show below two interesting cases which clearly highlight the 
potentiality and the limitations of such an analysis. 




A. Double power law Lagrangian 



As a first interesting example, we take 180 ] 



F(R) = R(1 + aR n + /3R~ m ) 



;i2.i) 



with n and m two positive real numbers. The following expressions are immediately ob- 
tained : 



-m 
L ' 



F(R ) = R (1 + aR n Q + /3R^ m ) , 
F'(R ) = 1 + a(n + 1)R% - 0(m - 1)R' C 
F"{Ro) = an{n + l)R^ + (3m{m - l)R~ (1+m) , 
F'"(Ro) = an(n + l)(n - 1)R^ 2 - (3m(m + l)(m - l)R~ {2+m) 



(12-2) 
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Denoting the values of F^\R ) determined through Eqs. flll.2ip - flll.3ip by <pi (with i 
0, . . . , 3), we can solve 



F(Ro) - 
F'(Ro) 
F"(R ) 
F"'(Ro) 



112-3) 



'3 • 



which is a system of four equations in the four unknowns (a, (3,n,m) that can analytically 
be solved proceeding as follows. First, we solve the first and second equation with respect 
to (a, (3) and we obtain 

(12.4) 



a 







n + m ' 
1 + n 



1 - 



R 



) 



n + m \ Rq 
while, solving the third and fourth equations, we get 

2 Rl~n [ 1+m+ (0 3 /0 2 ) jRo ] 



a 



n(n + l)(n + m) 

^ 2 Rl +n [1-71+ (fo/<p2)R ] 



;i2.5) 



m(l — m){n + m) 

Equating the two solutions, we acquire a system of two equations in the two unknowns 
(n, m), given by 

n(n + l)(l-m)(l-0oARo) = 1 

(12.6) 



cj) 2 Ro [l + m + {fal<h)Ro\ 
m (n + l)(m-l)(l-0o/i2o) 



1 . 



(j) 2 R [1-71+ (0 3 /0 2 )i?o] 
Solving with respect to m, we find two solutions. The first one is m = —n, which has to be 
discarded because it makes (a, (3) go to infinity. The only acceptable solution is 



m 



[1-n 



5 /0 2 )#o] 



;i2.7) 



which, inserted back into the above system, leads to a second order polynomial equation for 
n with solutions 



n 



, , 03 D , a/A/"(0 O ,02,03) 
1 + —Kn±- — — — ; — — — 



IRo) 



12.8) 
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Here, we have defined 

A/"(0O, 02, 3 ) = {Rl<Pl - 2i?o0O + R 4 ) 03 + 6 (#000 ~ + R 3 ) 0203 + 

+ 9 (<j>l - 2^00 + R 2 ) 0^ + 4 (R 2 (f> - R 3 ) <f>\ . (12.9) 

Depending on the values of (go, Jo, s o, ^o), Eq. (112.81) may lead to one, two or any acceptable 
solution, i.e. real positive values of n. This solution has to be inserted back into Eq. (112. 7p 
to determine m and then into Eqs. (112. 4j) or ( 112. 5p to estimate (a,/3). If the final values 
of (a, P,n,m) are physically viable, we can conclude that the model in Eq. (112. ip is in 
agreement with the data giving the same cosmographic parameters inferred from the data 
themselves. Exploring analytically what is the region of parameter space of (go, Jo, s o, ^o) 
which leads to acceptable solutions of (a, (3,n,m) is a daunting task far outside the aim of 
the present work. 



B. The Hu-Sawicki model 



One of the most pressing problems of F(R) theories is the need to escape the severe con- 
straints imposed by the Solar System tests. A successful model has recently been proposed 



by Hu and Sawicki 



i8i| (nsy 



F(R) =R-R C 



a(R/R c 



12.10) 



1 + f3(R/R c ) n 

As for the double power law model discussed above, there are four parameters which we can 
be expressed in terms of the cosmographic parameters (go, jo, $o, lo)- 
As a first step, it is trivial to have 



F(R ) - 
F'(Ro) 
F"(Ro) 
F"'(Ro) 



Ro — R. 



aR% c 

i + mt 

anR c R,Q C 



anR c R% c [(1 - n) + (3(1 + n)R% c } 



12.11) 



Rl(l 



(3R% ^ 



L 0c) 



anR c R^ c (An 2 + Bn + C) 

rW+prW 



4 Note that such a model does not pass the matter instability test and therefore some viable generaliza- 
have been proposed. 



tions 
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with R 0c = Rq/R c and 



A = -f3 2 Rl n c + 

B — 3(1- 2 R%) , 



12.12) 



denned as above, we could, in 



C = -2(1 -(3RD 2 . 

Equating Eqs. (112. lip to the four quantities (0o, 0i, 02, 03) 
principle, solve this system of four equations in four unknowns to obtain (a, (3, R c ,n) in 
terms of (0o, 0i, 02, 03)- Then, by using Eqs. f lll.2ip - flll.3ip we acquire the expressions 
of (a, fl, R c ,n) as functions of the cosmographic parameters. However, setting 0i = 1 as 
required by Eq. f II 1 . 2 2 j) gives the only trivial solution anR c = so that the HS model 
reduces to the Einstein - Hilbert Lagrangian F(R) = R. In order to escape this problem, 
we can relax the condition F'(Rq) = 1 to F'(Rq) = (1 + e) -1 . As we have discussed 
in Sec. XI, this is the same as assuming that the current effective gravitational constant 
G e s,o = G/F'(R ) only slightly differs from the usual Newton's one, which seems to be a 
quite reasonable assumption. Under this hypothesis, we can analytically solve the equations 
for (a, f3, R c ,n) in terms of (0 O , e, 02, 03). The actual values of (0o,02,03) will be no more 
given by Eqs. flll.2ip -f lll.24p . but we have checked that they deviate from those expressions 5 
much less than 10% for e up to 10% well below any realistic expectation. 
With this caveat in mind, we first solve 



F(R 



F"(R ) 



12.13) 



to acquire 



a 



P 



n(l +e) 

e 

n(l + e) 

e 



Rq 
R c 

Rc 



1-n 



<Po_ 
R 



Rq 



nil + £ 



Inserting these expressions into the equations in (112. lip , it is easy to check that R c cancels 
out, so that we can no more determine its value. Such a result is, however, not unexpected. 
Indeed, Eq. (112. 10p can trivially be rewritten as 

aR n 



F(R) = R- 



l + /3R r ' 



12.14) 



Note that the correct expressions for (tfio, <p2, fo) may still formally be written as Eqs. (|11.21|) — (|1 1 .24[) . 
but the polynomials entering them are now different and also depend on powers of s. 
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with a = aR]r n and $ = j3R~ n which are indeed the quantities determined by the above 
expressions for (a,j3). Reversing the discussion, the current values of F^ l \R) depend on 
(a,P,R c ) only through the two parameters (&,$). Accordingly, the use of cosmographic 
parameters is unable to break this degeneracy. However, since R c only plays the role of a 
scaling parameter, we can arbitrarily set its value without loss of generality 

On the other hand, this degeneracy allows us to get a consistency relation to immediately 
check whether the HS model is viable or not. Indeed, solving the equation F"(Rq) = 02, we 
find 

_ (0p/go) + [(1 + £)/e] (1 - 2 i? o ) - (1 - e )/(l + e) 

1 - 0o/ 

which can then be inserted into the equations F"'(Rq) = 03 to obtain a complicated relation 
among (0o, 02, 03) which we do not report for sake of shortness. Solving such a relation 
with respect to 03/0o and executing the Taylor expansion to first order in e, the resultant 
constraint reads 



^3 



1 + 6 



e Rn 



edn 1 / 2e 



12.16) 



1 + e \ 1 - o /i?o 

If the cosmographic parameters (qo, jo, $o, lo) are known with sufficient accuracy, one could 
compute the values of (i?o, 0o, 02-03) for a given e (eventually, using the expressions obtained 
for e — 0) and then check if they satisfied this relation. If this is not the case, one can 
immediately give off the HS model also without the need of solving the field equations and 
fitting the data. In fact, given the still large errors on the cosmographic parameters, such a 
test only remains in the realm of ( quite distant) future applications. However, the HS model 
works for other tests as shown in 1811 ] and thus a consistent cosmography analysis has to 
be combined with them. 



XIII. CONSTRAINTS COMING FROM OBSERVATIONAL DATA 

Eqs. (I11.21l) - (lll.3ip relate the values of F(R) and its first three derivatives at the present 
time to the cosmographic parameters (go, jo, s o, ^o) and the matter density £1$ • ^ n principle, 
therefore, a measurement of these latter quantities makes it possible to put constraints on 
F^(Rq), with % = {0, . . . , 3}, and hence on the parameters of a given fourth order theory 
through the method shown in the previous section. Actually, the cosmographic parameters 
are affected by errors which obviously propagate onto the F(R) quantities. Indeed, since 
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the covariance matrix for the cosmographic parameters is not diagonal, one has to also take 
care of this fact to estimate the final errors on F®(Rq). A similar discussion also holds 
for the errors on the dimensionless ratios 7/20 and 7730 introduced above. As a general rule, 
indicating with g(Q^,p) a generic F(R) related quantity depending on fij^ and the set of 
cosmographic parameters p, its uncertainty reads 



dg 



i=4 



a 



M 



i=l 



dpi p ' ^ dpidpi 3 



2 



dg dg 



where are the elements of the covariance matrix (being Cu = 0"^.), we have set 
{Pi,P2,P3,Pa) = (Qo, jo, so, lo)- and assumed that the error om on is uncorrelated with 
those on p. Note that this latter assumption strictly holds if the matter density parame- 
ter is estimated from an astrophysical method (such as estimating the total matter in the 
universe from the evaluated halo mass function). Alternatively, we will assume that tt^ is 
constrained by the CMB radiation related to experiments. Since these latter mainly probes 
the very high redshift universe (z ~ z\ ss ~ 1089), while the cosmographic parameters are 
concerned with the cosmos at the present time, one can argue that the determination of 
fijjj is not affected by the details of the model adopted for describing the late universe. In- 
deed, we can reasonably assume that, whatever the candidate for dark energy or alternative 
gravitational theory such as F(R) gravity is, the decoupling epoch, i.e., the era when we 
can observe through the CMB radiation, is well approximated by the standard GR with a 
model comprising only dust matter. Hence, we will make the simplifying (but well moti- 
vated) assumption that &m may be reduced to very small values and is uncorrelated with 
the cosmographic parameters. 

Under this assumption, the problem of estimating the errors on g(Q^ , p) reduces to the 
analysis of the covariance matrix for the cosmographic parameters given the details of the 
data set used as observational con straints. We address this issue by computing the Fisher 
information matrix (see, e.g., 1861 ] and references therein) defined as 

F - = (md (13 ' 2 > 

with L = —2 In C(9i, . . . , 8 n ), where C(9i, . . . , 9 n ) is the likelihood of the experiment and 
(#1, . . . ,9 n ) is the set of parameters to be constrained, and (. . .) denotes the expectation 
value. In fact, the expectation value is computed by evaluating the Fisher matrix elements 
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for fiducial values of the model parameters (6\, . . . ,9 n ), while the covariance matrix C is 
finally obtained as the inverse of F. 

A key ingredient in the computation of F is the definition of the likelihood which depends, 
of course, of what experimental constraint one is using. To this aim, it is worth remembering 
that our analysis is based on fifth order Taylor expansion of the scale factor a(t) so that 
we can only rely on observational tests probing quantities that are well described by this 
truncated series. Moreover, since we do not assume any particular model, we can only 
characterize the background evolution of the Universe, but not its dynamics which, being 
related to the evolution of perturbations, unavoidably need the specification of a physical 
model. As a result, the SNela Hubble diagram is the ideal test to constrain the cosmographic 
parameters. We therefore defined the likelihood as 



2( TJ \\ A/sNoIa r I . \ I .. TT \ 1 2 



C(H ,p) oc exp ^ y — J , x 2 (#o,p) = Yl 

where the distance modulus to redshift z is described as 



/^obs(^) - Hih{Zn, H , p) 



;i3.3) 



/i th (z, H , p) = 25 + 5 log 10 l-^-j + 5 log 10 D L (z, p) . (13.4) 

Here, Dl(z) is the Hubble free luminosity distance and it follows from Eq. (12 .9p that Dl(z) 
in the flat universe is expressed as 

Using the fifth order Taylor expansion of the scale factor, we find an analytical expression for 
Dl(z, p), so that in the computation of no numerical integration need (this consequence 
makes the estimate faster). As a last ingredient, we need to specify the details of the SNela 
survey giv ing t he redshift distribution of the sample and the error on each measurement. 
Following 183], we adopt 6 




(13-6) 

where z max is the maximum redshift of the survey, cr sys is an irreducible scatter in the SNela 
distance modulus, and a m is to be assigned depending on the photometric accuracy. 



6 Note that, in |l87l |. the authors assume the data are separated in redshift bins so that the error becomes 
a 2 = <J 2 ys /Nbin + Nbin{z/z ma x) w hh Nun the number of SNela in a bin. However, we prefer to not 
bin the data so that Nu n = 1. 
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In order to run the Fisher matrix calculation, we have to set a fiducial model, according 
to the predictions in the ACDM model for the cosmographic parameters. For Q$ = 0.3 and 



the Hubble constant h = 0.72 



68 



in units of lOOkm/s/Mpc, we acquire 



(<?o, Jo, s , l ) = (-0.55,1.0,-0.35,3.11) . (13.7) 
As a first consiste ncy check, we compute the Fisher matrix for a survey mimicking the 



recent database in [188j and thus set (AsNeia, &m) = (192,0.33). After marginalizing over h 
(which, as well known, is fully degenerate with the SNela absolute magnitude Ai), we find 
the uncertainties 

((7i, a 2 , a 3 , <r 4 ) = (0.38, 5.4, 28.1, 74.0) , (13.8) 

where we are still using the indexing introduced above for the cosmographic parameters. 
These values are reasonabl y wel l com pared with those obtained from a cosmographic fitting 



of the Gold SNela dataset 



189 



190j] 7 



q = -0.90±0.65, j = 2.7±6.7, s = 36.5±52.9 , l = 142.7±320 . (13.9) 

Because of the Gaussian assumptions we rely on, the Fisher matrix forecasts are known to 
be lower limits on the determination of a set of parameters, the accuracy to which a given 
experiment can attain. This is indeed the case with the comparison suggesting tha t our 
predictions are quite optimistic. It is worth stressing, however, that the analysis in 1891 . 190 1 



used the Gold SNela dataset which is poorer in high redshift SNela than the one in [188 ] 
we are mimicking, so that larger errors on the higher order parameters (sq, Iq) could be 
expected. 

Rather than computing the errors on F(Rq) and its first three derivatives, it is more 
interesting to look at the precision attainable on the dimensionless ratios (7720, Vso) introduced 
above because they quantify how much deviations from the linear order exist. For the fiducial 
model we are considering, both 7720 and 7730 vanish, while, using the covariance matrix for a 
survey at the present time and setting (Tm/^m — 10%, their uncertainties read 

(^20,^30) = (0.04,0.04) . (13.10) 



7 Actually, such estimates have been obtained by computing the mean and the standard deviation from 
the marginalized likelihoods of the cosmographic parameters. Hence, the central values do not represent 
exactly the best fit model, while the standard deviations do not give a rigorous description of the error 
because the marginalized likelihoods are manifestly non-Gaussian. Nevertheless, we arc mainly interested 
in an order of magnitude estimate so that we would not care about such statistical details. 
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As an application, we can look at Figs. [T] and |2] showing how (7720, ^30) depend on the current 
value of the EoS wq for F(R) models sharing the same cosmographic parameters of a dark 
energy model with its constant EoS. As it is clear, also considering only the la range, the 
full region plotted is allowed by such large constraints on {r]2o,V3o)- Thus, this means that 
the full class of corresponding F(R) theories is viable. As a consequence, we may conclude 
that the SNela data at the present time are unable to discriminate between a A dominated 
universe and this class of fourth order gravity theories. 



As a next step, we investigate a SNAP - like survey [191| and the refor e take (AsNeia, & m ) — 
(2000,0.02). We use the same redshift distribution in Table 1 of 1871 ] and add 300 nearby 
SNela in the redshift range (0.03, 0.08). The Fisher matrix calculation gives the uncertainties 
on the cosmographic parameters 

(<7i, a 2 , a 3 , a A ) = (0.08, 1.0, 4.8, 13.7) . (13.11) 

The significant improvement of the accuracy in the determination of (go, jo, sq, Iq) is trans- 
lated into a reduction of the errors on (7720, 7730)5 which is now given by 

(^20, ^30) = (0.007, 0.008), (13.12) 

where we have supposed that, when SNAP data will be available, the matter density param- 
eter would be determined with a precision a M /Oy ~ 1%. Looking again at Figs. [1] and 
[21 it is clear that the situation is improved. Indeed, the constraints on 7720 makes it possible 
to narrow the range of allowed models with low matter content (the dashed line), while 
models with typical values of Q$ are still viable for wq covering almost the full horizontal 
axis. On the other hand, the constraint on r/ 30 is still too weak, so that almost the full region 
plotted can be allowed. 

Finally, we examine an hypothetical future SNela survey working at the same photomet- 
ric accuracy as SNAP and with the same redshift distribution, but i ncreasing the nu mbe r 
of SNela up to A/Wia = 6xl0 4 as expected from, e.g., DES [192j . PanSTARRS [19 j |. 
SKYM APP ER 194J, while still larger numbers may potentially be observed by AL- 
PACA 



195| and LSST [196| . Such a survey can achieve 

(<7i, (7 2 , (J 3 , <r 4 ) = (0.02,0.2,0.9,2.7) 



so that, with o"m/^ 



(o) 

M 



0.1%, we obtain 

(^20,^30) = (0.0015,0.0016) 



13.13) 



13.14) 
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Fig. Q] shows that, with such a precision on 7720, the region of wo values allowed essentially 
reduces to the value in the ACDM model, while, from Fig. [21 it is clear that the constraint 
on 7730 definitively excludes models with low matter content further reducing the range of 
values of Wq to quite small deviations from the Wq = — 1. We can therefore conclude that 
such a survey will be able to discriminate between the concordance ACDM model and all the 
F(R) theories giving the same cosmographic parameters as quiescence models other than 
the ACDM itself. 

A similar discussion may be repeated for F(R) models sharing the same values of 
(Qo, jo, s oJo) as those in the CPL model even if it is less intuitive to grasp the efficacy 
of the survey being the parameter space multivalued. For the same reason, we have not 
explored what is the accuracy on the double power -law or HS models, even if this is tech- 
nically possible. In fact, one should first estimate the errors on the current values of F(R) 
and its three time derivatives and then propagate them on the model parameters by using 
the expressions obtained above. 

In conclusion, notwithstanding the common claim that we live in the era of precision 
cosmology, the constraints on (qo,jo, so, lo) are still too weak to efficiently apply the program 
we have outlined above. We have shown how it is possible to establish a link between the 
popular CPL parameterization of the EoS for dark energy and the derivatives of F(R), 
imposing that they share the same values of the cosmographic parameters. This analysis 
has led to the quite interesting conclusion that the only F(R) function, which is able to give 
the same values of (qo,jo,so,lo) as those in the ACDM mode l, is indeed F(R) = R + 2A. 
A similar conclusion holds also in the case of f(T) gravity 147] . If future observations 
will inform us that the cosmographic parameters are those of the ACDM model, we can 
therefore rule out all F(R) theories satisfying the hypotheses underlying our derivation of 
Eqs. (jll.2ip -( lll.24p . Actually, such a result should not be considered as a no way out for 
higher order gravity. Indeed, one could still work out a model with null values of F"(Rq) and 
F'"(Rq) as required by the above constraints, but non - vanishing higher order derivatives. 
One could well argue that such a contrived model could be rejected on the basis of the 
Occam razor, but nothing prevents from still taking it into account if it turns out to be both 
in agreement with the data and theoretically well founded. 

If new SNela surveys will determine the cosmographic parameters with good accuracy, 
acceptable constraints on the two dimensionless ratios 7720 oc F"(R )/F(R ) and 7730 cx 
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F"'(Rq)/F(Rq) could be obtained, and thus these quantities allow us to discriminate among 
rival F(R) theories. To investigate whether such a program is feasible, we have pursued a 
Fisher matrix based forecasts of the accuracy, which future SNela surveys can achieve, on 
the cosmographic parameters and hence on (7/20 > ^30) ■ It turns out that a SNAP -like survey 
can start giving interesting (yet still weak) constraints allowing us to reject F(R) models 
with low matter content, while a definitive improvement is achievable with future SNela 
survey observing ~ 10 4 objects and hence makes it possible to discriminate between the 
ACDM model and a large class of fourth order theories. It is worth emphasizing, however, 
that the measurement of should come out as the result of a model independent probe 
such as the gas mass fraction in galaxy clusters which is, at present, still far from the 1% 
requested precision. On the other hand, one can also rely on the estimate from the 
anisotropy and polarization spectra of the CMB radiation even if this comes to the price of 
assuming that the physics at recombination is strictly described by GR, so that one has to 
limit its attention to F(R) models reducing to F(R) cx R during that epoch. However, such 
an assumption is quite common in many F(R) models available in literature and therefore 
it is not a too restrictive limitation. 

A further remark is in order concerning what kind of data can be used to constrain the 
cosmographic parameters. The use of the fifth order Taylor expansion of the scale factor 
makes it possible to not specify any underlying physical model by relying on the minimalist 
assumption that the universe is described by the flat FLRW metric. While useful from a 
theoretical perspective, such a generality puts severe limitations to the dataset one can use. 
Actually, we can only resort to observational tests depending only on the background evo- 
lution so that the ra nge o f astrophysical probes reduces to standard candles (such as SNela 



and possibly GRBs 



1671 ]) and standard rods (such as the angular size - redshift relation 



for compact radiosources) . Moreover, pushing the Hubble diagram to z ~ 2 may rise the 
question of the impact of gravitational lensing amplification on the apparent magnitude of 
the adopted standard candle. T he magni fication probability distribution function depends 
on the growth of perturbations jl97 201 1 . so that one should worry about the underlying 
physical model in order to estimate whether this effect bia ses the estimate of the cosmo- 



graphic parameters. However, it has been shown |202H206| that the gravitational lensing 
amplification does not alter the measured distance modulus for z ~ 1 SNela significantly. 
Although such an analysis has been executed for models based on GR, we can argue that, 
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whatever the F(R) model is, the growth of perturbations finally leads to a distribution of 
structures along the line of sight that is as similar as possible to the observed one so that 
the lensing amplification can be approximately the same. We can therefore discuss that the 
systematic error made by neglecting lensing magnification is lower than the statistical ones 
expected by the future SNela surveys. On the other hand, one can also try reducing this 
possible bias further by using the method of flux averaging 2071 ] even if, in such a case, our 
Fisher matrix calculation should be repeated accordingly. Furthermore, it is significant to 
note that the constraints on the cosmographic parameters may be tightened by imposing 
some physically motivated priors in the parameter space. For instance, we can suppose that 
the Hubble parameter H(z) always stays positive over the full range probed by the data 
or that the transition from past deceleration to present acceleration takes place over the 
range probed by the data (so that we can detect it). Such priors should be included in 
the likelihood definition so that the Fisher matrix should be recomputed. This is left for a 
forthcoming work. 

Although the data at the present time are still too limited to efficiently discriminate 
among rival dark energy models, we are confident that an aggressive strategy aiming at a 
very precise determination of the cosmographic parameters could offer stringent constraints 
on higher order gravity without the need of solving the field equations or addressing the 
complicated problems related to the growth of perturbations. Almost 80 years after the 
pioneering distance- redshift diagram by Hubble, the old cosmographic approach appears 
nowadays as a precious observational tool to investigate the new developments of cosmology. 



XIV. CONCLUSION 



In summary, we have presented the review of a number of popular dark energy mod- 
els, such as the ACDM model, Little Rip and Pseudo-Rip scenarios, the phantom and 
quintessence cosmologies with the four types (I, II, III and IV) of the finite-time future 
singularities and non-singular universes filled with dark energy. 

In the first part, we have explained the ACDM model and recent various cosmo logical 
observations to give the bounds on the late-time acceleration of the universe. Furthermore, 
we have investigated a fluid description of the universe in which the dark fluid has a general 
form of the EoS covering the inhomogeneous and imperfect EoS. We have explicitly shown 
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that all the dark energy cosmologies can be realized by different fluids and also considered 
their properties. It has also been demonstrated that at the current stage the cosmological 
evolutions of all the dark energy universes may be similar to that of the ACDM model, and 
hence these models are compatible with the cosmological observations. In particular, we 
have intensively studied the equivalence of different dark energy models. We have described 
single and multiple scalar field theories, tachyon scalar theory and holographic dark energy, 
in which the quintessence/phantom cosmology with the current cosmic acceleration can be 
represented, and eventually verified those equivalence to the corresponding fluid descriptions. 

In the second part, as another equivalent class of dark energy models, in which dark 
energy has its geometrical origins, namely, modifications of gravitational theories, we have 
examined F(R) gravity including its extension to F(R) Hofava-Lifshitz gravity and f(T) 
gravity. It has clearly been explored that in these models, the ACDM model or the late-time 
cosmic acceleration with the quintessence/phantom behavior can be performed. 

Finally, it is significant to remark that there are a number of var ious dark energy 
models which we did not discuss in this review, such as F(Q) gravity 2081 ] . where Q = 
R 2 — AR^ V R^ U + R fl v pcr R fJ '' ypcr with R^ v and R^ v ^ a being the Ricci tensor Riemann tensors, re- 
spective ly, i s the Gauss-Bonnet invariant, F(R,Q) gravity 2091 ], scalar- Gauss- Bonnet dark 



energy 



2101 ]. k-essence da rk en ergy models 94|, ghost condensates scenario 



extension to inflation, see 



2111 ] (for its 



2121 ]). viscous dark energy 8J, [85] , non- minimal derivative dark 



energy models 2131 ] . G-essence dark energy models [214j . non-local gravity 56( produce by 
quantum effect, which is investigated to account for t he c oincidence problem of dark en- 
ergy and dark matter, and galileon dark energy models 215| (for its application to inflation, 



called G-inflation , wh ich has recently been proposed, see 



models, see, e.g. 



216j ) [as recent reviews on galileon 



2171 ]] . In particular, galileon gravity has recently been studied very exten- 



sively in the literature. The most important feature of the Lagrangian for the galileon scalar 
field is that the equation of motion derived from the Lagrangian is up to the second-order, 
so that the appearance of an extra degree of freedom with the existence of a ghost can be 
avoided. The galileon field originates from a brane bending mode in the Dvali-Gabadadze- 
Porrati (DGP) brane world scenario 2181 ]. and therefore galileon gravity might be regarded 
as an indirect resolution for the issue of a ghost in the self-accelerating branch of the DGP 
model. Since we have no enough space to describe the details of all these models, we again 
mention the important procedure of our approach to show the equivalence of dark energy 
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models to represent each cosmology. In all of the above models, it follows from Eqs. (I3.3p 
and (13. 4p that in the flat FLRW background the gravitational equations can be described as 
H 2 = (k 2 /3) pde and H = — (k 2 /2) (p DE + Pde)- In each model, the difference is only the 
forms of the energy density pde and pressure Pde of dark energy. Hence, the expression of 
the Hubble parameter H to describe the concrete cosmology, e.g., the ACDM, quintessence 
and phantom cosmologies, can be reconstructed by using these gravitational field equations. 
Similarly, by applying pde and Pde, the EoS wde = DE/pde in the fluid description in 
Eq. (13.11) with Eq. (I3.2p can also be presented. 

Finally, it is worth stressing the role of cosmography in this discussion. As shown, it is 
a fundamental tool because it allows, in principle, to discriminate among models without 
a priori assumptions but just laying on constraints coming from data. However, the main 
criticism to this approach is related to the extension of the Hubble series, the quality and the 
richness of data samples. In particular, observations cannot be extended at any redshift and, 
in most of cases, are not suitable to track models up to early epochs. However, the forth- 
coming observational campaigns should ameliorate the situation removing the degeneration 
emerging at low redshifts and allowing a deeper insight of models. 
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Appendix A: Inertial Force and u>de 

In this appendix, we check if the occurrence of inertial force in the universe with rip 
(Little rip or the finite-time future singularities) may somehow constrain the wde- 

As the universe expands, the relative acceleration between two points separated by a 
distance I is given by Id/ a. If there is a particle with mass m at each of the points, an 
observer at one of the masses will measure an inertial force on the other mass of 

-Fincrt = mla /a = ml (h + H 2 ^j . (Al) 

By using the deceleration parameter g dcc = —aa~ 1 H~ 2 in Eq. (I3.97p . we may express the 
inertial force -F inert as 

^inert = -mlH 2 q dec . (A2) 
The observational constraint of the value gdec = <7dec(o) i n the present universe is given by j^O | 

- 0.60 < g dcc(0) < -0.30 . (A3) 

Here, we mean the present value by the suffix 0. The present value of the Hubble rate 
H = H could be 

H ~ 70 km/ (s ■ Mpc) = 2.3 x 10^ 18 s" 1 . (A4) 

We now rewrite the expression in flA2j) by using the gravitational field equations in the 
FLRW space-time, 



2 H 2 = p, -—(2H + 3H 2 )=P J (A5) 



3 ^ - n 1 

2 

K z K 

in Eqs. (I3.29P and f)3.30p . Since the pressures of the usual matters and cold dark matter are 
negligible, if we neglect the contribution from the radiation in the present universe, we may 
identify the pressure P with that of the dark energy: P = Pde = ^dePde- Thus, we obtain 

^de^de = ~ (2H + 3# 2 ) , (A6) 

where £7de = Pde/ Perl w ith p„i t = (3/ k2 )Hq in the first relation in Eq. (12. 6p . Then, the 
inertial force -Fi ne rt has the following form: 

TfllH 2 

-^inert tt 1 C 1 + 3w D e^de) • (A7) 
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The galaxy group has a size of 10 23 m and the galaxy is moving with the speed of 
10 5 m/s. Hence, the acceleration by the central force can be estimated to be a ccn trai ~ 
((10 5 ) 2 /i0 23 ) m/s 2 = 10" 13 m/s 2 . On the other hand, IH% ~ 10 23 x (1(T 18 ) 2 = l(T 13 m/s 2 . 
Therefore, we find a centra i ~ Wq, which informs that the precise measurement of the sizes 
of galaxy groups and/or galaxy clusters and the rotation speeds of galaxies may give a 
constraint on the value of Wde^de- 

Conversely, provided that wbe ~ — 1, ^de ~ 0.73, Eq. (IA7I) implies 

F inert ~ 0.6m/# 2 ~ (3.2~ 36 s" 2 ) x ml . (A8) 

On the other hand, if we choose I as the size of a galaxy group or galaxy cluster and let v a 
rotational speed of a galaxy in the galaxy group or galaxy cluster, the central force by the 
gravity -F gravity minus the inertial force could be estimated as 

,2 



mv 



cent J 



(A9) 



Thus, since F cent = Fg ravit y — F{ neTt , by combining Eqs. ( 1A8I) and ( 1A9I) . we acquire 

Fgpoty _^ ^3_ 2 -36 s -2 (A1Q) 

lm r 



If F ccnt ~ -gravity with a difference by a factor, we may find 

^~10- 36 s" 2 . (All) 
I 2 

We may also consider the constraint coming from the energy density of galaxy clusters. 
By using the gravitational field equations (]A5j) . the inertial force ( 1A1I) can be rewritten as 
-Flnert = — (m/fi; 2 /6) (p + 3P) in the first equality in Eq. f)3.126p . Now, we assume the energy 
density p c iuster in a galaxy cluster is almost homogeneous. In this case, the total mass inside 
the sphere with a radius / whose center is the center of cluster is given by 

4-7T 

M = —l 3 p chlster . (A12) 
Accordingly, the Newton gravity which the point particle with mass m suffers is given by 



,mM 4nG , mln 1 
T 2 



Fgrav. — G — — — — - — rnlp c \ nstev — — - — Pciuster • (^13) 
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Therefore, if -Fi ncrt > -F grav ., that is, if — (p + 3P) > p c iustcr, the point particle is separated 
from the cluster. We here define wq = Pj p. Since p c iustcr = 200p, we find the bound for the 
EoS parameter w as w > —67. As a result, this is not so strong constraint. 
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